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1. INTRODUCTION 



The crucial steps in the evolution of any scientific discipline are the identification of the underlying degrees of 
freedom and the dynamics governing them. For the theory of strongly interacting particles these degrees of freedom 
are the quarks and gluons, and the elegant quantum chromodynamics (QCD) Lagrangian 



Cqcd = J2 ^'(-^ + "^»)^' + tr(F;,Aa)2 (1.1) 

i=l 

prescribes the dynamics. 

Here we would like to review how one can deduce directly from >Cqcd , and from its Hamiltonian counterpart (with 
possible additional assumptions), various inequalities between hadronic masses and/or other hadronic matrix elements 
(observables). 

The euclidean correlation functions of color singlet (gauge invariant) local operators Oq. (x) are given by the func- 
tional path integral llir 



d[A,] j d[i,] /d[^]0„,(xi)...Oa,.(x„)e-/'^'^^(^\ (1.2) 



with / A[A^] J d['0] / d['0] indicating the functional integral over the ordinary gauge field and fermionic (Grassman) 
degrees of freedom. By analytically continuing the corresponding momentum space correlations Wai...o„(pi ■ • -Pn) all 
hadronic scattering amplitudes can be determined. 

The simplest two-point functions are particularly useful. The spectral representation for such functions 

Waix,y) = {0\Ja{x)Jl{y)\0) = f d{^^')aa{^i')e-'^\^-y\ (1.3) 



yields information on the hadronic states in the channel with Jq quantum numbers, i.e. the energy- momentum 
eigenstates \n) with non-vanishing (0|Ja|"-) matrix elements. Thus a lowest state of mass mi'^^ implies an asymptotic 
behavior which, up to powers of |x — y|, is 

{0\Mx)Ji{ym e-™^"'!--''! . (1.4) 

The hadronic spectrum can also be directly obtained via the Schrodinger equation 

i^QCD* = , (1.5) 

with 4* a wave functional describing the degrees of freedom of the vale nce quarks and any number of additional qq 
pairs and/or gluons. The complexity of the physical states in Eq. (|l.5|) or the richness of field configurations in the 



functional integral equation (1.2) impede quantitative computations of hadronic matrix elements and the hadronic 
spectrum, a goal pursued over more than two decades, utilizing in particular lattice calculations |^-^. 

The QCD inequalities are derived by comparing expressions for different correlation functions (or the energies 
of different hadronic systems) without requiring explicit evaluation. We only need to assume that an appropriate 
regularization scheme and gauge fixing have been devised to make the path integral (or the Schrodinger problem) well 
defined. 

A key ingredient in deriving relations between correlation functions is the positivity of the functional path integration 
measure d^{A) obtained after integrating out the fermionic degrees of freedom. The bilinear X]i=fi + "ii)V'i part 
of >Cqcd then yields the determinantal factor 

Det = J|Det(^-|-m,) , (1.6) 

which for any vectorial (non-chiral) theory can be shown to be positive for any A^{x) (see Sec. ||) ]6|-p^. 

If the integrand in the path integral expression for one correlation function is greater than the integrand in another 
correlation function for all A^(x), then the positivity of the path integration measure guarantees that this feature per- 
sists for the integrated values. A rigorous inequality between the two correlation functions for all possible (euclidean) 
locations of the external currents will then follow. 
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For the particular case of two-point functions an inequality of the form 

{0\Ux)Jl{y)\0) > {0\Mx)jl{y)\0) 



(1.7) 



implies, via Eq. ( |l.4| ), the reversed inequality for the lowest mass physical states with the quantum numbers of 
Jaix), Jb{x): 



(1.8) 



Most of the inequalities (1.7) involve the pseudoscalar currents {Ja = ipil5'^j) s-nd the corresponding mass inequalities 
(1.8) the pion {i.e. the lowest pseudoscalar states in the U'^sd, u'^^u — d'^^d, d'^^u channels) 0: 



m'°''(any meson) > m^, 



(1.9a) 



m^'^)(any baryon) > m^, 



(1.9b) 



7Ti7r+ > rriTrO . 



(1.9c) 



The efforts to obtain inequalities in the Hamiltoni an a pproach follow a similar general pattern. Rather than 
attempting to solve the QCD Schrodingcr equation (1^) for a particular channel, relations are sought between 
baryonic and/or mesonic sectors of different flavors Bijk, Af^, and different spins. 

Flavor enters the Lagrangian (1.1) only via the bilinear, local mass term. Comparison of masses (or other features) 
of mesons or baryons differing just by flavor may be easier than ab initio computations. The additive form of the 
mass term implies a relationship between the Hamiltonians obtained by restricting the full Hqqb to different flavor 
sectors. Using a variational principle for the ground state masses and assuming flavor symmetric ground state wave 
functions, these relations imply the inequalities 



.(0) 



> 



(0) , (0) 



(1.10) 



with 



,(0) 



the ground state mass in the mesonic sector 



Physical color singlet states can be achieved in a variety of ways: the mesonic - ip^ipja] baryonic - eabc4'i''Pj''Pk'^ 
exotic - e'^^eabc4'ii'j4'k'^i (or hybrid - 'ipaipbGrX^^"'''^) configurations; or glueballs - G^Gr,frst (or drst) G^C^G* [Here 
Gr ~ Gfiu represents the chromo-electromagnetic antisymmetric field tensor. We use abc {rst) for triplet (octet) color 
indices; A''*^"'') are the Gell-Mann matrices; and ijkl refer to flavors. We reserve a/37 for spinor indices, and A/iK 
for vector Lorentz indices.] These configurations (mesonic, baryonic, etc.) differ dynamically by having different 
"color networks". However, the different (non-glueball or hybrid) sectors contain quarks and/or antiquarks which are 
sources or sinks of chromoelectric flux of the same universal strength. This suggests that the QCD Hamiltonians in 
the different sectors may be related and mass relations of the type |11| 



^baryon ^ ^ ^meson 



(1.11) 



can be obtained. 

The rigorous inequalities (1.9a) and (1.9b) derived via the euclidean path integral formulation amount to the well- 
known fact that the pion is t he ligh test hadron. These inequalities have, however, profound implications for the phase 
structure of QC D. Equation (1.9a) implies no spontaneous breaking of vectorial global (isospin) symmetries in QCD. 
Equation ( 1.9b ), along with the 't Hooft anomaly matching condition, proves that the axial global flavor symmetry 
must be spontaneously broken. 



It is important to note that the basic feature of positivity of the determinant factor (L6) and the functional path 
integral measure d/i(A) is common to all gauge, QCD-like vectorial theories with Dirac fermions. This has far-reaching 
implications for composite models for quarks and leptons. A basic puzzle facing such models is the smallness of the 
masses of the composite quarks and leptons: rrie — rn^ — ~ MeV in comparison with Ap > TeV, the compositeness 
("preonic") scale. A natural mechanism for protecting (almost) massless composite fermions is an unbroken chiral 
symmetry, 't Hooft and others p3| , using the anomaly matching constraint, formulated some necessary conditions 
for such a realization of an underlying global chiral symmetry in the spectrum of the theory. Together with the mass 
ineq ualiti es, th ese c onditions rule out all vectorial composite models for which fermion-boson mass inequalities like 
Eq. (1.9b) [or (1.11)] can be proven. 
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In this work we will mention, at one stage or another, most of the papers written on QCD inequalities, or on the 
related subject of inequalities in potential models. Particular attention will be paid to the seminal works of Weingarten 
1^, Vafa and Witten ]^ - | ^, a nd Witten - all utilizing the e uclide an path integral approach. Weingarten proved the 
inequalities ( 1.9a ) and (|r.9b[) and pointed out the relevance of ( 1.9b ) to spontaneous chiral symmetry breaking (SxSB) 
in QCD. Vafa and Witten directly used the measure positivity to prove tha t par ity and global vectorial symmetries 
like isospin do not break spontaneously 0,^. Finally Witten proved ( 1.9c ) and the interflavor relation ( l.lOj ) 
which holds rigorously for the case of pseudoscalars, with no need for the flavor symmetry assumption. 

To date, QCD inequalities have been mentioned in approximately 600 papers. Most authors were concerned with 
symmetry breaking patterns, the motivation for the Vafa- Witten paper , which is cited most often. Here we equally 
emphasize the other facet of the inequalities, which constitute useful, testable, constraints on observed (and yet to 
be discovered) hadrons. This is why we elaborate on the baryon-meson mass inequalities and related inequalities for 
the exotic sector despite the fact that we have not been able to prove them via the rigorous euclidean path integral 
approach; and on the inequalities between mesons of different flavors, which cannot be justified without the additional 
specific assumption of flavor symmetric mesonic wave function(al)s. 

In general we will follow a didactic rather than a chronological approach. We start in Sec. || by deriving Eq. ( I.IC ) 
using a simple potential model [l4[ and the flavor symmet ry as sumption. This is followed up in Sec. |^ by a potential 
model derivation of the baryon-meson mass inequalities ( 1.11 ) [p 15 1^], and a discussion of baryon-baryon mass 



inequalities in Se c. H. I n Sees. ^ and |6| and App. |^ we show that these relations may be valid far beyond the simple 
potential model |Tl[|l7|] . The elegant analysis of E. Lieb ||l^ of the potential model approach is also reproduced in 
some detail in App.^B[ Section ^ concludes the first part of this review which utilizes the Hamiltonian variational 
approach by verifying that the inequalities indeed hold in the approximately forty cases where they can be tested, 
and by presenting lower bounds to the masses of new, yet to be discovered, mesonic and baryonic states. 

The middle part of the review focuses on the rigorous euclidean path integral approach. We start in Sec. H by 
proving the positivity of the measure and the ensuing inequalities ( |1.6| ), stating that the pseudoscalar pion is the 

lightest meson. We illustrate the power of these inequalities in Sec. || where to|^j > m^^ is used as a shortcut to 
motivate the Vafa- Witten theorem on non-breaking of isospin (which is then presented in some detail). We proceed 
in Sec. |o| with Weingarten's proof of the pion-nucleon mass inequality [Eq. (|l.9b| )] and in Sec. ^ indicate how it 
can be utilized to prove SxSB and discuss its implication for composite models of quarks and leptons. Section |l^ 
presents Witten's proofs of Eq. ( 1.9c ) and of the flavor mass inequality for pseudoscalars. We present also an alternate 
derivation of Eq. (|l.9cD p9| ]. The beautiful Vafa- Witten argument |^ for non-spontaneous breaking of parity in QCD 
is presented in Sec. |13[ Section |l^is concerned with the inequalities in the large Nc limit of QCD [g^. Section |l^ 
is devoted to the glueball sector ]2l|], and Sec. discusses inequalities in the continuum meson-meson sector and for 
exotic states; in particular we also discuss extensions to two-point functions involving local quark combinations of a 
quartic degree ||2^,^. In Sec. w e discuss extensions to finite temperature, finite chemical potential, and external 
electromagnetic fields. In Sec. |18| we discuss the constraints implied by QCD for the chiral Lagrangian approach, 
and also discuss the utilization of QCD inequalities to constrain the QQ potential V{R) in heavy quarkonium states, 
quark mass ratios, and weak matrix elements. Finally Sec. |l^ discusses extensions beyond two-point functions. 

Towards the end of the review we adopt a more heuristic approach, applying QCD inequality-like relations to four 
and five particle states in App. ^ and to electromagnetic effects on scattering lengths in App. We discuss some 
applications of the inequalities in atomic, chemical, and biological contexts in App. |^. Also in Sec. |l^, we use the 
large N (planar) limit to extend the interflavor meson mass inequalities which were ri gor ously proven only for the 
pseudoscalar channel jl^] to other cases. Sections [ij, |l^, |l6|, and also Apps. |^, ^, Q and |^ constitute mostly 
new, unpublished material. Section includes a short summary. We also present two new conjectures concerning 
the possible utilization of the ferromagnetic character of the QCD euclidean Lagrangian, and a possible monotonic 
behavior of mass ratios with the number of quark flavors Nf , as well as inequalities for quantities other than hadronic 
masses. 



2. DERIVATION OF FLAVOR MASS RELATIONS IN A SIMPLE POTENTIAL MODEL 

We first discuss the inequalities in a simple potential "toy model" JT4| , pT|jT^ , |2^ -p6t , which contains some features 
of the full-fledged QCD problem. Specifically the interactions ~ represented here by the potentials ~ are flavor 
independent, and all flavor dependence is manifested only via the masses in the kinetic term. 

Let us consider a two-body system described by the Hamiltonian 

Hi2^Ti+T2 + Vi2. (2.1) 
For a nonrelativistic Schrodinger equation, the kinetic terms are 
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Ti^^. T2^^ (2.2) 
2mi 2m2 

with mi 2 the masses of particles 1 and 2. We assume that the potential V depends only on the relative coordinate 
r = Ti — r 2 (translational invariance), and we also take V — V{\r\) — V{r) only to ensure rotational invariance. 
We can separate the motion of the center of mass: 

V> = e^^-^V^i2(r) 

^ mifi + m2r2 

JX = 

mi + m2 

r = r 1 - r 2 , 

and write _ffi2 as 

with P = pi + P2,P = Pi — P2,M = mi + m2, and u = . For the subsequent discussion we will specialize to 

P = 0, i.e. to the center of mass system (CMS). 

We will be interested in the bound states of hi2 satisfying 

^12V'12 = ei2V'l2 

hi2^^ + n2ir) = ^ + ^ + Vi2(r) + + v{r) , (2.4) 

2/i 2mi 2m2 2mi 2m2 

with V'12 = ilJi2{r) a normalized (J d'^r'0i2(r) — 1) state. (We assume that such bound states exist.) 

We can next consider two additional systems with identical potentials Vii(r) — V22(?') — V{r), but made of two 
particles with the same mass mi (or m2) 

\2mi 2mi y 

,2' 1 1 



To mimic the QCD problem, with quark-antiquark bound states, we still take the particles as non-identical so as to 
avoid issues of statistics. 

Let £12^ , £11^ , £22'' be the ground state energies for the three Hamiltonians 

h.M?ir)^^^,^4?ir)- (2.6) 

We wish to derive the relation 

e^2'>l{e^^'+eg')- (2-7) 



Equations ( p.4[ ) and (2.5) imply the operator identity 

hi2^^{hn + h22) . (2.8) 

Let us take the diagonal matrix element of both sides of this equation with ■012'' i the ground-state wave function of 
hi2- We then have 

(<2^|/^i2|<2^) = ^ \ ((<2Vn|<2^) + {i^f^\h22Wf^)) - (2.9) 

By the variational principle |2^ each of the expectation values on the right-hand side exceeds eii"* (or £22'') respectively, 
which are minima of {ip\hii\^) and {ip\h22\'4>) with ip — tp'f^ and if) — ip^2 ■ Thus Eq. (2.7) is obtained. 
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The previous discussion has been carried out in the CMS frame, Ptotai = 0. In this frame, upon adding the rest 
etc., the inequahty (2^) translates into an inequality for the total masses of the bound states: 



masses toi 



m tnf(°) 
- 1712 to 7 



(0) ^ W (0) 

"^12 ^ 9 r'li 



,(0) 



(2.10) 



Also in this frame total angular momenta are the spins of the composites. 

Let us now make the following observations which will also be very useful for the discussion of Sec. ||: 

1. The addition of the two Hamiltonians hn and ft.22 and their comparison with hi2 may appear to present some 
formal difhculties. In nonrelativistic physics we have a super-selection rule for different masses. Hence strictly 
speaking /in, /i22, and /112 operate in different Hilbert spaces consisting of states with the (1,1), (2,2), and (1,2) pairs 
of particles respectively. However, the crucial observation is that these are identical spaces. All the Hamiltonians 
hii, /i22j ^12 (or any other hij) which are written in terms of f andp can, therefore, be made to operate on the common 
generic Hilber t sp ace of wave functions "0(0 ■ The mass dependence is explicit in hij. It happens to be additive and 
therefore Eq. ( p.8|) is a meaningful and useful operator identity. 

2. The above generic space can be block diagonalized by using the symmetries of the Hamiltonians hij. Thus for the 
case of central potentials, we consider separately subspaces with given total angular momentum I and solve 



Ihj tPij 



13 



(2.11) 



The projection operator Pi used in Pihij{Pi)^ = hfj can be expressed in terms of r and ponly, and is independent of 
the masses. Thus Pi simultaneously projects into the / subspace every term in Eq. ( 



'-12 



'22 



(2.12) 



Using the same arguments as above, we can deduce that 



'(0) > 

tin ^ 



2 V'" 



(0) , m 

~r ^22 



(2.13) 



holds for the ground states in each /-wave separately. 

3. To the extent that heavy quarkonium systems can be treated via a nonrelativistic Schrodinger equation ap- 
proximation, with negligible spin effects, the above derivation would indeed directly suggest mass relations such 



(0) ^ 
as m T > n 

CO — 2 



,(0) 



4. T he inequalities (2/7) and ( 2.13 ) are obtained by taking just one particular matrix element of the operator relation 



(2.8). Evidently the latter contains far more information. In particular one may wonder if the inequalities hold for 



excited states in each given I channel {i.e. for radial excitations). 



The simple generalization, e.g. cj^2 

',(1) ,/,(!) o.^ ,/,(!) 



> 



1 fed) 

2 y^ii 



for the first radially excited state is, in general, incorrect. 



The point is that ^^"^ ^'12'' should minimize the expectation values (-01^^11 IV'): ('01-^221^) and (?A|i?i2|0) 

respectively, subject to different constraints: 0^^'' should be orthogonal to 0j'J'' etc. and the three ground state wave 
functions are in general different. However, such constraints can be avoided if we consider the two-dimensional space 
V2 spanned by the first and second excited states. The variational principal (for hij) tells us that (?/;*^°) |/iy j?/;'^'')) -|- 

('0'-^'' \hij IV'*'^^) = ('^ij' is the minimal value that tr^^/iy can achieve when we consider all possible two-dimensional 

.(1) 



subspaces e|°^ 
states: .g)+e«- 



An-l) 



min^jtr/iy . More generally we have for the sum of the ground state and the first n — 1 excited 
min„^tr/iy , where the trace is now to be minimized over all n-dimensional subspaces 
Considering then the operator relation /112 — \ (/in -|- /122) and taking the trace of both sides in the space v^'^ 
which minimizes tr^^/112, we conclude that 



,(0) 

^12 



.(1) 



1 r 



> - 

- 2 



(4? + 



.(1) 



.(1) 



(2.14) 



with the summation extending up to any one of the radially excited states. 
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5. It is amusing to see how the inequaHties work in famiUar cases. For the Coulomb potential the ground state bindings 
are 



Obviously, B^f < \ [b, 



(0) 



(0) 
(0) 



1 2 

1 /TO, 



1 / TOiTOj 



a 



2 V TOi + TO^ 



(2.15) 



so that the masses of the corresponding states satisfy the desired inequality 



(0) , d(0) ^ 1 / (0) , (0)^ 

my = TO, + TOj - B^j. > - \ m\^ + to}/ 



(2.16) 



Next, consider the harmonic oscillator potential, 5^'^?'^, as a prototype of a confining potential. The ground state 
energies, measured with respect to the minimum of the potential well, are 



,(0) ^ ^^(0) /JL = ^ k( — + — 

2 2 \/ u„ 2 V I TO. TO, 



Wc therefore have 



(^ji )^ + (4j -'^ which implies by simple algebra that 



u - 2 



2 V " y 



(2.17) 



and upon adding the rest masses: 



(0) ^ 
< ^ 2 



,(0) 



(2.18) 



We note in passing that the energies for these two systeins happen to have quantization rules of the form e^") ~ 

and (2n + 1) for the two 



f{n)e^'^\ where /(n) is independent of fi, the reduced mass.|^ Specifically, /(n) 

cases. Hence all the above inequalities happen to hold separately for each of the excited states e-/', e-"'', and e 



An) 



n ■ 



6. While the above discussion was in the framework of a nonrelativistic Schrodinger equation, the result (2.17) holds 

->2 [ — \ 

in a far more general context. The particular form of the kinetic energy Ti = did not play any role in deriving the 

operator relation ( ^.8| ). Hence Ti could have an arbitrary p dependence. In particular, we can take Ti 
the expression appropriate for relativistic motion. 

The key ingredient was that the potentials V12, Vn, and V22 arc all the same, i.e. that we have "flavor independent 
interactions" . However, besides the requirement of translational (and preferably also rotational) invariance, the 
"potential" is not restricted. Thus, V need not depend on r alone, but could have arbitrary dependence on r and p. 



Vij so that hij 



l{hu + hjj) 



The only aspect we have to preserve is that hij — Ti + Tj 

Nonrelativistic quark models have potentials with some explicit flavor (quark mass) dependence, e.g. in the "hy- 
perfine interaction" p8[: 



HF 



(A, - A,) 



(2.19) 



In Sec. I we show that in the full-fledged theory we still have an operator relation 



Hi 



Hki - H-ii + iJ, 



k] 



(2.20) 



'^This feature stems from the fact that these potentials have a power law behavior, r", in which case e^"-* = f{n) is clearly 
manifest in the semiclassical WKB limit. 
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with Hij being the QCD Hamiltonian restricted to a particular flavor sector. With i ~ k,j = I the last relation 
reduces to 

Hij + Hji — Hit + Hjj , (2-21) 



which is very reminiscent of Eq. (2^) and, to the extent that the wave functions i/'ij are symmetric in flavor, again 
leads to the same conclusion. 

One may also avoid the mass dependent si ■ S2 ot s ■ L interactions by considering the "centers of mass" of the 
various multiplets, e.g. (3mp + TOot}/4, representing the mass values prior to hyperfine splittings and applying the 
inequalities to these combinations llSlJlq ]. 

We note that the hyperfine interaction is strongly attractive for the spin singlet p seud oscalar meson case. Since 
furthermore m^^ + to^^ > 2(mim2)~^, this extra binding enhances the inequality ( ^.7[ ) derived for the spin inde- 
pendent part of the interaction. Indeed, the only case for which relation ( ^.71) was proved by utilizing the rigorous 
euclidean correlation fimction approach, is that of the pseudoscalar mesons ||10||. 



3. BARYON-MESON MASS INEQUALITIES IN THE GLUON EXCHANGE MODEL 



The simple "quark counting rules" relating meson and baryon total cross-sections 1 29 30|] and masses [ pl[ were early 
indications for the relevance of the quark model js^jQ. If the mass of a ground state hadron is just the sum of the 
masses of its quark constituents, then 

2-s=-iI^+-ff+™2^^ (3-1) 

with TOj-^], the mass of the ground state baryon consisting of quarks qi,qj,qk; and m^^'' the mass of the lowest lying 



qiQj meson. The following discussion motivates a related QCD inequality |ll| , |l5| , |l8y34 

2r4l>-^+rr^+-^- (3-2) 

The particular variant > was derived by Weingarten Q from the correlation function inequalities, and will 
be discussed in Sec. 

If qq (qq) interactions are generated via one gluon exchange, then the color structure is V12 oc Ai • A2, with A a 
vector consisting of the iV^ — 1 Gell-Mann matrices of the fundamental SU(N) representation. In a meson the quarks 
qiq2 couple to a singlet and, using (• • •) to indicate expectation values, 

0= ((A1+A2)') meson 2A2 + 2(Ai-A2) meson 5 

(3.3) 

so that 

(Ai • A2)m = -A^ . (3.4) 

Here A^ = J2n=o^ ('^n)^> is a fixed diagonal N x N matrix proportional to the unit matrix. The SU(N) baryon is 
constructed as the color singlet, completely antisymmetric combination of N quarks 

ea,...a„g'^^- "" . (3.5) 

Thus 

0=/(^AJ \ =iVA2+^(A,A,)B=7VA2+iV(iV-l)(Ai-A2)B, (3.6) 



where we used the fact that all the N{N — 1) expectation values (A,; • \j)B are equal. Thus using Eq. (3^) we find 
that 

This implies that the strength of the one gluon exchange interaction in the meson qiqj is, for QCD (i.e. = 3), 
precisely twice the corresponding gluon exchange interaction for a qtqj pair in a baryon. 
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Since (Ai • A2) b = ^{^1 ■ ^2}m is just an overall relation of the color factors, we have the corresponding ratio of the 
pairwise interactions in the meson and in the baryon: 

%(---)m = 2%(---)b, (3.8) 

as long as the qiqj are in the same angular momentum, spin, radial state, etc. as the qiqj in the baryon (or in the 
appropriate mixture of states). 

For N = 3 and a general two-body interaction, we write the Hamiltonian describing the baryon as: 

Hb = H^jk ^T,+ T, +Tk + + Vjl + Vi , (3.9) 

with Ti the kinetic, single-particle operators. 

The qiqjqk baryonic system can be partitioned into (2 + 1) subsystems in three different ways: 

{{ij), k), {i, (jk)), {j, (ki)). Let us consider the three mesonic qq systems corresponding to these partitionings: 
Mij, Mji^, Mkj, composed of {qiqj), {qjqk), (lk<ii)- The Hamiltonians describing these mesons are 

Hk, = Tk+Ti + V^ . (3.10) 

Prom V^'^ = etc. it is easy to verify the key relation 

2i?ij/s = + Hji. + ilfcj (3.11) 

between the Hamiltonian describing the baryon Bijk and the three mesonic Hamiltonians corresponding to its diquark 
subsystems. Let us next take the expectation value of the last operator relation in the normalized ground state i^fj^ 
of the baryon at rest [to simplify the notation we suppress Lorentz {J^'^) quantum numbers]. The left hand side 
yields 

H4fk\Hijk\4jl}=^rnfk, (3-12) 
with m|°]. the mass of the ground state baryon. The right hand side is 

{>pSl\Hij\^Sb + + {4jl\Hk,\4jl} ■ (3-13) 

The matrix elements of the two-body operators in the three-body wave function, e.g. {^iljl\Hij\ipflll) . are evaluated 
by considering Hij as a three-body operator which is just the identity for quark k. Thus we view 'tjjijk, for fixed 
coordinates of the quark qk, as a two-body (qiqj) wave function V'^^^, and compute (-(/'Ij'' l^/'^'). Finally, this is 
integrated over all values of the coordinates of qk. 

The key observation is that "0(0^ the two-body wave function prescribed by V'^jli is, in general, different from V-'(^j): 
the ground state wave function of the meson qiqj. By the variational principle the latter wave function minimizes the 
expectation value of Hij, the mesonic Hamiltonian: 

min^^^_(Viij|J?ij|4-) = i^i? iHijli'l?) = , (3.14) 

where 'm\j\ the energy of the state at rest, is the mass of the ground state qiqj meson. Hence we have 

(:4>l^\HiSl?}>r4? (3.15a) 

and likewise 

{^^\H,,\i>^)>rr^ (3.15b) 

(^E^|iffcdV!E^)>4?- (3-15C) 
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Since each of the se inequahties persists after the normahzed integration over the coordinates of the third quark [e.g. 
Qk for Eq. ( 3.15a )], we have, 



{^l°i\H.j\i^'ub>m^, etc. 



(3.16) 



Equating ( 3.12 ) and ( 3.13 ), the matrix elements of the left hand side and right hand side of the original operator 
relation, we arrive at the desired inequality (^^) 0,^. The following remarks elaborate on these inequalities and 
the conditions for their applicability. 

1. We have not displayed the quantum numbers of the baryonic tpijk or the mesonic ifjij states. As indicated by 
the above construction of "trial" mesonic functions from the baryon wave function, tpij must be in the state (or 

in general, in a mixture of states) prescribed by the original baryonic wave function t/;^^^ . 

Consider first the A"*""*" baryon. In the approximation where L 7^ components in the ground state wave function 
are ignored, the state with Jz{A'^^) = 3/2 consists of three up quarks with parallel spins: u "f it f u t, with the arrow 



indicating spin direction. Each of the will be, in this case, u "f it |, in a spin triplet state p or lu. The inequality 
reads 



TOA++ > -^mp (or m^) 



(ijk = M 1 u I u "f) 



(3.17) 



Likewise, 



(ijk = s t s t s t) 



(3.18a) 



2771=0 > 2mK* + rn^ {ijk = s t s T T) 



(3.18b) 



{ijk = M t u t s t). 



(3.18c) 



The situation is different for the J ~ 1/2 {i.e. S* = 1/2 in this L = approximation) baryons such as the nucleon (see 
App. ^). In this case the diquark systems are, on average, with equal probability in the = and = 1/2 states, 
as can be readily verified from ((si + S2 + ss)^) = (s^) = 3/4. [Strictly speaking, the uu diquark is pure triplet and 
the ud diquarks are in a triplet (singlet) state with probabilities 1/4 (3/4), respectively.] 



To see the effect of having a mixed rather than a pure trial mesonic state, we revert back to Eq. (B.15a). Instead 



of a single matrix element of Hij in a specific V'(ij) state, we have, in general, a weighted sum of matrix elements 
corresponding to the mixture of two-body states generated from the three-body baryonic wave function: 



(3.19) 



Here n labels the different {J^) states and are the normalized weights (^„ c^ — 1). 

The variational argument can be applied separately for each of the (trial) ^[n) states and to their matrix elements 

{'4^{n)\Hij\^(n)) to obtain 



i]{n) 



(3.20) 



For the specific case of the nucleon the weights of the qq singlet (triplet) configuration are c^^^^^^.^ = c^^^^-j^-^^g = 1/2- 
The corresponding qq states are the p{lo) and the tt, and the inequality ( 3.20 ) (plus similar ones for the matrix elements 
of Hjj. and Hki) yields: 



Stoat > -(777^ 



■777p) 



(3.21) 



Likewise by considering the A hyperon we find 

2777A > i [777^ 



3777 



K 



(3.22) 
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2. A simplified version of the inequalities applied just to the spin-averaged ("center of mass'') multiplets was suggested 
by Richard Note that unlike for the case of the meson-meson mass relation in Sec. |2[ the baryon-meson mass 
inequalities do apply even when we have mass dependent qiqj (or qiqj) interactions such as the hyperfine and/or 
spin orbit interactions generated by one gluon exchange. The point is that the same flavor (mass) combinations 
(jj), (jk), (ki) appear in both the diquark subsystems and in the mesons Mij,Mj^, and Mki- Hence unlike the case 
of the meson-meson relation there is no motivation to consider only the above simplified version. 

3. The above discussion utilized the nonrelativistic quark model with a one gluon exchange potential. We show next 
(and in Sec. ||) that the inequalities are valid in a vastly larger domain [0. 

First we note that any vertex or propagator insertions into the one gluon exchange diagram leave the Ai • A2 color 
structure intact (see Fig. |^). Such insertions generate a running coupling constant a(q^) or a(r). Even if the non- 
perturbative effect of the propagator insertions can generate a confining potential V ~ ar, this potential still 

has the Ai • A2 color structure, and the derivation of the inequalities still applies. 

4. The most general two-body interaction [due to any number of gluon and qq exchanges; see Fig. || in the qq or qq 
system, has, from the t-channel point of view, a3(8)3 = 8®l color structure. The octet part corresponds to the 
Ai • A2 structure; the singlet to 1-1. Thus we need not assume that the important interactions are due to one gluon 
exchange. Rather, we have to assume that the octet, Ai • A2 part of the full two-body interaction is dominant. We 
note that in the large Nc limit we expect the Ai • A2 part to dominate over the 1 • 1 part. 

5. The trilinear gluon couplings are a potential obstacle to having the operator relation ( ^.11 ) in a general setting. 



These couplings could yield genuinely nonseparable interactions in the nucleon of the type shown in Fig. ||. However 
it turns out that the trigluon diagram vanishes in the baryon state. The relevant color factor is 

frst \aa \bb' \cc 
J ^abcf-a'b'd , 

with rst = 1, . . . , 8; ahc {a'b' c') = 1, . . . , 3 the adjoint and fundamental color indices, respectively. This color factor 
vanishes since an exchange r— !-s,a— >6, a'— s-fe' changes the sign of the summand. This cancellation is not modified 
by any further dressing of this diagram with more gluon exchanges. 

6. The above discussion did not require any specific form of the one-body kinetic terms T^. We could have the 



nonrelativistic Ti relativistic Ti = \/pf + rnf form, or, by considering the operator relation (|3.1l|) as a 

matrix in spinor space, also a Dirac Ti + rrii form. 

7. Just as for the meson-meson relation, we could use the variational principle in the space of the lowest n states to 



obtain a relation between masses of radially excited baryons and mesons analagous to Eq. (2.14) 



A diquark in a baryon cannot annihilate into gluons. This is not the case, however, for the uu or dd in a meson. 



These annihilations are avoided by choosing 1 = 1 [i.e. p rather than lo in Eq. (3.17)] combinations. 

9. A simple intuitive feeling for the deviation from equality expected in the baryon-meson relation has been offered 
by Cohen and Lipkin |^ (in a paper which predated the QCD inequalities). The point is that the diquark systems 
in the baryon are not at rest but recoil against the third quark with a typical momentum p of a few hundred MeV/c. 
Thus the trial meson wave functions correspond to a meson moving with p 7^ 0. Consequently 

2'mijk > Eij + Ejf. + Eki 

~ +Pk + \/^% +Pt + \J^M +Pj > + ^jk + rriki . 

A more comprehensive investigation of the pattern of deviations from equality as a function of the quark masses, 
and in particular for logarithmic interquark potentials was recently performed by Imbo ||3^ . 

10. The baryon-meson relation can be easily extended via Eq. ( |3.7D to any number of colors N: 

mf > - 



4. BARYON-BARYON MASS INEQUALITIES 

The meson-meson mass inequalities above (see Sec. ^ follow in potential models only if the various two-body 
quark-antiquark potentials are independent of the quark masses. Baryon-baryon mass inequalities motivated by 
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similar convexity arguments were suggested |17 originally to also hold when all two-body quark-quark interactions 
are flavor-independent. It turns out however, that inequalities of the form 



) + (m, M, M) < 2E'-"^ (m, m, M) 



(4.1) 



do not hold in general |p^ , p5| , ^ . As will be clearly indicated in Sec. the key requirement for proving either meson- 
meson or baryon-baryon mass inequalities is the flavor symmetry of the ground state wave function. In the two-body, 
pure potential case, the relevant, relative coordinate wave function '0(*'^(r),r = ri — r2 is guaranteed to have the 
T^2 "flavor" symmetry. This however is not the case for three-body potential systems 



H = Ti{pi) + T2(p2) + 73(^*3) -I- Vi2(ri - f2) -I- V23(r2 - fs) + Vsiif^ - n) 



(4.2) 



Indeed even for flavor-independent potentials V12 — V23 = V31, and simple nonrelativistic kinetic terms Ti = pf /2mi, 
the ground state wave function il>''^\f 1^x2^^^) is not (flavor) symmetric under interchange of 1 ^ 2, etc. As various 
counter examples show |l^,^ (particularly the simplest one due to Lieb which we present in App. ^), for 
certain, rather "extreme" potentials Vij{r) = V{r), the kinematic assymetry due to the different quark masses 
(mi = 7713 = M, 7712 = "i) is strongly enhanced s o th at Eq. (4.1) fails. This notwithstanding, the elegant work of 
Lieb [|l8| has shown that the conjectured equation (4T) does hold for a wide class of two-body potentials V(fi, fj) for 
which the operator 



Lp{.x,y) = e 



(4.3) 



is positive semidefinite. The latter is equivalent, for V = V{x — y) and Lp = Lp{x — y), to a positive semidefinite 
Fourier transform of Lp. A sufflcient condition for this is that V{ri, fj) — V{fi — fj) — V{r) satisfies 



V'{r) > (monotonically increasing) 



(4.4a) 



V"{r) < (convex) 



(4.4b) 



and 



V"'{r) > 0. 



(4.4c) 



Intere sting ly for the case of heavy quarks (where potentials can be derived via the Wilson loop construction), 
Eqs. (4.4a) and (4.4b) can indeed be prov en, as we will show in Sec. |l^. Indeed all potentials used in quark model 
phenomenology to date satisfy all of Eqs. (4^) and hence (4.3). 



We present in App. fq Lieb's proof of Eq. (4.1) for positive semidefinite exp[—pVij{x,y)]. We do this in some 



detail, even elaborating somewhat beyond the original concise paper, since the baryon-baryon inequalities are indeed 
born out by data. Also, Lieb's line of argument invoking the full three-particle Green's function serves as a "bridge" 
between the Hamiltonian, variational, largely potential model motivated, first part of our review; and the more formal 
Lagrangian correlator inequalities proved via the path integral representation which is the approach of the second 
part. 



5. RELATING MASSES IN DIFFERENT FLAVOR SECTORS 



In the following we investigate the interfiavor mass relations ( 2.1^ ) in a lattice Hamiltonian formulation of QCD 



The conservation of quark fiavors allows us to break the QCD Hilbert space into flavor sectors. Each flavor sector 
([/, D, S,C, . . .) consists of a net number U of up (u) quarks (U is negative for u excess), D of down (d) quarks etc. 
with U + D + . . . = (mod 3). We will be interested in the following "low" sectors: 

1. M(o): The flavor vacumir U = D = ■ ■ ■ = 0. It consists of states with an arbitrary number of gluons and {qiqi) 
pairs. 

2. -M(jj), i j: The meson sector, with an excess of one quark flavor of type i and a different antiquark flavor j 
together with any number of gluons and (qiqi) pairs. 

3. B(^ijk): The baryon sector with a net excess of three quarks, g^, qj, and g^. 

4. M f^^rjj^^{i ^ i,i ^ l,k ^ j^k ^ I): The exotic meson sector with a net excess of two quark flavors and two 
(different) antiquark flavors. 
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We will also discuss in Sec. 16 and App. ^ other sectors such as the pentaquark and hybrid sectors. 

We can now define Hij to be the QCD Hamiltonian restricted to the sector M^j; Hijk the Hamiltonian restricted 
to Bijk etc. The meson spectrum will be given by 

H^^.j) ^m,0^j) , (5.1) 

with 5* a wave functional belonging to the Af^ sector. Similar equations hold in other sectors. 
We show in App. ^ that 

H,, + H,i=Ha + Hkj. (5.2) 
Particle masses are given by the Schrodinger equation, e.g. 

H^j\'^^J) ^ (my + Ao)|^'»j) , (5.3) 

with Aq an additive, common constant representing the vacuum energy, i.e. the lowest energy obtainable for function- 
al in M(°l To ensure that Aq is finite we restrict ourselves to finite lattices so that Aq = T^eiac, with V the volume 
and eiac the vacuum energy density. In writing Eq. (5.3) we assumed that the systems have no net translational mo- 



tion.0 Additional symmetries (rotations, parity, and for certain states, charge conjugations) could be used to project 
any desired state of given J^'-^ quantum numbers. To simplify notation, the J^*^ labels are omitted. The lattice 
formulation reduces the symmetry from the full rotation to the cubic subgroup. However, the operator relations most 
likely also hold in the continuum limit where the full rotation symmetry is regained. 
The l^'ij) are wave functionals: 



E 4onf)l(^°nf))- (5-4) 

[ifif^uratioiis 



Even for a finite lattice the (discrete) summation includes infinitely many configurations characterized by the locations 
and spinors of all quarks and/or antiquarks and by E'i at all links where the latter are constrained by Gauss' 



law [Eq. (C4)] (recall that arbitrarily high SU(3)c representations are a priori allowed). The specific flavor (ij) 
dependence enters only into the probability amplitudes ^('Joj^f^ for finding a given configuration in The generic 

configurations used in App. ^ have been defined in such a way so as to make operator relations such as Eq. ( |5.2| ) 
manifestly true. is normalized: 

l4onf)l' = l- (5-5) 

configurations 

A key observation is that the variational principle applies to wave functional solutions of the Schrodinger equation 



(5.3) much in the same way as it does to wave functions. In particular, the ground state l^*!^') in any specific channel 



(say a mesonic ij channel with given J ) is given by the requirement that: 



Ao + m,, = {^f\H,, = min(vl'|i/,, |v|/) , (5.6) 
with the minimum sought in the space of all (normalized) functionals with the given quantum numbers: 

= ^Aconf|(conf (mesonic))) . (5.7) 

conf 

We would like next to argue that we can set i — l^j — k in Eq. ( |5.2| ) and still obtain a meaningful operator relation 

H,j + Hfl - Ha + . (5.8) 



■^To ensure P — i.e. translational invariance, we need to sum over all locations of the centroid of the wave functional. For 
finite lattices of size L, only P < 1/L seems achievable, yet for periodic boundary conditions a discrete version of translational 
invariance persists. 
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In the light (ud) quark sector |m,i — m^l ^ Aqcd, and since oem ^ ctQCD, isospin is a good flavor symmetry. The 
1=1 \uu — dd + gluons + pairs) states do not then mix with the 1 = 0, M^"^ sector, so these states should then be 
used in_deriving the inequalities. For heavier states ss, cc, bb there is a strong "Zweig rule" p2| , ^j4l| suppression of 
ss, cc, b b an nihilation. He nce mcsonic sectors Mi^ distinct from M^^'^ can be defined for which the i = l,j = k version 
of Eq. (U), namely Eq. (U), holds. 

Unfortunately, we cannot use this relation to obtain mass inequalities without further assumptions. In general, 
^E'l^'' , the ground state wave functional for (for example) a heavy quark qi and a light antiquark qj , may be different 

Indeed one of the inequalities, namely 

(0) 



(0) 



from 



even though ^'y and ^ ji are related by charge conjugation.] 



2m K' < + TOp, is not manifest. Symmetry is automatically guaranteed in the simple potential model where ^' 
depends only on r = fi —fj, the relative coordinate of qiqj. It is also plausible in a large Nc limit where the important 
degrees of freedom are the gluonic ones. In the following we will assume ? <-> j symmetry. 



By taking the expectation value of the left hand side of Eq. (5^) in ^Ij'', we have (vf^^ 
and the variational argument for deriving the inequalities proceeds exactly as in Sec 
2Ao term representing vacuum energy we obtain 



(0) 



H. 



2mg^+2Ao, 



After cancelling a common 



o (0) ^ (0) 



(5.9) 



We also have the relations for the sum of the first n excited states in any M'HjP'^) channel [Eq. ( pls] )] 
For the baryonic sector we readily find a relation 



Hklr — H. 



ilr 



H 



jkr 



(5.10) 



Eq. ( 5.10 ) is analagous to relation ( p.2| ) and is obtained by adding a common "spectat or" q uark r ^ i,j,k,l. 

Assuming a flavor symmetric baryon ground state, we obtain from variants of Eq. ( [3.1C ) with some of the flavor 
indices set equal to each other, analog convexity relations for baryonic states 



.(0) 



,(0) 



(5.11a) 



(0) ^ 1 



,(0) 



(0) 



'-fefei 



'■kkj 



(5.11b) 



It should be emphasized ^] that for the baryonic case we have two relative CMS coordinates, and flavor symmetry 
of the wave function is not guaranteed even in the framework of a potential model. Indeed, as di scus sed in Sec. ^ 
and App. |^, for certain extreme type of potentials one can show that the inequalities such as Eq. (5^) are violated 
- though as elaborated in App. ^ they do hold [|l^ for the class of potentials which are of interest in QCD. 



6. BARYON-MESON INEQUALITIES IN A NON-PERTURBATIVE APPROACH 



Inequalities relating masses of baryons and mesons were motivated in Sec. |^ via a potential model. It was suggested 
there that such inequalities may persist beyond the (dressed) one gluon exchange approximation. In this section we 
will employ a nonperturbative framework which still allows the application of variational arguments and the derivation 
of the baryon- meson mass inequalities [po| . 

It has been argued that confinement in QCD is the electric analog of the Meissner effect in a superconductor. The 
nonperturbative QCD vacuum develops a condensate of color monopole pairs and/or of large loops of magnetic flux 
p3|-^5[. More recently, lattice and other approaches have made this much more concrete, particularly in the context 
of "center vortices" p^ , ^ . In this vacuum, the chromoelectric flux emitted from a quark and ending on an antiquark 
is localized along a thin tube |^] - the analog of the magnetic vortex in the ordinary superconductor. In the — > oo 
limit the chromoelectric flux tube may become infinitely thin and the original dual string model |]50|-p^ for 
hadrons could emerge. Also in the strong coupling limit of Hamiltonian lattice QCD a single set of minimally excited 
links of total minimum length connects the quarks and antiquarks in a single hadron. 



^This important point (missed in Ref. ||17|) will be further elaborated in the Summary section. 
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A more general approximation for the meson wave functional is a sum over configurations of strings [or chains of 
lattice hnks, see Fig. |j(a)] connecting q and q (indicated by the symbol ^ below), with a probability amplitude for 
each configuration: 

|^M12) = ^A(^)I^) . (6.1) 

The end points 1,2 can also vary. The kinetic {B^ and 'tpDip) parts of the QCD Hamiltonian move the string and 
the fermions at the end points, respectively; while the "potential" term gives a diagonal contribution proportional 
to the length (number of links) of the string multiplied by the string tension (weighted by the values of the second 
Casimir operator for the representation residing on each link). 

In the same approximation, the baryon's wave functional is a sum over configurations where the three E vortex 
lines (symbolized for convenience by Y in the following equation, though the vortices need not be straight) emanate 
from quarks 1, 2, and 3, and join together at a common "junction point" x [see Fig. ^(b)]: 

|^B123>=E^(^)I^)- (6-2) 
Y 

Just as in the case of the nonrelativistic potential model, we would like to extract trial wave functionals for the 
ground states of the mesons qiq2, 9293, and q^qi from the ground state baryon wave functional. This is achieved in 
the following way (see Fig. ||). We consider quarks 1 and 2, together with the string connecting them in the baryon as 
a possible configuration in the wave functional of the meson M12. To this end we need to color conjugate the quark 52 
and reverse the chromoelectric flux in the section a; — 2. A detailed analysis using a lattice formulation indicates this 
can be done for the case of minimal flux strings, without effecting the matrix elements of the Hamiltonian. Likewise 
quarks 2 and 3 with the string 2 — x — 3 are viewed as a possible configuration for the meson A/23 and similarly for 
M31. 

This suggests the operator relation 

2Hijk = Hij + Hjj. + Hki , (6.3) 

with ijk the flavors of 91 (7293- Indeed by applying iJy, H^^. and H^i to each of the above mesonic configurations we 
see that each part of the baryonic Hamiltonian Hij^ is encountered twice. Thus the kinetic term and mass term 
{iPi + mi) associated with the motion of qi{— qi) occurs in both Hij and i/fej, and so does the kinetic and potential 
energy associated with the motion, and total length, of the string bit ^ connecting quark qi with the junction point 
X. A similar argument can be applied to the 2 — x an d 3 — x string bits. 



Let us next take the expectation value of Eq. (6.3) in the ground state wave functional of the baryon at rest. As 



in Sec. ^, we obtain one one hand simply 2to-^]. and on the other hand the sum of the expectation values of Hij, Hj]^ 
and Hki in the mesonic wave functionals extracted, in the manner described above, from the baryon's ground state 
wave functional. In analogy to fixing the coordinates of the "third quark" (say qk) in the potential model discussion 
of Sec. I we fix in each case the string bit {x — qk for Hij) which should be ignored in order to form a mesonic (Mij) 
trial wave functional. The expectation values of Hij in these trial wave functionals 

Y 

are then integrated with the weight implied by the original baryonic amplitude A{Y) over the "ignored" section as 
well. 

N ote that because of the special character of the baryonic state, the configurations in the trial wave functional 
( |6.4[ ), used here for the three mesonic qiqj, qjqk, IkQi ground states, are in fact correlated. Specifically, the strings (or 
fiux lines) associated with all three of these meson states are forced to have a common junction point x (which again 
is integrated over at the end). 

This extra constraint only reinforces the general result of the variational principle, namely that 

(^Slm^^l}>m^ (6.5a) 
{4-l\Hfk\4jl)>m^jS (6.5b) 
(V'I-l|i?fe.|V'I-l)>mi°). (6.5c) 
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Specifically these inequalities state that the mesonic trial wave functionals extracted from the subsystems of the 
baryon's ground state wave functional are not optimized so as to minimize the energy of the mesonic subsystems. 
Rather ipl^jl is constructed so as to minimize the expectation value of Hijk and 

(^Sl\Hvk\^Sl^=rr^- (6.6) 
Thus we have from Eq. (|6.3| ) the required inequality 

2m^>m^+n^+n^. (6.7) 

It is amusing to see how this inequality is realized for the strong coupling limit. The flux then proceeds in straight 
lines in both the mesons and the baryons so as to minimize the total string length. This yields the potentials 
Vi2 — cr|ri — ^^21, V23 = (j\r2 ^ ?^3|, V31 ~ crjra — fi\ for the meson and the genuine three-body interaction ||53|| 

V123 = minj;cr(|fi - x\ + \r2 - x\ + jfs - x\) 

(with the total length of the Y-shaped string configuration minimized over the choice of function point x) for the 
baryon. The inequality reduces in this case simply to 2V123 ^ ^12 + V23 -I- V31 which is just a sum of three triangular 
inequalities: |ri — x\ + |r2 — x\ > jri — r2| etc. (see Fig. ^). Since the potentials are linear for both systems, the virial 
theorem implies that the total energy is given by 2(y), with {V) the expectation value of the potential energy, and 

2-El23 > ^'12 + E23 + £^31-3 

We have so far considered the approximation in which the quarks are connected by a single set of minimally excited 
links on the lattice or by a single Y-type configuration for the baryons. We still have many paths of different lengths. 
All of these configurations can be generated by repeated application of the kinetic {B^) term in the Hamiltonian which 
generates a closed flux line around a plaquette and shifts the initial flux line as indicated in Fig. ^ The problem of 
finding the mesonic or baryonic ground state wave functionals, even in this approximation, is intractable, and at best 
we could hope for some numerical results. We would like to point out, however, that the baryon-meson inequalities 
can be proven in an even broader context when qq pair creation, bifurcation of flux lines, and excitation of links to 
higher SU(3)c representations are all allowed. 

If we consider the procedure for extracting mesonic trial wave functionals from the baryonic wave functional, we 



can pinpoint the crucial ingredient for deriving Eq. (6.3). It is that the network of links in any configuration in the 
baryonic wave functional includes only one junction point x, where we can separate the network into three "patches" 
Pi,P2, and P3 connected to the external quarks qi,qj, and qk respectively (see Fig. ||). We consider P1P2, where P2 
is the patch with all flux flows reversed and qj — > qj, as a configuration in the trial meson functional \ipMij) with an 
amplitude A{Pi, P2, P3), (P3 fixed) inferred from the baryonic wave functional. Likewise P2, P3 and P3, Pi offer trial 
wave functionals for the 9293(9j9fe) and qsqiiqkQi) mesonic systems. We can easily verify that the operator relation 
2^123 = H12 + H23 + H31 still holds for this class of baryonic and corresponding mesonic functionals. The kinetic 
and potential parts of the full QCD Hamiltonian operate on each patch Pi, P2, P3 separately and those contributions 
are counted twice in H12 + -^^23 + -^31 ■ 

The variational argument is applicable in the larger class of baryonic and mesonic trial wave functionals and it 
yields the required baryon-meson mass inequalities 2mijk > rriij + rrij^. + m^j. 



The fact that the same inequality, Eq. (3.2), motivated by the color factor for the perturbative one-gluon exchange, 
can be rederived in a strong coupling, string-like limit (and some generalizations thereof), suggests that these in- 
equalities may indeed be a true consequence of the full-fledged QCD theory. Indeed a particular relation involving 
the nucleon and pion masses has been rigorously proved by Weingarten g|, who utilized the correlation function 
techniques, and is reproduced in the second half of this report.^ 

Finally, we note that a construction of mesonic trial wave functionals from the baryon's ground state functional 
can be carried out for SU(N) with general N as well. The configuration contributing to V'(n...ijv) consists of 
N strings joined at a common junction point. It can be formally separated into N(N-l)/2 mesonic subsystems with 



■*The virial theorem (with massless quarks and linear potentials) was used in [gij to motivate the equipartition of the light cone 
momenta between quarks and gluons, and more recently in connection with the mass dependence of Bose-Einstein correlations 
in multi-particle final state s [pSi . 

^Despite many eff'orts |^,^jpto extend the correlation function technique to more detailed, flavor-dependent baryon-meson 
and baryon-baryon inequalities, no fully convincing results have been obtained. 
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each part of the baryonic Hamiltonian counted (N-1) times (since a segment or patch Pia connecting the quark qi^ 
and the junction x appear in the (N-1) qi^Qi^, a — b systems), and we have: 

(iV-l)i/,,..,„ = ^ H,^,^, 



leading to the baryon- meson inequaUty ( 3.23| ) for general SU(N) 



The case of TVc = 2 is rather special. The lightest diquark baryon is in fact a 0++ mesonic state. The inequality 
rriB > Tt^r is in this case most likely an equality: 

(6.8) 

Indeed the giuon couplings inside the meson and "baryon" are the same here, so that to all orders in perturbation 
theory we expect the S-wave qq (or qq) states to be degenerate. To obey the generalized Pauli principle, the S-wave 
spin and color singlet qq state ought to be a "flavor" antisymmetric ud combination which is to be compared with 
the wy^d pion. The pseudoreality of the SU(2) group implies that in the nonperturbative string or flux tube picture, 
the flux emitted by one quark in the bosonic diquark can readily end on the other quark. There is no junction point 



in this case, the wave functionals and Hamiltonians for the q^^q and 0"*" qq systems are identical, and Eq. (6.S) 
follows 



7. COMPARISON OF THE INEQUALITIES WITH HADRONIC MASSES 

We proceed next to list the various inequalities and compare them with available particle data [^ . 

First consider the meson-meson relations. As emphasized above, these inequalities do rely on an additional as- 
sumption of a flavor-symmetric wave function. In all testable cases (with one possible exception) these inequalities 
are satisfied. The inequalities could then be used as a rule of thumb to restrict the masses of as yet undiscovered new 
particles. 

There are six relevant rriij, ^ flavor combinations: us, uc, sc, sb, ub, and c5.^ Because of the small violation 
of J-spin symmetry (via radiative electromagnetic corrections and the effect of |m„ — md\ — 4-5 MeV) we have not 
separately considered members of an /-spin multiplet (obtained by u ^ d substitutions), but rather averaged over 
these states. Also in order to minimize the effects of the annihilation channels we consistently choose the I = 1 uc? 
rather than the 1 = Q uu states {e.g. p and not w, a2 and not /o etc.). 

The comparison with measured masses is summarized in Table |, which indicates the specific particles and masses 



relevant for each inequality. In some cases all three masses in Eq. (2.18) are fairly well known or reliably estimated 



In other cases lower bounds are predicted for certain M^^^^"^ . All masses are listed in MeV. 

For the pseudoscalars the inequalities have been derived by Witten |l0[ using the euclidean Green's function method 
to be discussed at length in the second part of this review, and indeed hold with a fairly wide margin. The light 
pseudoscalar mesons tt, K, rj can also be viewed, to a good approximation, as pseudo-Goldstone bosons. Current 
algebra methods jso) then lead to rripg.. cx (wg^^ -I- rnf.^). Hence in this approximation 'i.m'^^.. = ni^^.^ + m^^.. and 

2mpsi3 > "^psii + '"T'psjj is guaranteed.!] 

For the 7 = pseudoscalars the strong coupling to the gluonic M° channel (which in particular accounts for the 
massive r]') |^,^ suggests strong mixing between the two 7 = states that are made of light quarks, -^{uu + dd), 

and the ss state. Hence \ss) = a\r]) + I3\ri'). The reasonable expectation [M that the eigenstates 77,77' are the SU(3) 
flavor octet -^(uu + dd — 2ss) and singlet -^{uu + dd + ss) implies a — \/2/3, f3 = \/T/3 and 

mir^ = + PW\) H (a|77) + f3\r^')) = a^rrir^ + (3^m^, . . . 



PC 



^Since the decay width of the top quark exceeds a few GeV, it decays before tu, etc. states can form. 

''The last argument is reminiscent of that made earlier for harmonic potentials. It was noticed that the "dual" Lovelace- 
Shapiro- Veneziano formula for vr — vr scattering has a remarkable tendency to conform to soft pion theorems. Considering 
the harmonic string origin of "dual" amplitudes, these may perhaps be related features. 
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and yields m^j ^ « 706 MeV. This indeed satisfies m^j + m||^ < 2to^°^ ■* [i.e. m^j + m^r < 2TOi<-] with a large 

margin. In general 2m k > rriT^ + rnfg •* holds as long as fP' jo? < 1.5, a constraint satisfied by all mixing schemes 
suggested for the 0~ mesons . 

In the relation 2m^^^ > ttiuu + nif^i we use the mass of the vectorial state, tot, instead of the as yet unknown 
TO^j^. Following Weingarten we show in Sec. ^ that, when mixings with are neglected (which seems justified 
by asymptotic freedom in the case of heavy quarks), the lowest state in any Mij sector is indeed pseudoscalar and 
Tot > fTT-rjh^ so that the original inequality is a fortiori satisfied. 

The vector meson mass inequalities hold with smaller margins than the corresponding inequalities for the pseu- 
doscalars. This reflects the hyperfine splittings (as indicated in the conclusion of Sec. H), which tend to weaken 
(enhance) the inequality for the vector (pseudoscalar) mesons. We find it impressive that despite the very large spin 
splittings, e.g. nip — uItt — 650MeV ~ Stott, the inequality may still marginally hold for the vector mesons. For the 
specific case of K*, p, and (p, the inequality seems to fail, although only within 0(1%) of the widths of the states 
considered. 

The inequalities for the tensor mesons can be tested in the sectors not involving the b quark. In the b quark sectors 
we can use the inequalities to predict lower bounds for the masses of the i?tonsor particles. The -0^2+ is not known for 
certain to be a 2++ state, but this is the quark model prediction. The inequalities for the axial 1++ mesons all hold 
at the level of a few percent. Finally, the inequalities for scalars can be tested (and verified) only for the case of the 
uc and ub flavor combinations.^ 



The deviations from equality 



g^Sm ^ - (1/2)(to,, + TOJJ-) 



TO TO. 



are often small [0(1% — 2%)]. This may be traced to the variational principle. Since {ip\Hii\ip) , {ilj\Hjj\ip) are 
extremized at ip = i'^Ki' = ip^f , respectively, the deviations A,; = {'iPif\Hii\ipif) - (V^i^^l-ffiilV'ii''') ^^d = 



(V'i^l-ffjjIV'i?^) — (V'7-?l-ffjjlV'7-?) are quadratic in the wave function shifts, e.g. 



5m^{5^f^{^t;^-^t:^f. (7.2) 
The 5^1) in turn are expected, on the basis of flrst order perturbative estimates, to be 

(7.3) 

Ato Ato Ato 



(0) (1) 

considered. 5H is solely due to the quark mass differences. Taking i ^ u, j — s we have for the simple model of Sec. f^: 



with Am ~ 1 GeV the typical splitting m\j — mlj ' between the ground and flrst excited states in the specific chann_el 



5H ~ ^m^ + (p2) - + (p2) , (7.4) 

with (p2) the average (momentum)^ in the wave functions. If we then use the bare quark masses ~ 0, to^ ~ 150 



MeV and (p^) ~ (SOOMeV)^ then \5H\ ~ 7° ~ 50 MeV, and Eqs. (7.1) - (7.4) yield 5 ~ 1% 



For heavy-light quark combinations the difference in the kinetic parts of the Hamiltonian are larger. However, most 
of the masses are given in this case by the heavy quark mass itself, so that fractional deviations again remain small. 

Approximate mass equalities 2mii = ma + rriij or 2r nfj = to^^ + rnfj were suggested quite a while back on the basis 
of the SU(3) (Gell-Mann - Ne'eman) flavor symmetry |Q. The underlying assumption was that the flavor symmetry 
breaking part of the Hamiltonian behaves like H^, the neutral isoscalar member of an octet (a feature which is 
manifestly true for the QCD Hamiltonian), and that mass splittings can be obtained by using (0'^°^ with t/)'*'^ 
the SU(3) symmetric wave functions. The present discussion suggests that these relations can be transformed into 
inequalities for the linear mass combinations. 

When the particles in question are broad resonances, the mass inequalities need not apply. This is not just due to 
the technical difficulty of precisely defining, for example, rrip (different methods based on the tttt mass distributions 



®A. Falk tried using heavy quark symmetries to fix some of the unknown . The QCD inequalities nicely complement 
this program. 
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or the Argand diagrams for phase shifts yield somewhat different values for the masses of resonances [Q). Strictly 
speaking, the inequality 2m^5^^ > m^ul^ + mf]^ holds only for the lowest 1^ states in the M„s, M^g, M^j sectors 



which are KTr^inr, and (neglecting Zweig rule violations, i.e. mixing with the gluonic M'^ sector) KK P-wave states 
at threshold. The inequality would then appear to degenerate into the trivial statement on the kinematical threshold: 
rriK + < (l/2)(2TO7r + 2mK). As we will show in Sec. 16, the operator relation enables us to go beyond this and 
deduce relations for the phase shifts. To the extent that the threshold physics is completely dominated by narrow 
K* , p, and (f> resonances the mass inequality 2mK* > mp + can be regained from the phase shift inequalities. 

In the large Nc limit, resonance widths and Zweig rule violations v anis h like and l/N"^ respectively |67j and the 
distinction between the relation (5.2) and the approximate version (5.8) is lost. Also the flavor symmetry assumption 
may be on better footing as the gluonic degrees of freedom dominate. 

We proceed next to the baryon-meson inequalities, with the results listed in Tabl e pi[ The specific flavor-spin 
combinations appearing there are fixed in the manner explained in some detail in App. |D|.^| At the present the of 
A+ and S^, are not established and the choice made here (J^ — 1/2+ , rather than 3/2+) is essential for the inequalities 
to be satisfied. In general the baryon-meson inequalities are satisfied with a higher margin than the meson-meson or 
baryon-baryon inequalities. This could be attributed to the weak basis for the interflavor inequalities whose derivation 
required, beyond the rigorous operator relations, also an assumption of predominantly flavor symmetric ground state 
wave functions (or func tion als) . A lso the mesonic wave functions for the different flavors may indeed be very similar 
to each other [see Eqs. ( [7.l| ) - (7.4)], whereas the baryonic wave functions are intrinsically different from the mesonic 
ones. In the simple potential model picture the diquark subsystems are not at rest. Also the Y-configuration in the 
baryon's functionals are different from the configurations in the corresponding mesonic functionals. The baryonic 
ground state wave functions are therefore a poorer approximation for the mesonic wave functions than the wave 
functions of mesons with different flavors. 



Baryon-baryon inequalities are lis ted in Table [II. From Eqs. (5.11) we have also an inequality version of the equal 
spacing for baryons in the decuplet ( |3.11a ) and of the Gell-Mann - Okubo (GMO) mass formula for the baryon octet 
p^ . For the latter we use the inequality (5.11b): muds > (l/6)(mx;+ -I- m^- -I- nip + m„ -I- toh -|- m^o). Identifying 
uds with [1/ ^/2)[hP + TP)^ this becomes 3mA0 + Smj^o > (7715]+ + mj]- + niip + nin + + m^o), which in the limit 
where /-spin splittings are neglected becomes the GMO relation. While both the linearity in the decuplet and the 
GMO formula are very accurate it is gratifying that the small deviations are consistent with the inequalities. We 
should emphasize however that the flavor symmetry of the three quark baryonic ground state wave functions which 
underlies these inequalities is strictly an additional assumption. In spin independent quark-quark potentials ^(H), 
the inequalities hold for the fairly large class of potentials discussed in App. 

We have not listed most of the many mass relations involving radially excited states, Eq. ( 2.14 ). While we have 
several known radially excited states in the heavy quark QQ = cc, bb systems (particularly in the 1 channel), there 
are very few known radial excitations in the Qq or qq systems to allow useful comparisons. An exception is the case of 
radial excitations in the D (cq) system, for which the continuous experimental effort at Fermilab [|68[ keeps providing 
a relatively elaborate charmed meson spectrum. We have therefore indicated the relevant charm radial excitations in 
Table ffvl. 



TABLE I. Meson-meson Inequalities 





us sector: mus 


> 


2 (^^uu 


+ mss) 






pseudoscalar 


rriK > 1 (m^ + |m,, -I- |m^/) 






495.009 


> 


411.08 


vector 


mx* > ^ {mp + m^) 






893.14 


> 


894.7 


tensor 


niK; > \ (m„2 +mi^) 






1427.7 


> 


1296.55 


axial 


rriKi > \ {rriai + ruf^ ) 






1273 


> 


1256.0 


scalar 


niK* > \ (mao +mfg) 






1429 


> 


982 




uc sector: rriuc 
^ 1 / , \ 


> 


2 i^'^uu 


+ rricc) 







®The specification of the combinations of spins and flavors of the mesons appearing in these inequalities has used the notation 
and level assignments of the nonrelativistic quark model. This is, however, mainly done for convenience and does not detract 
from the degree of rigor of the derivation. Thus let us consider the comparison between the polarized A^"'" and p masses. 
Assume that the ground state has important L 7^ components. Then we could have a uu S = 1 diquark with L = 2, 

or 5 = with L = 1, . . . , and we would then simply need to construct, from the A^"*" wave functional, the corresponding trial 
wave functional for the lowest J = 1 ud state. The latter is still the p - no matter what orbital or other components it may 
have. 
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vector 
tensor 
axial 
scalar 



rriD' > \ {mp + mj/^) 
mo* > I (mas +"i-xc2) 
muo > 5 (mai + m^ei ) 

Xco) 



u6 sector: m^j > 



2008.9 
2458.9 
2422.2 

n 



> 
> 
> 
> 



1933.4 
2437.1 
2370 
2200.4 



pseudoscalar 
vector 
tensor 
axial 



TUB > I (m^ + mr) 
ms- > i (mp + mT) 
msj > i (ma2 + mxi,2 ) 
m^o > i (rua, +m^„ j 

'"^'o I ("'"1, + ) 



5279.0 
5324.9 
5739 

rugo 

"ID,. 



sc sector: m^c > i (rriss + nice) 



> 
> 

> 
> 



4799.2 
5115.2 
5615.65 
5560.95 

512i.() 



pseudoscalar 
vector 

tensor 
axial 

scalar 



niD, > ^ (f m^ + im^' + m^„) 
tud* > 5 (m^ + mj/y,) 

mo „+ > 5 (m/2 +mxc2) 



m,D^j 



> 7 (™ /i + m xci ) 



1968.5 
2112.4 

2573.5 
2535.35 



I (mss + m(,s) 



> 
> 
> 

> 



1832.0 
2058.15 
2415.6 
2396.2 

1709 



pseudoscalar 


mgO > i 


(§m 


+ ^rriri' + mr) 5369.3 


> 


5072.2 


vector 


ms. 




(m^ + mr) 


5416.3 


> 


5239.89 


tensor 


mB,+ 




fm/2 +mxi,J 


ms^+ 


> 


5594.1 


axial 






(m/i + m^.^ I 


ms^^ 


> 


5586.9 


scalar 


mB„+ 


>ll 


(m/o +m^,J 


ms„+ 


> 


5421.9 








c& sector: 


mc6 > 1 (mcc + m^i) 






pseudoscalar 


m^± 


>^ 


(m^^ + mx) 


6400 


> 


6219.6 


vector 


mB^^_ 


>i 


(mj/^ + mr) 




> 


6278.6 


tensor 






mxe2 + m;^,^ 


) "*Sc2+ 


> 


6734.7 


axial 


"*Bel + 




mx,, + mx,^ 


) + 


> 


6701.2 


scalar 








) "^Soo+ 


> 


6638.6 
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TABLE II. Baryon-meson Inequalities 



^ 3/2+ 



niA (uuu) > (3/2)mp 
TTiE (suu) > (l/2)(mp + 2mK' ) 
ms (ssu) > (l/2)(m^ + 2mK*) 
m^- (sss) > (3/2)m^ 
m„++(cuw) > (l/2)(mp + 2mi3*) 
m-He (csu) > (l/2)(mif* + + m^j) 
rriQO (ess) > (l/2)(m0 + 2mD*) 
(ecu) > (l/2)(mj/^ + 2mD') 
mj^+ (ccs) > (1/2) (mj/^ + 2mr,*) 
m^++(ccc) > (3/2)(mj/^) 
rriE^ (6uu) > (l/2)(mp + 2ms*) 
mH(, {bsu) > (l/2)(mK* + ms* + mg^) 
mn^ {hss) > (l/2)(m0 + 2mB* ) 
ms,,, [hbu) > (1/2) (mT +2ms.) 
»7in,6 (bbs) > (l/2)(mT +2msj) 
mo,,, (bbfc) > (3/2)(mT) 




1232.0 
1384.6 
1533.4 
1672.45 
2519.4 
ms, 

mE, 
ms, 
mn, 
"is,. 




IV IV IV IV IV IV IV IV IV IV IV IV IV IV IV IV 


1155 

1278 

1402.85 

1529.12 

2394 

2507.2 

2518.6 

3557.3 

3660.8 

4645.32 

5710 

5793.7 

5879.0 

10055.1 

10099.5 

14190.56 


J'^ = 1/2+ 


rriN (uud) > (3/4)(m-^ + mp) 
TTiE (suu) > ( 1 /4) (2mp + 3mi(: + thk* ) 
ms (uss) > (l/4)(2m0 + 3mx + rriK*) 
(cuu) > (l/4)(2mp + 3mo + mo') 
mn^ (ssc) > (l/4)(2m0 + Smn, + m^j) 
(ecu) > (l/4)(2mj/^ + Smu +mi3.) 
mf^+ (ccs) > (l/4)(2mj/^ +3mDs +m,i3j) 
msfc (buu) > (l/4)(2mp + Sms + rriB') 
rriA (uds) > (1/4) (2m^ + 3mK + mx* ) 
m.+ (udc) > (l/4)(2m.n- + Smo + mo*) 
m^o (udb) > (l/4)(2m^ + 3ms +ms*) 
ms^ (cus) > (1/4) (mx + mx*) + ("lu^ + m^) + (m^,* + mu*)] 
mst (bus) > (1/4) (mx + mx-) + {mB, + mu) + (m^j + ms*)] 
mf,cs (bcs) > (1/4) [(mu, + mo*) + (ms^ + tub) + {niB^, + ms*)] 




938.919 
1193.15 
1318.1 
2452.9 
2704 
ms, 

ms, 
1115.683 
2284.9 

5624 
2469.0 

ms, 

mtcs 




IV IV IV IV IV IV IV IV IV IV IV IV IV IV 


681 
980 

1104.25 

2288 

2514.1 

3451.5 

3552.9 

5675.5 

663.56 

1973.0 

5359.5 

2336.4 

5682.7 

6355.9 




TABLE III. Baryon-baryon 


Inequalities 










for the octet (Gell-Mann 




Okubo) 








3mA + 1 (mj;- + ms+ + m^o) > mp + m,i 


+ m=- + m=o 




4731.6 


> 




4514.1 


for the decuplet 


mE > (1/2) (mA + ms) 
ms > (1/2) (mn + ms) 






1384.6 
1533.4 


IV IV 




1382.7 
1528.5 


TABLE IV. Inequalities involving charm 


meson radial excitations 








pseudoscalar rn^^^ + m^' > ^ 


(0) , (1) , (0) 1 (1) 

mi,^ + rurfj + mir + mir 




mW 




> 


2138.2|^ 


vector m^'i + m^i > i 


{(>) 1 (1) , (0) , (1) 


4645.9 




> 


4508.95 



^°We can predict a lower bound on the mass of the as yet undiscovered D^^\ since the other masses in the inequality are 
known. 
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8. QCD INEQUALITIES FOR CORRELATION FUNCTIONS OF QUARK BILINEARS 



The set of all euclidean correlation functions, i.e. ordinary Green's functions continued to the euclidean domain, 
contains the complete information on any field theory |]69[| . In particular, two-point euclidean correlation functions are 
closely related to the spectrum of the theory, and indeed play a key role in most attempts to compute the spectrum 
of QCD. Following Weingarten 0, we will prove in this section inequalities for mesonic correlation functions. Let 



Faix,y) ^ {0\Ja{x)Ji{ym , 



(8.1) 



with Ja{x) a general, local, gauge invariant (i.e. color singlet) operator with the index a indicating Lorentz and/or 
flavor indices, be such a two-point function. Note that in the euclidean case all a; — y intervals are spacelike. The 
usual time ordering stating that the creation operator J\{x) should act prior to the annihilation Ja[x) is redundant. 

With an eye to the original Minkowski configuration, we will still assume that x^ —if^ > 0. In particular, we could use 
rotational (Lorentz) invariance to make x'^ — ip = \x — y\ = t by choosing y — (0, 0), x = {t, 0). Inserting a complete set 
of physical energy momentum eigenstates and using the (euclidean) time translation operator e~^* jQ(0)e'^' — Ja{t), 
we obtain a spectral function representation 170] 



Fa{x -y) ^ dfi aaifJ- )e 
for the correlation function. The spectral function is given by 



2^p-A'k-y| 



0'a(M^ 



■2) 



(8.3) 



with pn the four momentum of the state n. 

T he a symptotic behavior of Fa{x — y) (as \x — y\ 
Eq. (U): 



oo) is controlled by the state of lowest mass contributing in 



F^{x-y)^e-t'°\-y\ , 

where ~ means equality up to a residue factor 7q and powers of \x — y\~^ . 
two particle threshold. It could also be an isolated contribution p^(na) = 
numbers of the current J^, i.e. a state for which (riaj J^|0) 7^ 0. In this case. 



(8.4) 

Here /io could correspond to a physical 



,(0) 



Fa{x-y) 



>\x-y\ 



la 



of a particle with the quantum 



(8.5) 



We note that due to color confinement in QCD the finite energy physical spectrum consists of color singlet states 
and therefore it is sufficient to consider only correlation functions of color singlet operators.^ 

Eq. ( |8.5| ) is used in QCD lattice Monte Carlo calculations of the spectrum [^,^. An exponential form is fitted to 
the large (many lattice spacings) distance behavior of the appropriate correlation function. The latter is estimated 
by numerical sampling via the functional path integral expression for the correlation functions: 



, . wtr ^im ^ /d[A4:r)] /d[V;(x)] / d[^(x)] J,(:r) jt(x)e-^(-) 
^"^^ ^ /d[A,(x)] /d[^(x)] /d[^(a:)]e-^(-) 



with the action 



•S'qcd = j d^x £qcd ^ J 



J2MPA + nii)^,+tT{F;,^\ 



a\2 



i.6) 



■7) 



and pA ~ 7^(9^i + iA°yla^). For the purpose of lattice calculations, and also for the derivation of the QCD inequalities, 
we should eliminate the Grassman variables "0(3;) and (corresponding to the anti-commuting fields) in the 

functional integral. Since ipt and V'i appear only in the bilinear kinetic and mass term in the action, this integration 



We will not address the interesting suggestion that the inequalities mandate specific patterns of color symmetry breaking. 
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can be done in closed form Q. In the partition function in the denominator of Eq. (8.6) this results simply in the 
additional factor Y[i=li Det(9'A + mi), which modifies the integration weight from 

d^l{A)^e-''^'^'^^U[A^{x)]... , 
where Sym{A) = tr(F^j,A'')^ is the pure Yang-Mills part of the QCD action, to 

d^l{A) = e-«^*'(^) Yl Bet{pA + m,)d[A^{x)] . (8.8) 

The same change occurs in the numerator if Ja are functions of gluonic {Ffj_^) degrees of freedom only. 

The currents of interest - for the non-glueball part of the QCD spectrum - are, however, bilinears of quark fields, of 
the form = tp^Tipja- We use the scalar, pseudoscalar, vector, and pseudovector currents and bilinear expressions 
with extra D^j, derivatives: 

Jl = ^P^^ja. J%=i^tlBlt.^ja, (8.9) 

in order to extract the spectrum of the 0^,0^, 1+, 1^, and 2+ mesons made up of quark flavors qi and qj. 

In all of these cases we need to contract the extra (V')')/' at (x) and {'ip)ijj at (y). One chain of consecutive contractions 
can involve both ipix) and ipiv)- This gives a "connected" flavor structure with the quark qi propagating from y to 
X and antiquarks qj from cc to y [see Fig. ^(a)]; or, if i = j, we could also have the qi{x)qi{x) and qi{y)qi{y) in 
two separate contractions, yielding a (flavor) disconnected structure [see Fig. §(b)]. The expression for the two-point 
correlation function then becomes 

j A^i{A)tr \TS\{x,y)TS\{y,x)\ j J d^i{A) (8.10a) 

for the flavor connected case, or 

d^i{A) tr [rS\{x, x)] tr [TS'Av, v)] / / M^) (8.10b) 

for the fla vor di sconnected case. 

In Eqs. ( 8.10| ) SAx,y) is the full euclidean fermionic propagator (for flavor qi) in the background field A^{x); T is 



the Dirac matrix in the expression for the current ( p.9D ; and the trace refers to both spinor and color indices (which 
are suppressed). 

The quark propagator solves the equation 

{pA + m^)S\{x,y) =S{x-y) 

and formally is given by 

S\ix,y)^{x\{pA+m,)-'\y) (8.11) 

with^* ~ 1tj.Dfj,{A) involving the covariant derivatives Dfj_{A) = S^ + jAa^^ in the given background field configuration 
Afj_(x). The euclidean p is purely antihermitian and the 7^ are all hermitian and satisfy {7^,71/} — 2(5^^. 

It is well known that the singular behavior of tr[75S'^(x, y)] when x ^ y [ p9|j7^ , |7^ , or, alternatively some subtleties 
in the fermionic integration [ [74[ , may induce an anomaly ~ FF in the fiavor disconnected parts for the pseudoscalar 
case. Thus for much of the following discussion we will take i ^ j and avoid altogether the flavor disconnected 
contribution. 

An important element for the derivation of the inequalities is the positivity of the determinantal factor _]"[ Det ( / 
Da + mi), and hence also of the integration measure dfi{A) in Eq. (8^), for any given A^(x) configuration j|J^. The 



euclidean QCD action Ce = E'^ + B'^ is real so long as we do not have an iOE ■ B term, and hence exp [—Sym{A)] is 
positive. To show the positivity of the determinant factor, let us consider the eigenmodes tpA^^A) of the hermitian 
operator i Pa satisfying 

I PAi^^ ^ X{A)^P^ (8.12) 
with real eigenvalues X{A). The 75 anticommutation 757^75 = — 7^ implies the relation 



24 



l^A = -75 ^A75 • (8.13) 
This in turn forces all nonzero eigenvalues oil^A to appear in complex conjugate pairs. Indeed from ( S.12 ) and ( 8.13| ) 

« ^a(75^a) - lli pAlsi'x = -75(« Pa)^x = -X{A)-/5ipx , (8.14) 

so that "f5^x is the eigenfunction corresponding to the eigenvalues —X{A). Thus we have explicitly positive determi- 
nants 



Det(^A + m,) = J] [iX{A) + m,] [-zA(A) + m,] = ]J [X'^iA) + m^] 



(8.15) 



and hence 



d^j.{A) > 



(8.16) 



for all A^(x). 

It is important to notice that the above argument applies only for vectorial, non-chiral models with the quarks being 
Dirac fermions. If 75V'a is not a distinct new spinor (as is the case for chiral fermions) then the positivity argument 
breaks down. 

The measure positivity all ows us to prove inequalities between correlation functions if we can show that the corre- 
sponding integrands in Eq. (8.10a) satisfy, for any A^{x), the same inequality. Thus, let us assume that by simple 
algebra one can prove 

tY[raS\x,y)TaS^{y,x)]>tY[nS\x,y)rbS^{y,x)] . (8.17) 
The same inequality continues to hold after we integrate over the normalized, positive dfj,(A), yielding 

{0\Ja{x)Jl{y)\0) > {0\Mx)jl{y)\0) (8.18) 

with 

Ja{b){x) = '0i(a;)rQ(f,)-0i(x) . (8.19) 



We will next show that a relation like Eq. (8.17) holds with = 75. To this end we use Eq. (8.11) defining the 
propagator 6*^(0;,?/). We then have 



l5SAix,y)^5 = {xhsipA+m) ^j5\y)^{x\-{pA + m) '^\y) 



-nt 



= S\iy,x) 



(8.20) 



where the last dagger refers to the conjugate matrix in color and spinor space. 

It is instructive to see how the same result is obtained in various more explicit expressions for the propagator. Thus 
consider the contribution of any fermionic path on the lattice to the propagator S\{x^y) in the hopping parameter 
expansion ||7|]. It is [Q 



U{n, n + fi) [a + 67^) 



.21) 



where the product extends over any set of links connecting the initial point y to the final point x. Multiply- 
ing by 75 on the left and right and commuting the 75 factor through, we have Yin n+fi^i^^''^ + A)(o ~ ^7m) — 
On+zi n l^i^ + A' ^)('J — ^7ai)]^ (tli6 euclidean 7^ are hermitian), which is the same as the contribution of the reversed 
path of links (connecting x to y) to [5*^(2/, a;)] ^. Summing over all fermionic paths, we reconstruct S\{x,y) and 
S\{y,x) and Eq. ( 8.20| ) is satisfied. 
The key observation now is that if 



SA{x,y) = S\{x,y) 



.22) 



then Eq. (8.2C) implies that the integrand in the pseudoscalar correlation function is a positive definite sum of squares: 
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tr 



tr 



(8.23) 



with the sum extending over color (a) and spinor (r) indices. 

The integrand for any correlation functions of the other, non-pseudoscalar currents in Eq. ( ^.9| ) involves again sums 
of bilinear products of the same matrix elements {S'^)aa' ,tt' , but in general with alternating signs. Thus the integrand 
for the pseudoscalar correlation functions is larger than all the other integrands. This then yields the desired result: 



l^i(^)(j;Dny)|o)>(o|J^^(x)j^(2/)|o), 



(8.24) 



with Jr any non-pseudoscalar current. 



The requirement S\{x,y) = ^^^(x,?/), with i ^ j, is made in order to avoid the flavor disconnected contribution 



( B.10b| ) which could invalidate the derivation. It can be satisfied in the /-spin limit - with rriu — rud <C Aqcd 
and EM effects neglected - by taking i = u,j — d. Since S^{x,y) and S'^{x,y) then satisfy identical equations, 
^4 (x , y) = S'^{x,y) as required. The lowest mass particle in this i'ul^i'd channel is then the pion. The inequality 
(3.18) between the co rrel ation functions implies, by going to the asymptotic |a; — ?/| — > oo region and utilizing the 
asymptotic behavior ( ^.4|) , the reverse inequality between the lowest mass particles in the corresponding channels: 



m(0) < to[°) . 

Thus we obtain the result that the pion is the lightest meson in the ud channels 

TOtt < mp (1~) 



(8.25) 



.26a) 



< (1+), nia^ (2+) 



.26b) 



.26c) 



We note that these level orderings are indeed expected in any qq potential model |Q where the singlet S-wave state 
- i. e. the tt meson - is the lowest lying state. 

Could the various renormalizations which are required in order to render the correlation functions and the local 
products ^l)'^{x)T%l)^ {x) finite invalidate the derivation? The regularization required in order to make the path integral 
/ d[A^(x)] well-defined will not in general interfere with the derivation of the inequalities. Thus we could use a lattice 
discretization with any spacing a and with any action - the minimal Wilson action or with terms in the adjoint or 
higher representation. The inequalities will continue to be satisfied at every step of the procedure of letting a — > and 
V , the volume of the lattice, to infinity. Thus, to the extent that any sensible regularization exists, the inequalities 
will be satisfied. 

Instead of the local (non-pseudoscalar) current J^{x) — q{x)Tq{x)^ we could use a non-local gauge invariant version 



J^{x,x') = (7i(a;)rexpi f A^dx^qj{x') 



qi{x)TUcqj{x') 



'i.27) 



which creates an extended state - a quark at x, antiquark at x', and a connecting flux string along a curve c. The 
path-ordered Wilson line factor Uc is a unitary matrix. The corresponding correlation function 



J^{x,x') [J^{y,y')]^ 0)= f d[^l{A)]tr[TUcSA{x,y)UtrSA{y',x')] 



can readily be shown, using a Schwartz inequality, to satisfy: 

it 



J^{x,x')[j^{y,y')y O) < (0|J,(x)jt(y)|0)(0|Jp(x')j;(2/')|0), 
from which it follows that the mass of any particle x created by J^{x,x') still satisfies > mjr.f^ 



(8.28) 



(8.29) 



The advantage of the more general construction using a nonlocal is that x states of any spin can be created. 
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The derivation of the inequahties still applies regardless of the choice of "gauge fixing" which limits the allowed 
Afj^{x) (or Ufi,n+fi in the lattice version) which should be summed over, but does not affect the positivity of the measure. 
Also, the derivation is not sensitive to the issue of how we put th e ferm ions on the lattice. We can use Kogut-Susskind 



fermions ||76[ [which would correspond to a = 0, 6 = 1 in Eq. ( 8.21 )], Wilson fermions ||77| (a = 1,6 = 1), or any 
intermediate procedure.^ 

The inequalities are independent of the issue of quark and color confinement In fact the inequalities also apply 
in non-confining vectorial theories like QED and in the short distance \x — y\ —^ limit where perturbative QCD is 
adequate. They would amount to the statement that the pseudoscalar spin singlet "para-positronium" type state is 
the lightest (this happens in QED even with i = j = electron since the eflFects of the FF anomaly are weak for rrie ^ 0) . 
Likewise, the euclidean pseudoscalar correlation function dominates, as a; ^ (or g — > oo in momentum space), all 
the other correlation functions in the zero and one gluon approximation respectively. These inequalities first derived 
by Weingarten are the simplest of the exact QCD inequalities and yet, as will be made clear in the next section, 
are extremely useful. The positivity of the contribution of each A^{x) (or Uft.H+fi) configuration to the path integral 
defining the pseudoscalar propagator {0\Jps{x)J^g{y)\0) is a novel feature specific to QCD-type theories. It is quite 
distinct from the positivity (by unitarity) of the contribution of each physical intermediate state |(0| Ja|7i)p(5(p^j — /i^) 
to the spectral function which holds for an arbitrary current Ja in any theory, vectorial or otherwise. We cannot in 
general appeal to both types of positivity at the same time, since, in order to have the correct physical states, we 
have to sum over all the background An(x) configurations first. Otherwise, we do not respect even the translational 
invariance used in order to derive Eq. (^.31) . 

Anishetty and Wyler and Hsu [pl|| noticed that the measure positivity and ensuing inequalities hold also for 
SU(2) chiral gauge theories with an even number of fiavors (required in order to avoid global anomalies |^^. Indeed 
in this case, we can define 

V;' = + (^'+^//2)C, i = l,...,iV//2 
to be Nf/2 vectorial, interacting Dirac fields. 



9. QCD INEQUALITIES AND THE NON-BREAKING OF GLOBAL VECTORIAL SYMMETRIES 



The symmetry structure of a field theory, e.g. whether it spontaneously breaks vectorial and/or axial global 
symmetries present in the original Lagrangian, is closely tied in with the zero mass sector. Thus spontaneous breaking 
of global symmetries implies, via the Goldstone theorem |^3|, the existence of massless, scalar Goldstone bosons. 
Likewise spontaneous breaking of axial global symmetries implies Nambu- Goldstone massless pseudoscalars [ p4| . 
Finally, an unbroken axial symmetry may require for its realizati on m assless fermions in the physical spectrum (or a 
parity doubled spectrum). Evidently mass relations such as Eqs. (L9) between baryons (fermions) and scalar and/or 
pseudoscalar masses could restrict some of these possibilities and dictate the patterns of global symmetry realization 
in QCD and other vectorial theories. 

As a first illustration we will use the inequality (B.26c) - due to Weingarten - to motivate the Vafa-Witten theorem. 
The theorem states that in QCD (and in vectorial theories in general) global vectorial symmetries do not break down 
spontaneously. This will then be supplemented by some of the more rigorous discussion in the original Vafa-Witten 
paper. 

To be specific let us first consider a nonabelian global vectorial symmetry such as the SU(2) isospin symmetry for 
QCD. This symmetry is generated by 



Since the derivation of the inequalities may formally apply to any dimension and any shape of lattice, it should also hold 
for the new domain wall fermions [ |78[ . 

^''in the limit of unbound massive quarks we could have additional relations between the correlation functions. To leading 
order, the asymptotic behavior of all correlation functions is now the same: exp[—{mi + mj)\x — y\]. However, the P-wave 
intermediate states in the scalar or axial current correlation function cause a stronger (power) vanishing of f(/i^) at = mt+mj; 
hence the pseudovector correlation functions have an inferior power behavior |Q when \x — y\ ^ oo, as compared with the 
vectorial correlation functions. The inequalities {JY'^{x)JY'^{y)) > {JA^{x)J^{y)) and are true to all orders of perturbation 
theory and thus may be true even in the confining phase. In this case we expect additional mass ordering relations such as 
mp (1") < mai- 
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i+ = I d'x^Ux)M^) 

I' = \J d^a; [^i{x)Mx) - i^l{x)Mx)\ ■ (9-1) 

In the limit of ml'^-' ~ mf^ and no electromagnetic interactions it is an exact symmetry of Hqqd . If the symmetry 
breaks down spontaneously (either completely or to a U(l) subgroup generated by I3) we expect massless Goldstone 
scalars O^j and O^j. The essence of the argument against spontaneous breaking is that such massless scalars can be 

ruled out if m^"^ = m^^'^ ^ 0. Alternatively an exponential falloff as \x — y\ 00 of all correlation functions can be 

proven for ml°'' = m^^-* > 0. This conflicts with the power falloff expected due to the intermediate zero mass 0^^ state 



contributing via Eq. (B.5) to the spectral representation of the scalar correlation function. 



Let us present first the more heuristic argument which is analogous to the first argument of Vafa and Witten. The 



inequality (8.26c) m^^^ ^ > m^^^ ' states that the lowest scalar state in the ud channel must have a mass equal to or 



larger than the mass of the lowest ud pseudoscalar. If the /+ symmetry is spontaneously broken, then m ^ = 

ud 



and hence to' j ^ must vanish as well. If the axial global symmetry is not spontaneously broken, such a vanishing 



could only be accidental and hence most implausible. We next make a more precise and specific argument why to^j ^ 
should not vanish. 

In the limit ml"'' — rn '^^ = the QCD Lagrangian possesses an extra global SU(2) symmetry generated by the 
axial analogs of Eq. ( |9.l[ ): 

1+ = [ d^xijiix)j5Mx), etc. (9.2) 



The spontaneous breaking of could then yield a massless pseudoscalar ud Nambu- Goldstone boson and m|^j ' > 

m|^j ^ would then be trivially satisfied as > 0. However, following Vafa and Witten we keep ml"-' = m^^*^ ^ 0. No 

exact chiral symmetry then exists and to^^^^ ^ should not vanish. The Weingarten inequality implies m|^j ^ > m|^j ■* > 0, 
which as we argue next, negates the spontaneous breaking of symmetry. 

As we approch the mi°'' = m^^'' — limit, the lowest lying pseudoscalar particle becomes a pseudo-Goldstone 
particle whose (mass)^ is given by [ p9| , p5[ 

mlf^^(m<;^^+m^P)m (9.3) 

and vanishes linearly with the explicit chiral symmetry breaking (ml"'' + m^^'^) in the Lagrangian. Similar 
manipulations utilizing the analog of the soft pion theorem and current algebra yielding Eq. ( |9.3D suggest an analogous 
relation for the mass of the tentative O^t Goldstone boson y: 

ud 

If 7^ 0, i.e. the Goldstone 0+ boson does not decouple, then Eqs. (^^) and ( ^.4| ) are inconsistent with the 

QCD inequality "* > ''^Ij ^ long as m^'^ — rn'p <C rri"u^ + which we can maintain even as we approach 

(0) (0) n " 

TO,, — m. — (J. 



The inequality ( 8.26c ) was derived in Sec. §by using S\(x, y) = S\(x, y). If, however, /-spin {i.e. u ^ d) symmetry 
breaks spontaneously, this last equality could be violated as well! Is our argument then circular, and we proved /-spin 
symmetry only after assuming it at an intermediate stage'ij^ We do not think that this is the case. The rationale 



^^A similar argument was made recently in a more forceful manner in connection with the second Vafa- Witten theorem, 
regarding the non-breaking of the discrete parity symmetry in QCD. 
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behind our argument was alluded to in Sec. |^ above. Spontaneous symmetry breaking is a collective long-range 
phenomenon manifest in the V oo limit. 

The inequalities are proved first for finite lattices. The proof is valid for any volume V (and also for any UV 
momentum cutoff), no matter how large. We thus expect the inequalities to hold in the V oo limit. 

Vafa and Witten present also a more sophisticated and compelling argument. They show that the propagator 
S^{x,y) of a massive quark between two regions of size A and around x and y is bound for any Af^{x) by 



-mo|a;-iy|g2moA 

moA4 



(0) 



,(0) 



Too 



(9.5) 



The square of this bound applies to tr [r5^(a;, y, mQ)rS^{y, x, toq)] . Using the weighted normalized averaging with 
/ d/i(yl) ,d/i(A) > in Eq. ( 8.10a ), the same bound is also derived for the "smeared" correlation functions for any 
current: 



\JAix)my)]^0) < 



-2mo\x-y\p4m„A 



(9.6) 



The fini te sme arin g around x and y does not affect the asymptotic |a; — ^ A behavior, which is still given by 
Eq. (B.5). Eq. (9.6) then implies that the lowest state in any mesonic channel qiQj has a mass satisfying 



to|°^ > 2mo 



/ (0) ^ (0) 
(or ml-' > ml 



.(0) 



for 



,(0) 



(9.7) 



Thus for non-vanishing bare quark masses there are no massless bosons (and by a simple extension also no massless 



fermions) in the physical spectrum. In particular for m. 



(0) 



,(0) 



^ we could have no scalar Goldstone bosons. 



We now review the argument in detail. It is easy to show, first, that the propagator for a colored scalar in a 
background field is always maximized by the free propagator {x\KQ^\y) , Kq = —d^d^^ + m^ ||8^]. The latter can 
be cast in the form of a path integral expression |8q| with a positive definite integrand 



Do{x,y,mo) = (x 



1 



-(l/2)Tifo| 



dT 



x{T)=y 



dXf^{t) exp ■ 



x{0)=x 



dX 



At TO^T 



(9.8) 



The effect of any external gauge field is to introduce just the additional path-ordered "phase factor" for each path 
connecting x{Q) = x and xiT) — y: 





nx(T) 


P |exp 


I / dx^,{t)Al\a 
Jx{0) 



(9.9) 



This extra unitary matrix, when integrated over the positive measure in Eq. (9.8), decreases the functional integral 
/ d/i(A) . Consequent ly D A{x,y,mQ) < DQ{x,y,mo) (and in particular Z?yi(x, y, too) < exp [—TOo|a; — i/|]). 

The motion in Eq. ( |9.8[ ) is formally that of a particle in four space dimensions and a proper time. Vafa and Witten 
utilize an analogous expression for the fermionic propagator 



Si 



{a\iD- 



(9.10) 



where D is to be interpreted as the Dirac Hamiltonian in a (4-1-1) formulation, and a and f3 are smeared states of 
spatial extent A around x and y respectively. The norms {a\a) = (/3|/3) = A^* (norm / = max/ is required for the 
purpose of the subsequent discussion) are finite, and unitarity of exp [ir^iD)] implies 



(q;| exp [iT{iD)] 



< A" 



(9.11) 



The minimal distance between the supports of a and j3 is R = \x — y \ — 2A. Since the motion of the Dirac particle 
in 4-1 -1 M inkowski space is causal, the minimal "tim e" required for propagation bet wee n a and (3 is r^in — R- Using 
Eq. (9.11) and the lower limit Tmin = i? in Eq. (9.10) we obtain the desired bound (|9.6|) for the smeared propagator 



Si< 



dr e 



-mo\x-y\^2maA 
TOoA4 



(9.12) 
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The smearing of the states a and /3 can be made gauge invariant by a refined procedure (analogous to that used by 
Schwinger ||69[) without affecting the basic argument. 

For the abehan case the current, e.g. Jg = '4>i'y^'4'i carries no net baryon number. However, Vafa and Witten 
argue that if spontaneous baryon number violation is to occur we will have a vacuum expectation value of some baryon 
number carrying operator such as T{x) — Bijk(x) ~ tp^{x). This, in turn, would imply that {TA{x)TA{y)) cannot 
have an exponential falloff. Such an exponential falloff is implied by Eq. ( |9.5| ) and {T{x)T{y)) ~ J d/i(A) [5*^(2;,?;)]^. 

The need for smearing the fermionic propagator arises from the existence of zero modes of the Dirac operator Da 
p9| in any external B field Bz{x, y), with J B^ix, y) dxAy > 0o ~ [h/e). This lowest Landau level corresponds to the 
classical spiraling motion of particles along B with a Larmor radius r ~ B~^/^ . If B is aligned with {x — y) the particles 
originating in x do not diverge geometrically with a resultant \x — y\~'^ factor in the four-dimensional propagator. 
Rather, particles can be funneled from x to a spot of size ~ B~^ at y (see Fig.|lO|) yielding an enhanced propagator 
Sa — _Be^™°l^~^l (instead of Sa — e"'"°l^"^l/|a; — yp). The explicit B dependence prevents the derivation of the 



-mo|a:- 



^1 bound directly for 5*^1(2;, y) and requires the spreading over a region of geometric size A^. Once B 3> A 



further focusing of the motion to regions smaller than A^ will not enhance the a ^ (3 propagation and the uniform 
{B independent) bound Eq. ( |9.6| ) is obtained. 

The same type of spiraling motion also occurs for a charged scalar particle in a magnetic field. However, in 
this case the motion has no corresponding zero modes. The contribution of the spiraling trajectories to the path 
integral expression for the bosonic propagator DA{x,y) will decay as ~ e^^'^ at large proper times and is therefore 
unimportant, explai ning the fact that DA{x,y) can be bound without utilizing the smearing procedure. 

At first sight Eq. (9.7) seems somewhat surprising. In QED we have positronium states with nonzero binding and 
thus we might expect Eq. (9.7) to be violated in a weak coupling limit. However, recall that mo refers to the hare mass. 
If we have an electron and positron (or quark and antiquark), then besides the attractive interaction between and 
e~ there are the self-energies of e+ and Due to the vector nature of the gauge interaction this self-interaction of a 
smeared (or otherwise regularized) charge distribution is repulsive. Since 



(lrAr'pi{r)K{r — r')p2{r') < 



1 



(9.13) 



for any positive kernel K, the positive contribution exceeds the coulombic binding and Eq. ( p.7| ) does hold. Scalar 
self- interactions can be attractive and hence this reasoning is specific to gauge theories . Indeed for the scalar case 
"^0 ~* 5"/* + 'TiQ, and by choosing (j) ~ ~mQ/g, we can generate propagators D'^ which do not have the characteristic 
Q-mg\x-y\ bei^avior and the above proof of the Vafa- Witten theorem fails. For a 75 Yukawa coupling, toq ruQ + ig^^, 
and we even lose the positivity of the determinantal factor and of the measure. 



10. BARYON-MESON MASS INEQUALITIES FROM CORRELATION FUNCTIONS 

We next extend the discussion of Sec. p| to baryonic correlation functions of currents trilinear in quark fields: 

mM^)BlM\0) (10.1a) 



Bijk = ll^i{x)'lplix)lp'^{x)Tabc, 



(10.1b) 



with Tabc indicating some matrix in spinor space plus the eabc color factor. Eq. (q^) implies that the asymptotic 
behavior of F^i,{x,y) is prescribed by the lowest-lying baryon in this channel (say the nucleon for i^j = u,k~d 
and a coupling to overall spin 1/2): 



F^^A^,y)^e-^-\^-y\, |:,_y|^oo. 
Also in analogy with Eq. (B.lOa) we have a path integral representation of F^{x, y) 

Fi%ix,y) = / dp{A) SA{x,y)aa'S\{x,y)bb'SA{x,y)cc''['abcra'b', 



(10.2) 



(10.3) 



The strategy for deriving the inequalities is to compare this trilinear expression with the positive definite quadratic 
expression for the correlation function of the pseudoscalar in the rrii = rrij limit 
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If an inequality of the form 



i/p 



can be derived, with (1/p) > 1/2, then Eqs. ( |10.2[ ) and (B.5) imply the inequality 

TUN > {l/p)m^ . 



(10.4) 

(10.5) 
(10.6) 



[If (1/p) < 1/2 the last relation still allows for m,r > 2mAr, in which case the two nucleon threshold is the lowest 
physical state in the J^j ■* channel and the inequality derived from Eq. ( 10. 5| ) becomes the trivial uin > {'2/p)mN-] 
The baryonic correlation function in Eq. (10.5) is readily bound by 



i^ffe(a;,y) < / d^^{A) 



3/2 



In order to bound the right hand side by [i^j^(x, y)]^, Weingarten |^ appeals to the Holder inequality |9C 



< 



i/p 



|p/(p-i) 



[(p-i)/p] 



valid for p > I, and chooses 



(10.7) 



(10.8) 



\ aa' 

9 = 'Y\S'A{x,y)aa' 



(10.9) 



so that J d^fg is the right hand side in Eq. (10.7). For p ~ (n — 2)/(n — 3) the general Holder inequality (lO.S) 
implies 



d^,{A)Y,\S\{x,y)c 



dfi{A) 



(10.10) 



Note that the first term in curly brackets is simply F^j{x, y) which is then raised to a 1/p = (n — 3)/(n — 2) power. 
Let the number of ligh t degenerate flavors be larger than n which we take to be even. The second term in the square 
brackets in Eq. ( 10. IC ), in which the integrand is raised to the n/2 power, then has a direct physical interpretation. 
It is the two-point correlation function of products oi n/2 pseudoscalar currents with all flavors ii,im different: 



K{x,y) = (0 



t 









(10.11) 



There is no possibility of contracting a quark and an antiquark emanating from x (or from y), and we also have a 
unique pattern of contracting quark lines from x with those from y [see Fig. |ll|(a)] so as to form n/2 loops with 
separate color and spinor traces. 

Each of these traces then yields a factor X)aa' \S\{x,y)aa'\'^ and altogether we have (X)aa' y)aa' P)"^^ as 

required. If we now make the ra ther weak assumption that the lattice regul arized correlation function K{x^ y) is 
finite, we deduce from Eq. (10.10), that, up to a numerical constant, Eq. (10.5) is indeed valid. Finally, we conclude 
from Eq. ( |1G.6| ) that 



n — 3 

niN > nij, 

n — 2 



(10.12) 
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with n the number of hght degenerate quark flavors, which should be even. Thus we can obtain ttitv > (3/4)m7r 
if n = Nf = 6. This large number of light degenerate flavors required is associated with the particular derivation. 
Thus, if we assume only two light (u and d) flavors we can still construct Eq. (10.11) by taking all ii . . . i„/2 — u and 
Ji ■ • ■ Jn/2 — d and avoiding any flavor "disconnected" xx or yy contractions. We have, however, several contractions 

i- , -1 ^ n/2 

tr S\{x,y)SA{xTy) \ but also terms like tr{S\SAS \Sa) [see Fig. |ll|. 

We could improve the bound (10.12) if, instead of K{x, y) < constant, we appeal to the asymptotic behavior 



K(x,y)~e-'""%l"-^l 



(10.13) 



with TO|j*y2 lowest mass state created by the action of (</^j)"^^ on the vacuum. The bound derived from Eq. ( lO.lOj ) 



would then read 



,' n — 3 ^ 
mN > I rn„ 

n — 2 



n~2j "/2 



(10.14) 



If there were no bound states in the exotic channel with 7r+...7r+ quantum numbers then the lowest state is at 



threshold and 



{n/2)m.^. In this case Eq. (10.14) becomes 



TOjv > (3/2) m^r , 



(10.15) 



the result suggested by the more heuristic discussion in Sees. || and ||. To our knowledge the nonexistence of exotic 
7r+7r+-type bound states has not been proved. This can be done in the large limit (as will be discussed in Sec. p^ , 
in which case 



niM > 



(10.16) 



Weingarten Q makes the simple observation that if SA{x,y) can be uniformly bound by an A^{x) independent 
constant. 



SA{x,y) < constant , 



then from Eqs. (10.3) and (10.4) one can directly show that 

F^(x,y)<F^{x,y), 

implying 



(10.17) 



(10.18) 



(10.19) 



We can utilize for Eq. (10.17) the Vafa-Witten bound on SA{x,y) [Eq. (9.5)]. The requisite smearing of the points 
X and y into two regions of size A around x and y respectively was shown not to effect the asymptotic [a; — ?/| — > oo 

behavior and the ensuing mass relations. 

Weingarten [^ has independently motivated the bound ( 10.16 ) by considerations of lattice QCD: (-^)iatticc +™o — 
mo + R + il (with R,I being Hermitian, and R having a non-negative spectrum |9^). Thus for mo > 0, p + toq 
has no non-vanishing eigenvalues which implies that SAix, y) = {x\(Pa + fno)~^\y) is regular and bounded for all 
configurations. For mo = 0,1/ a does in general have zero modes which could be important for the issue of spontaneous 
breaking of chiral symmetry. 



11. MASS INEQUALITIES AND SxSB IN QCD AND VECTORIAL THEORIES 

The question of whether global axial symmetries are spontaneously broken in QCD (or in other field theories) is of 
utmost importance. In the ml°^ = m)^"^ = hmit, QCD has the global axial SU(2) symmetry generated by 1^,1^, 
and /|. If we can show that this symmetry is spontaneously broken {i.e. S^SB occurs), then QCD is quaranteed 
to reproduce Goldstone pions and the successful rich phenomenology of soft pion theorems and current algebra |5|] 
developed in the late sixties. 
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Conventionally, SxSB is achieved via formation of a qq condensate in the QCD vacuum, (ip-ip) 0. (ipip) is the 
order parameter for a phase transition from an axial SU(2) symmetric phase at high temperature to the broken phase 
at low temperature f^^-J";/^) (at zero temperature) serves, along with (F^), as one of the nonperturbative inputs 
in the QCD sum rules |9^,Q. Many attempts have been made to prove {ipi') 7^ in lattice QCD, to estimate its 
magnitude, and its finite temperature behavior. Evidently 



(11.1) 



and {ipip)A was shown by Banks and Casher to equal 7rp(0), with p(A) the spectral (eigenvalue) density of the 
Dirac operator i PaiI-^a ~ ^''Pa- These authors also argued that the same set of "fluxon" configurations of Afj_{x) are 
responsible for both confinement and the requisite dense set of zeroes of Pa- This supports earlier, more heuristic 
arguments 
to SxSB. 

The issue of S^SB is also crucial for the case of composite models for quarks and leptons |)7[ - |99f| ■ In such theories 
one assumes that the axial symmetry is not broken. The existence of practically massless fermions (me, mu.md — & 
few MeV and even mf, — 5 GeV are very small when compared with the compositeness scale Ap > TeV [^,100|) in 



96 1 suggesting that confinement, or binding, of a massless fermion in a vectorial theory necessarily leads 



the physical spectrum is then believed to be a manifestation of these axial symmetries. However, the baryon-meson 
mass inequalities such as mAr > m^r provide an additional strong argument for SxSB in QCD and similar vectorial 
theories. 



Let us assume that the axial isospin symmetry (which holds in the limit m 



(0) 



,(0) 



0) is not broken sponta- 



neously, i.e. /^|0) = 0. This symmetry can then be realized linearly in the massive spectrum by having degenerate 
parity doublets. However, one can show via the 't Hooft anomaly matching condition jl^ that we must also have at 

least one 7=1/2 massless physical spin-1/2 fermion (the nucl eon) . 

Anomaly constraints have been discussed extensively ]l^ , |l3 ,101|. Following Ref. [101|, we consider the vertex 



r^.A (919293) = j A^xi J d4x2e^('?^"^+«^"^)(0|r[J^(xi)J.(x2)JA(a;3)] |0) (11.2) 
of three identical fiavor currents, which, in terms of chiral right and left combinations are 

J"" = C7'' (1 + 75) + (1 - 75)] i^aj , (11.3) 

with A± matrices in flavor space. We take tiA^ = trA_ = to avoid the 75 anomaly involving SU(Nc) gluons 

df.J'' = acitrA+-tTA-)tTFF = . (11.4) 



It was shown |[72|,[73 , 102 1 that Tfj^^x satisfies the anomalous Ward identity 



93%.A = ^tr(A^-A^_) 



(11.5) 



The right hand side is given by the lowest triangular graph involving the massless fundamental fermions of the theory. 
Equation (11.5) is true to all orders |101 lOS] and most likely is a genuine nonperturbative result. Equation (11.5) 



implies that some invariant amplitudes in the decomposition of T^^\ have a pole at (93) = of prescribed residue. 



given by "the anomaly", i.e. the right hand side of Eq. (11.5). Such a singularity can arise either by 



(a) having zero mass (pseudo-) scalar boson states with appropriate nonvanishing (n*^! J+|0) matrix elements; 
or by 

(b) having a multiplet of massless spin 1/2 physical fermions ("nucleons") such that in this basis 

J" = Viiv7^ [S+fe)(l + 75) + - 75)] i^^M . (11-6) 

with a "matched anomaly" 

ti{Bl - B^_) = tiiAl - At) . (11.7) 
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Possibility (a) implies Q|0) ^ with Q the charge J d^a:Jo(a;), i.e. a spontaneous breaking, which we ass umed does 
not occur. This leaves us then with possibility (b), i.e. the need to have massless fermions satisfying Eq. ( 11.7 ). The 
last algebraic anomaly condition constrains the various composite models of quarks and leptons where the latter are 
the massless physical states iV^ [|l2|j99[ . 

Let us now appeal to a new element, namely to the nuclcon-pion (or fermion-boson) mass inequalities 

mAT > (jnp > "t-b) ■ (11-8) 

Possibility (b) of massless fermions can then be ruled out. Eq. (11.8) implies that if the nucleons are massless so 
are the pions. However, unless the matrix element (0|J^|7r) ~ j-nl^i accidentally vanishes and the pion "decouples", 
we then have spontaneous breaking of the axial symmetry which is precisely the possibility that we were trying to 
exclude. 



It has been shown |104| that in the ~^ oo limit, that anomaly matchin g via possibility (b) is ruled out. Amusingly 
in this limit we have, as indicated in Sec. O below, the inequality ( |l0.16[ ): M^r > (A^c/2)tot- It excludes, for infinite 
iVc (and some very small non-vanishing rriu and rn'^'' and consequently finite m^r) even finite baryon masses, so that 
alternative (b) cannot a fortiori be realized. 



symmetry via niu 



(0) 



,(0) 



It is interesting to see how the inequality (ll.S) is satisfied if we introduce a small explicit breaking of the axial 



7^ 0. In this case we have no strict arguments for the vanishing of either mjv and/or m^, 



From Eq. (9.3), we have ttItt ~ VmP. If the nucleon mass is generated exclusively via the explicit chiral symmetry 
breaking term mPijjqipq, then this term and mNipNip^ both have to flip sign under the discrete exp[i7rQ5] transfor- 
mation. This suggests that — am^ < m.,^ ~ ^/m° for small m°, violating the inequality > rn,^. This in turn 
rules out the possibility of explicit chiral symmetry breaking only and SxSB is again indicated. 



An alternate approach to proving SxSB [lOS] would be to start from the Weingarten inequality ( B.24 ) for the 
pseudoscalar and scalar two-point correlation functions: 



and show that the inequality is strict 

(0|J^i(cc)j;i(y)|0) > (0| (a;)jf (y)|0) . 
If chiral symmetry remains unbroken (QsjO) — 0), we can show, by utilizing 



(11.9) 
(11.10) 



that 



\J%{x)J%{ym = (0|jf (a;)jf (y)|0) 



(11.11) 



Thus, proving Eq. (11.10) would amount to proving SxSB. 

The fermionic propagator in a background field can in general be decomposed in terms of the Dirac 7 matrices: 



Sf{x, y\ A) = S{x, y- A)I + -f^,v^'{x, y; A) + af.^f'^x, y; A) 
+ lbliia^{x, y; A) + -i^p{x, y; A) . 



(11.12) 



Substituting this last expression into Eq. ( ^.lOa ) with P = 75 and P = /, and utilizing Eq. ( |8.2C| ), we find that v^^ 
and contribute equally to (JpsJps) and {JgJs)- However, s,t^'^, and p contribute |sp + |pp to the first 

IpP to the second. Thus, if we can pinpoint any set of gauge field configurations for 



correlation and —\s\''^ — {f^ 
which 



!/|2 



dM(A)(|s|^ + rn' + ipr)>o 



(11.13) 



then those configurations could generate SxSB by making the strict inequality hold. Note that the positivity of 
expression (11.13) and of d^{A) insures that the effect of such configurations cannot be cancelled by some other 



configurations. 

The topologically invariant relation 



d*xp{x,x;A) — / dL^xtv [S(x, x; A)j^] 



mo 



d'^x Ff_,^{x)Ff^^{x) 



(11.14) 
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suggests that if we have an instanton |106| - anti-instanton [107| gas, the contribution of the region near any pseu- 
doparticle is 



J d'^xp{x, x; A) ~ ±1, and J dfi{A) \p{x, x;A)\'^ >0 . 



(11.15) 



Indeed, it has been speculated by s everal authors [107, 108| that in stantons may be the source of SxSB. The detailed 



109 lie]. 



quark propagator structure in Eq. (11.12) was utilized in 

One may try showing directly the existence of a zero mass state in the QCD spectrum in the m 
Thus, if the representation for the pseudoscalar propagator 

(^pf(a:)j;/(y)> = / d^iiA)tT\s\ix,y)SA{x,y) 



(0) 



,(0) 



hmit. 



as a sum of positive definite contributions could be used to show a power falloff of the correlation function, then the 
existence of massless hadronic states with 0~ quantum numbers (massless pions or massless nucleons) would follow. 

The above program attempts to achieve the opposite goal as compared with the original work of Vafa and Witten 
[0 described in Sec. ^ where upper bounds on the euclidean correlation function were found. It has not been realized 
in the form described above. Vafa and Witten have, however, succeeded by combining the measure positivity and 
various index theorems to prove a closely related result which we will briefly describe next. They considered the 
fc-point function 



S{xi, . . . ,Xfc) = {0\lpl1p2ixi)lp2'(p3{x2) ■ ■ .1pk1pl{xk)\0) 



(11.16) 



With all fermionic flavors distinct, there is a unique, cyclic pattern of contractions. The fc-point function then has 
the path integral representation 

(11.17) 



S(xi,...,xk) = J dfi{A) tT [Sa{xi, X2) ■ ■ ■ SA{xk, xi)] 
Integrating over all Xi and using translational invariance defines 

S{k) ^ J d'^xi . . . J d'^Xk-i S{xi,X2, ■ 



,Xk-i,Xk = 0) 



d'^Xk S{xi,X2) . . ■ S(xk-i,Xk) 



The idea is to prove that (after appropriate regularization) 

S{k) > Cfc5o(fc) 

with 



5o(fc) = Ad{R) 



d^p 
(27r)4 \p 



(k/2) 



(11.18) 



(11.19) 



(11.20) 



[here d{R) is the dimension of the (color) representation of the fermions], the corresponding expression for the free 
massless fermions. This implies that for fc > 4, S(k) has, like So(k), infrared divergences which in terms of the 
physi cal st ates can be understood only if we had zero mass hadronic states. S{k) can be written, using Eqs. ( |11.17| ) 
and (11.18), as 



S{k)^ / d^i{A)SA{k) 



(11.21a) 



SA{k) 



V 



d'^xi ... d'^Xk tr [Sa{xi,X2) . . . SA{xk,xi)] 



V 



1 



(11.21b) 



wher e the l ast trace and operator multiplication refer not only to color and spinor space but x space as well. To prove 
Eq. ( 11.19D , it is clearly sufficient to show, due to the measure positivity, that 
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SA{k) > CkSoik) 



(11.22) 



holds for each configuration separately. SA{k) can be expressed in terms of a Dirichlet sum over the eigenvalues 
of Da: 

iDA^ = M (11.23) 



SA{k) 



- OO 



V 



(11.24) 



i=l 



with Ai, A2, . . . , A^v the eigenvalues in ascending order, 
scales with the size of the box considered as 

C 



Vafa and Witten then prove that the lowest eigenvalue Ai 



1j 



(11.25) 



with C independent of the gauge field. This implies, for even k, that S{k) > C~^V~^ and the required infrared 
divergence ensues when k > A (the li miting case of A: = 4 can also be handled by a more refined discussion). We 
will not present here the proof of Eq. (11.25) which involves topological considerations and tracing out the flow of 
eigenvalues under gauge deformations and refer the reader to the original paper 

All the above disc ussion, while strongly suggesting SxSB in QCD, fails to demonstrate that is indeed non-zero. 
Recently Stern 1 1 1 1 ] suggested a novel pattern of SxSB in QCD with the pseudoscalars still being the Nambu- Golds tone 
bosons associated with this spontaneous breaking, but where {ij^ip) = 0. Indeed, as Kogan, Kovner, and Shifman |112] 
noted, there could be some residual "custodial" discrete axial symmetry which allows only higher (-0?/')^ order 
parameters to have non-vanishing VEVs. The new scheme is phenomenologically interesting. In particular, since now 



(0) ^(0) 



and a more symmetric rriu 



(0) (0) (0) 



4 = ©(m^), larger values of m„ 
However, as pointed out by Kogan, Kovner, and Shifman |112], the inequality 



mass pattern would be implied. 



Ca{x) 



JlAx)J^.A{m) < (0 1 4, (x) J,, (0)1) ^ Cp,{x) 



(11.26) 



can be judiciously utilized to rule out this pattern. Both correlators have their asymptotic behavior controlled by the 
lightest 0^ pion states, i.e. Cyi(a;), Cps{x) ~ e"™"'^' for x 00. However, the vacuum-to-pion axial matrix element 
has the conventional form 



(01410 

with _FV scaling like Aqcd and not vanishing in the niq 
element behaves in this scheme as 



ab 



limit; whereas the vacuum-to-pion pseudoscalar matrix 



If M ~ ©(Aqcd) is the mass of the first massive (non-pion) state in the channel, then for m^^ ^ x ^ M ^, we 
have [recalling that the pion intermediate state still dominates Cps{x)] 



Cps {x) 



-m,^ \x\ 



and 



Ca{x) 



and the inequality (11.26) is violated 



The analog of Eqr(|TT 



holds also for composite models based on an underlying vectorial gauge interaction [[18[ . 
This interaction should confine the massless hyperons at a scale (Ap)~^, leaving us with the quarks, leptons, Higgs 
particl es, an d conceivably also VF*, (and even the photon in some models!) as the low lying physical states. 

Eq. (11.8) suggests that global axial symmetry, assumed to protect the ( almost!) massless fermions from acquiring 
large masses ~ Ap, does spontaneously break down. The analog of ( ll.S| ) here is > raH and the fact that we 
have no Higgs particles lighter than the lightest quarks and leptons is bothersome. (A scenario with only three 
Higgs states, which disappeared in the process of SU(2) x U(l) breaking, appears to conflict the existence of several 
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fermionic generations.) Finally, we note that the m^^ ^ > m^^ ^ inequality suggests that any composite v ector particle 



(0-) 



should be massive, and gauge symmetries, and in particular the exact U(1)em (with < 10 eV! [ 113 |) should 
not arise "accidentally" due to the existence of almost zero mass vectorial bound states. 



All the constraints stemming from mass inequalities do not apply to chiral preonic theories [ 114 - 116 |, and/or 
prconic theories with fcrmions and bosons with Yukawa couplings |118]. In this case the measure positivity d/i(yl) is 
lost (inequalities cannot be proven in supersymmetric models). Also in chiral models the pseudoscalar current "075 V' 
which was used extensively above vanishes identically. The mass inequalities, together with the study of the anomaly 
constraints, shifted the focus of research for c omposite models from the early work on vectorial models ||l8,117| to 
chiral and/or fermion-boson composite models |114,116,119 120 suggested earlier Indeed it was observed [122] that 
in gauge and scalar (not pseudoscalar!) field theories, one can then prove the inequality 



(O|0^(x)0(x)0t(o)0(O)|O)(O|^(x)75V(a;)V^(O)75V'(O)|O) > |(O|0^(x)0(x)0^(O)0(O)|O)| 



(11.27) 



by integrating over ip, -0. Using the positivity of the d[A^]d[(/)]e"'^["^f"'^] integration and of the pseudoscalar correlator, 
the desired relation follows readily as a Schwartz-type inequality. However, we cannot now infer the mass inequality 
W-F > (l/2)(r?T,7r + ^TT-x) between the lightest particles in the (j)^ip, i/j-fs'il', and channels, simply because the 
{(f)^ (f){x)(f)^ (f)(0)) correlator always has the constant ((O|0^(/)|O))^ contribution due to the vacuum state, which cannot 
be avoided without spoiling the positivity. This negates the claim made by Nussinov in |l9| . 

Hsu suggested that QCD, and in particular the Vafa-Witten inequalities along with the measure positivity for 
SU(2) chiral with an even number of flavors, ca n be use d to exclude the possibility that such strongly interacting 
chiral theories underlie the standard EW model [124 jl25 |. Since his analysis relies on Vafa-Witten upper bounds for 
both fermionic and scalar propagators, one needs to choose a regularization point where the X(f>^ coupling vanishes 
(otherwise the quadratic integral defining ^^[^^(a:, 0)] cannot be performed). 



For a long period it was not clear, in view of difficulties encountered in putting the theory on the lattice | 126 ], if 
cons istent chiral gauge theories could be defined. However, the recent domain wall fermions [T^ and overlap formalism 
[127 1 put the (lattice) regularization of fermionic theories with gauge interactions on a much more solid foundation. 
These de velopments also allowed for lattice calculations incorporatin g ch iral symmetries in a more explicit manner 
[128,129|. Most likely proper regularization of chiral gauge theories |13C] will also soon be feasible. Our preceding 
discussion then suggests that chiral composite models for quarks and leptons may revive. 

We will not pursue here the possibility that this novel approach could "rigorize" the derivation of (some) of the 
QCD inequalities. 



12. INEQUALITIES BETWEEN MASSES OF PSEUDOSCALAR MESONS 



In this section we conclude the discussion of mass inequalities in the non-exotic qqq (baryonic) and qq (mesonic) 
channels, by proving that pO[ |: 



1 

> - 
- 2 



(0" 



(12.1a) 



with 



the mass of the lowest lying pseudoscalar meson in the qiqj channel; and JlOyi 



(12.1b) 



The derivation of both inequalities relies on the positivity of the integrand in the functional integral expressing 
correlation functions of pseudoscalar currents. We have, however, to make the additional assumption (a) that the 
qq annihilation ch annels generating the flavor disconnected part [Eq. (fe.lObl)] make negligible contributions. For the 

7r / 



derivation of Eq. ( 12.1b ) we could assume instead (b) the validity of the soft pion expression [131| for ( 



in terms of the (axial) vector spectral functions 
functions 



0|. In order to derive Eq. ( 12.1a| ) we compare the three correlation 



IJ 



\ix) [4(2/)]^|0) = I d^,iA)tr{[S\ix,y)ySiix,y)} 



(12.2a) 



^''Aharom; 
(see App. p 



et al. 121 1 foun d that QCD-type inequalities were still useful in SUSY theories, though in the mnig 



oo limit 



Nishino [L22| has shown, using the Vafa-Witten inequalities, that in SUSY theories parity is conserved. 
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(0|JS(^)[JS(y)f |0)= / d^iA)tr{[S^A^,y)]^S^A^,y)} 



(12.2b) 



(12.2c) 



where we have left out, following our assu mption (a), the flavor discon nected con tribution, e.g. 
J dAi(A)tr [756'^ (x,x)] tr [j5S\{y,y)] in Eq. ( |l2.2bD . The expressions in Eqs. ( |l2.2bD and ( |l2.2cD have the form 
of perfect squares of vectors S\^^ (or S-^^^) with A^(x) and a, a the spinor, color indices viewed as a joint index. 
Equation ( 12.2a ) is the corresponding scalar product S]^^^ ■ S\^^. The Schwartz inequality {S^)'^{S^)^ > \S^ ■ S^\'^ 



therefore implies that 



[JUy)]^\0) < {0\JU^) [J!M]^ \0){0\/^,{x) [J^,{y)]^ |0) 



and Eq. (12.1a) immediately follows from the last inequality by using Eq. (8.5) and going to the \x — y\ 
where 



(12.3) 
00 limit. 



(0^)1 I 

^~ml '\x-y\ 



(12.4) 



Eq. ( 12.3 ) is the lowest in a hierarchy of inequalities stating that all the principle minors of the matrix My = 5'* • 
are positive. However, as can be readily verified by considering the 3x3 minors, no useful new mass inequalities 
follow. 

We have noted in our earlier more heuristic discussion of the inequalities the need for making a "Zweig rule" 
assumption (a). It amounts to suppressing the annihilation of qiQi into gluons and allows us to consider the mesonic 
sector rriij = qiqi to be distinct from m^^\ the flavor vacuum sector. This is valid particularly in the heavy flavor QQ 
sector. 



(0) , (0) 

TO,- + rrij 



The inequalities are also automatically satisfied in a soft pion limit with 

We recall that in the actual comparison with particle data in Sec. the pseudoscalar inequalities were indeed satisfied 
with a larger margin t han th e other inequalities. 

We next prove Eq. (12.1b) along similar lines p^. We compare the propagators 



(J^A^KAy))^ f d^^{A)tr{[Sl{x,y)]^Si{x,y)} 



and 



with 



iWJ^oiy)) = \j di,{A)i,[[Sl{x,y)]^ Sl{x,y)+[S'i{x,y)\^ S'^ix^y)} 



(12.5a) 



(12.5b) 



(12.6a) 



Jla = ■y|(V'n75V'« - i^dlS'fpd), (12.6b) 

i.e. the pseudoscalar currents with the quantum numbers of the tt"^, tt" mesons respectively. We can take into account 
electromagnetic effects, which are the source of the 7r+ — 7r° mass difference]^, by considering the gauge group to be 
SU(3)cX U(1)em- The path integral measure then becomes 



It was realized early on that the Cottingham formula expressing the EM contributions to the mass differences, converges to 



the correct vr — vr mass difference 
the dominant contribution. 



132 



133], and fails for the K+ - K" and n — p difference - for which m 



(0) 



mif ' makes 
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N, 



,EMip-SYM(A;:)g-SEMK") 



and the propagator in the joint external fields is 



(12.7) 



(12.. 



Since the EM interaction is vectorial, the arguments leading to the measure positivity, d/i(A) > 0, and t he po sitive 
definit e form of the integrand for the pseudoscalar correlation functions remains intact. Comparing Eqs. (12.5a) and 
( 12.5b| ) we can then show, from + 6^ > 2a • b, that 



(12.9) 



from which > m^o readily follows. 

Th e starting point for the original derivation of > m^o by Witten jl^ is the soft pion - current algebra relation 

P3 



7r~ 



(12.10) 



expressing the tt"*" — tt^ mass difference as an integral over the momentum space vector and axial vector (isovector) 
correlation functio ns. T he latter are Fourier transforms of the euclidean space correlation functions which can be 
expressed via Eq. ( B.lOa ), so that we have 

{V'{k)V^{-k)) ^ {Al{k)Al{-k)) 



d ye 



-iky 



2 
V 



d ye 



-iky 



X Jdn{A)tr[yf,SAix,y)jf,-jf,-f5SAix,y)jf,-f5]SA{y,x). (12-11) 

To this lowest, oem order, there are no flavor disconnected contributions to the isovector correlation function. The 
factor of 2 comes from = ■07^T^?/) with tr [(T^)^] = 2. 

We note that in the present case - unlike in all previous inequalities - the actual sign of the expression in (12.11) is 
of crucial importance. A minus sign coming from the fcrmion loop has been cancelled by the fact that the euclidean 
currents are — i-ip^^ip (and A^ ~ *'07;i75'0) ^^id the supplies the extra minus sign. [7*^ — > ijjf is essential since 
the euclidean 7^ are hermitian and satisfy {7^,7,f} = S^u, whereas ^ jff, 7, ^1 = 

The difference [7^S'a(x, y)jfj, - 7^755'yi(a;, 2/)7^75] occurring in Eq. ( |l'2.1l| ) is simply ^^Ea{x, y)-f^,, with Ea{x, y) = 
(a; |mo(— I?^ + mp)"^! y), the 75 even part of SA{x,y) = (x\{Da + mo)^^ \ y) — EA{x,y) + OA(x, y). T he odd 
OA{x,y) behaves like a product of an odd number of 7 matrices. Thus tr(7p£^^7^0yi) = 0, and Eq. (12.11) can be 
rewritten as 



{V'{k)v;!{-k)) {Al{k)Al{-k)) = 



V 



d^xe''^^ / d^e 



-iky 



X / <in{A)tr[-f^EA{x,y)'yf,EA{y,x)] 



(12.12) 



Viewing = e*'^^7^ = e''^^(7^)cQ/ as a matrix in spinor {aa') and coordinate space jointly, we can rewrite this last 
expression as 



- {Al{k)Ali~k)) - - / dM(A)tr [M,EaM;Ea] , 



(12.13) 



where the last trace and matrix multiplications refer also to the coordinates x,y as indices. 
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The operator Ea is positive definite, as can be seen by going to the basis defined by zD^jA^) = A^IAyi). In this 
basis 



tT[M^EAM;EA] ^ I(M^)AaAa' 



is manifestly positive. The measure positivity then impHes that / d/i(A) tr [M^£'^M*i?^] is also positive (non- 
negative). We finally arrive at 



{V^{k)V^{-k))-{Al{k)Al{-k))>0. 



(12.14) 



and m7r+ > WttO follows from Eq. ( 12.10|) 



It is amusing to note that the Eq. (12.14) complements the asymptotic chiral symmetry |134|| (V'.f (fc)V^(— fc)) 



(A^(fc)A^(— fc)), as fc ^ oo. Such an asymptotic equality is indeed requi red for the d^k integration in Eq. ( 12.10| ) 
to converge. It motivated the spectral function sum rules of Weinberg |135| som e time ago. Using unsubtracted 
dispersion relations for {V^ {k)V^ {— k)) and (A^(fc)yl^(— fc)), the inequalities ( |12.14 ) imply 



> 



d^2 aA{fi^){fi^ + k^y 



with (Tv'((Tyi) the vector (axial vector) spectral functions. 
The configuration space analog of the inequality 



{A^{x)A^m < {V^{x)VM 



(12.15) 



(12.16) 



has been discussed at length |136[ but no convincing proof exists at present. One difficulty is that the one pion 
contribution to the "longitudinal" part of the axial correlator Ca{x) ~ e"™"'^' dominates at large distances and an 
appropriate transverse projection of {Af^Ai^) is required. (This is not the case in momentum space since the large 
pion pole contribution at small k is suppressed by fc^fc'^ "derivative coupling" factors.) 

Efforts to prove such an inequality are motivated not only by the experimental fact that the lightest hadr ons in th e 
1~ , 1^ sectors satisfy > nip. Precision tests of the weak interactions indicate that the "5 parameter" 137.13^ 



S^ — [Fv{k^)^FA{k^ 



(12.17) 



with Fv , Fa the covariant (transverse) parts of {A^{k)A^{—k)) and (y^{k)V^{—k)) , is negative. The second moment 
of the c onjectur ed inequality ( 12.16 ) will imply that S' > in all vectorial theories such as the original technicolor 



models |T3|JT4 



thereby ruling those out as viable mechanisms for dynamical breaking of the EW symmetry of the 

Standard^odel. 

The previous derivation of > m^o, based on Eq. (12.9), did not utilize the soft pion limit (tantamount to 
letting TTiu'^ + rri"p ~* 0), but only the weaker assumption of isospin \m^^ — — )■ 0. [In particular isospin was 

utilized to argue that there are no disconnected purely gluonic intermediate state contributions to {0\Jp'-'^° Uy)\n) .] 

At first sight the previous derivation appears to lead to a stronger result. This however is not the case [ p_41[ . The 
point is that there are intermediate states with one photon which make 



{0\{ipul5'4'u - V'd75'/'<i)|muhigluon state + 7) 7^ . 



(12.18) 



and which contribute to the same order ©(ckem)- (Because of charge conjugation parity and color neutrality, the 
lowest perturbative state of this type consists of at least three gluons and a photon.) This contribution could be 
neglected if we appeal to the analog of the Zweig rule hypothesis (a) used above. Alternatively, we could go to the 
soft 7r° limit in which case all t he cou plings of the Goldstone particle to the photon and any neutral hadronic system 
[the "multigluon state" in Eq. (12.18)] must vanish. 



Since the EM interaction conserves J^^ , soft pion theorems imply that the EM contributions to a massless neutral 
Goldstone particle vanish to all orders. Thus, if m^o(a = 0) = then also m^o(a ^ 0) = 0. The result to^+ > m^o, 
and its analog in other vectorial gauge theories, is therefore of great importance. Otherwise < to^o = and the 
TT"*" becomes tachyonic. A condensate of charged pions could then form, leading potentially to spontaneous violations 

(T) (T) 

of EM charge conser vation. S imilar results like m + > 771 could help fix th e pa ttern of symmetry breaking or 
"vacuum alignment" [142, 143 1 in other vectorial theories such as technicolor [ 144, 145 1. 
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The inequality m„+ > m^o is one aspect of the general property mentioned in Sec. || [see in particular Eq. 
that vectorial interactions make positive contributions to the masses of physical particles. We can motivate this result 
(see Sec. ||) by considering the EM self-interactions and mutual interactions between arbitrary charge distributions 
(representing e.g. the extended constituent quarks and including effects of qq pairs as well). This leads to the 
conjecture that the EM contributions to hadronic masses are always positive. 

In order to extract genuine EM contributions to hadron masses we need to form A/ = 2 mass combinations such as 

(52[S] = + m^- - 2m^o 

S2[p] ^ 2{mp+ - THpo) , (12.19) 

in which the effects of {m'^^ — rn'p) ^ have been cancelled to first order. (Next order QCD effects are smaller 
than the observed splitting.) In terms of a naive quark model the E isotriplets are xuu, xud, and xdd states with 
X = s,c,b . . some heavy quark. The EM contribution to (52(S) comes from mutual qi,qj,i =^ j interactions and we 
find 

62^a\Qu-Qd\'' I Pn{r)pn{?')/\r-r'\>Q, (12.20) 



with pn referring to the density of the light u ov d quark, which in the /-spin symmetric limit are equal. Our 
conjecture is then that this positivity is not an artifact of the simple model but would persist in the full-fledged 
theory. Experimentally |]5^ , we find 

(52[S(1190)] = 1.5 ±0.18 MeV (12.21a) 
^2[S(1380)] = 2.6±2.1MeV (12.21b) 
(52[i;c(2455)] = 2.0± 1.6MeV (12.21c) 

<52[p] = "0.2± l.SMeV. (12.21d) 

The positivity of (52[S(1190)] and (52[E(1380)] are statistically very significant (in view of the fairly small width). 
Estimating ■mp+ — nipo is difficult because of the large widths Tp ~ 150 MeV.pl We also note that the positivity of 
the A/ = 2 energy shift should apply not only to the ground state, but also to any excited states. 

13. THE ABSENCE OF SPONTANEOUS PARITY VIOLATION IN QCD 

A beautiful result also proved by Vafa and Witten |^ via the QCD inequalities technique is that there is no 
spontaneous breaking of parity symmetry in vectorial theories. Such a breaking should manifest via a nonzero 
vacuum expectation value of some parity odd operator A. The simplest candidates for such operators are the local 
quantities constructed from the gauge fields only 

A^FF^ e'^-^-F^.F^,, e^'-^^i^.F^.Z^^F^^, . . . (13.1) 

Alternatively we could use quantities involving fermions such as ^ = 'tp'fstjj, with 75 — e^"^'^ "f^i"ivl\'^a- The dis- 
tinguishing feature of any parity odd operator is that it includes an odd number of g^'^-*"^ tensors. While the 
matrices, fields A^^ , and metric gp,i, are hermitian and real in the euclidean continuation, e^"^"' , which transforms like 
da;^da;^dx"^da;*, picks up a factor of i and becomes purely imaginary. Consequently Acuc corresponding to all parity 
odd operators is imaginary. 

If (A) ^ 0, then from ordinary first order perturbation theory it follows that a (hypothetical) theory with a modified 
Lagrangian 



^ Also for the special case of the vectorial p^p p multiplet there is a relatively important one photon annihilation contribution, 
which the argument presented above misses. 
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Cx = Cqcd{X = 0) + \ Jd^xA 



has, for small A, a shifted vacuum energy density 

E{X) = E{0) + \{A) , 



(13.2) 



(13.3) 



with (A) = {A)q still computed for the A = vacuum. For a A with an appropriate sign, X{A) < and £'(A) < E{0) 
- just as the energy of a spontaneously magnetized ferromagnet is lowered by adding an external B field anti-parallel 
to the magnetization jl. 

The euclidean path integral representation of the vacuum energy density for Cx is 



^(A) = -^ln / d[A^] I d[^] / d[^]e 



-J d^xC, 



(13.4) 



We will next show that E{X) > E{0), nega ting the possibility that E{X) < E{0) and forbidding (A) ^ 0. 
If A is of the type indicated in Eq. ( 13.1 ), the / d[i/'] / d\^\ integration can be carried out and 



E{\) = -lln J d[A,]e^'^ " "^^^/ ^[^^^ = ^^^^ 



(13.5) 



since the oscillatory factor e*'*'^ always decreases the integral. The same is true if (A) = {'ip'fsip) ^ 0. Adding XA 
to the Lagrangian is equivalent to introducing complex masses in the determinantal factor J^jDet(-p + nii) {e.g. 
rrii —f rrii -MA) which again destroys the positivity of the determinantal factor, reduces -^(A), and increases E{X), so 
that Eq. ( pA^ ) holds. 

The fact that Z{X) is minimal for A = allows also the proof that of all the different QCD vacua (which can be 
generated by adding the topological density OFF term to the Lagrangi an) , the one with the lowest vacuum energy 
is that with — 0. (Also 6 = tt was excluded by similar arguments |146|, a result suggested by dilute instanton 
calculations |107| and other considerations |H3l). 



14. QCD INEQUALITIES AND THE LARGE Nc LIMIT 



We have commented on several occasions on the inequalities in the Nc oo limit. We will next show Ig^l that in 
this limit we have indeed the stronger version of the baryon-meson inequality 

Nc 

TUB > —rriT, . (14.1) 

We will also indicate that the interflavor mesonic mass inequalities could be extended to non-pseudoscalars. 

The Schwartz inequality implies that the baryonic correlation function {B{x)B{y)) — J d/i(A) [rSA{x,y)'^r] with 

S'a(x,?/) the common fermion propagator, is smaller than J d^{A)tr"S\{x,y)SA{x,y) . This last expression 

represents the correlation function of products of n pseudoscalar currents, (0 | [JP*(x)]" [>/^'*(2/)]"| 0). We can show that 
in the large Nc limit this joint correlation function effectively factors into a product of [( J^"'(x) JP'*(y))]". Specifically, 
gluon exchanges between different qq bubbles in Fig. ^ do not modify the energy A^c™7r by more tha n (1), of the 
system viewed as separately propagating Nc qq pairs. Indeed, the color trace counting argument [37 147 1 indicates 
that the interaction energy between any {qq){qq) pair of bubbles is 0{1/Nc). Since we have N^2 pairs, which 
can interact via planar diagrams of gluon exchanges, the total effect is 0(1). The interaction of triplets (or a 
higher number) of bubbles is 0{1/Nc), or 0{1/Nc), k > 2, and is less important. Thus in the large |a; — y| limit 
Jd^i{A)tT-S\{x,y)SA{x,y) 



and since {B{x)B{y)) ~ e-'"^!^-?'!, Eq. ( [14T| ) follows 



To extend the interflavor mass relations consider the following three correlation functions of four flavor currents 
(see Fig. IH): 

^^1 = {Q\J^+{x^)JK+{x2)J.-{Vl)JK-{y2m 

F2 - (0|Jk-(xi)Jx+(x2)Jk-(2/i)Jk+(2/2)|0) 
^3 = (0 I J^-(xi)J^+(a;2)J',7-(2/i)^7r+ (2/2)10) , 

with the euclidean points xi,X2,yi,y2 forming a rectangle of height T and width r in the, for example, 3 — 4 plane. 
In the large Nc limit, quark annihilation is suppressed and we have the planar quark diagrams indicating the possible 
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contractions. In the case of F2 we specifically choose the contraction with ss exchanged in the vertical (T) direction, 
so that in all three cases we have uu exchanged in the horizontal (r) direction. We can use the euclidean Hamiltonian 
to evolve the system in T, e.g. 



Pi ^^{0\J7v+{xi)JK+ix2)\n) {n\J^-{xi)JK-ix2)\0) e' 



^|(0|J,-(xi)Jk+(x2)| 



\|2 -Er,T 

n)\ e " 



(14.2) 



with (0 I Jjr- {xi)Jx+ (2^2)1 n) independent, by translational invariance, of the specific location in T of xi, X2. 
We can also evolve the system in r to find analogous expressions: 



(14.3a) 



F2 = Y.\{0\JKAx2)JK-{y2)\n')f e-^'^'^ 



(14.3b) 



E 



\{0\J^+ix2)J.-iy2)\n)f 



(14.3c) 



where, due to the choice of F2, we have the same states 
inequality we have 



(with uu flavor) in all three cases. From the Schwartz 



|Fi|2<F2F3 



(14.4) 



which, by going back to the representation like Eq. ( 14.2 ) for Fi,F2, and F3, yields inequalities between the masses 
of the lowest intermediate states for the three cases 



(0) ^ 



1 



(0) , (0) 
^uu ' ^ss 



(14.5) 



since by construction the intermediate "T-channel" states for F2 are ss type states (and uu for _F3). As emphasized 
earlier in the large Nc limit th e m„a , mss sector s are distinct from m^^\ the flavor vacuum, and from each other. 
To get perfect squares in Eqs. ( 14.3b|) and (14.3c), we need to have identical currents at xi,X2 and yi,y2. Thus the 
quantum numbers of these uu, ss (and us) states are those of JttJtt, C = +, even G-parity, e.g. 0"*""*" states. 

Indeed, similar unitarity motivated inequalities were suggested earlier in the context of dual models and were applied 
to intercepts of Regge trajectories on which 0+"'" states lie [148|, and it has been conjectured that in the large Nc 



limit QCD a pproa ches a dual, string-like model. The present derivation pinpoints, however, the crucial element for 
deriving Eq. ( 14. 5| ). It is that in the large Nc limit we can meaningfully separate the various planar contributions to 
the correlation function and each evolves separately under the appropriate Hi^ (or Hij) Hamiltonian. 

We note that the last argument did not utilize the measure positivity in the euclidean path integral, but rather the 
unitarity based, more general, spectral positivity. QCD inequalities have also been applied in the large Nc limit in 
(1+1) dimensions [149|. 



15. QCD INEQUALITIES FOR GLUEBALLS 

QCD also predicts "glueball" states co nsist ing of gluonic degrees of freedom and no quark flavors. Experimental 



evidence for such states exists at present [150| - but mixing with the many qq and/or qqqq states in the 1 - 1.5 GeV 



region (where the glueball spectrum is estimated to start [151|) complicates the analysis. 

The local density F^ — F^^F^^^ becomes, in the euclidean domain, E'^ + and maximizes all other F bilinears 

such as FF ^ iE ■ B. Thus, as pointed out by Muzinich and Nair [PJlS^, we have F'^{x)F^{y) > FF{x)FF{y). 
Due to measure positivity this translates into an inequality between the corresponding correlation functions 



\F'ix)F'iym= J dfi{A)F'ix)F'iy) 

> J dfi{A) FF{x)FF{y) ^ {0\FF{x)FF{y)\0) . (15.1) 
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F'^{x),FF{x) create from the vacuum 0"'""'" and + states respectively. Thus, following the standard argument, we 
derive the inequality 

< m%-'' (15.2) 
between the masses of the lowest glueball states in the r espec tive channels. Confinement and the existence of a mass 



gap in the gluonic sector prevents interpretation of Eq. ( 15. 2| ) as a statement about gluon thresholds. 

Any polynomial in the comp onent of F^^, used to create a glueball state of arbitrary J^*^, is bound by an appropriate 



power of F^. Thus Eq. ( 15.2 ) can be generalized to 

'Tlib < (mass of any glueball) . (15.3) 



All these correlation function inequalities are automatically satisfied if F^ has a nonvanishing vacuum expectation 
value: (0|i^^|0) 7^ 0. In this case, correlation functions involving F"^ [or (F^Vjjhave constant terms as \x — y\ 00, 
correspon ding to the obvious zero mass 0"'""'" "state" : the y acuur n, an d Eq. ( |15.3| ) need not hold. 

In Ref. |l53|] . West suggested that the inequahties ( 15. 2| ) and ( 15.5 ) can be recovered if we use the time derivative 



of the correlator 

^(0|F2(0)F2(i,0)|0) = {0\F\0)[H,F^{t,6)]m 

= (0|F2(0)i7F2(t,0)|0) , (15.4) 

and the corresponding expression for FF. The offending constant term in the intermediate summation would be 
absent in the new inequalities. The idea then is to prove that 

{0\FF{0)HFF{t,d)\0) > {0\FF{0)HFF{t,6)\0) . (15.5) 

This could be accomplished by using the positivity of the Hamiltonian for each A^(x) background configuration 
separately in a path integral representation of the last two correlators. There are however subtleties in this argument 
which West tried to address. 

In general, while the Hamiltonian is positive when operating on physical states, the latter states are obtained only 
after the path integral has been carried out, and the pointwise positivity of H for every Afj,{x) configuration is not 
a priori guaranteed. It is the case for the pure Yang-Mills theory if H = J dx{E^ + i?^);|^ however, in this case 
H\Q) = 0, which has been assumed in the above argument, holds only after an appropriate subtraction has been made. 



Such a common subtraction will then modify Eq. (15.5) to 

{0\FF{0)HFF(t,6)\0) - c{0\FF{0)FF{t,6)\0) > 

{0\FF{0)HFF{t,6)\0) ~ c{0\FF{0)FF{t,d)\0) , (15.6) 



with c positive. Due to Eq. (15.1) this may invalidate the desired inequality. 

In passing we remark that (0_F^|0) 7^ has been conjectured to be the driving mechanism for confinement. Its 



value controls the bag constant |154], the string tension in extended hadronic states |155], and enters QCD sum rules. 

The glueball sector can be divided into "even" and "odd" parts consisting of states with quantum numbers of two or 
three gluons with appropriate orbital and spin angular momentum. Thus the J^'~^ states 0+"'',2++,0 ^,2 ^ belong 
to the even part and 1 to the odd. If we adopt a simple "constituent" description of the lowest even (odd) states 
in terms of two (three) gluons bound by one gluon exchange potentials, then the arguments of Sec. ^ leading to the 
baryon-meson mass inequality ms > (3/2)mM can be repeated here to prove the inequality 

> Im^^J (15.7) 

between the masses of the lowest odd and even states. Just as in the meson (baryon) 3 (X) 3 ^ 1(3 (g) 3 ^ 3), we have 
in the 2g (3g) state 818)8^ 1(8®8^8) and we can repeat the argument by replacing Ai . A2 in the fundamental 
representation by Ai . A2 in the adjoint representation of SU(3)(7. 



Actually it is the euclidean Lagrangian which has this positive definite E'^ + B'^ form, and the euclidean Hamiltonian then 
has the form J{B'^ - E^). A non- vanishing (0|ii/|0) VEV then amounts to the nontrivial dynamical assumption of the dual 
Meissner efi'ect in the QCD vacuum - namely the preponderance of large magnetic fluctuations. 
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Evidently, as the previous discussion suggests, a constituent valence model for glueballs is on much weaker footing 
than that for mesons and baryons. In particular, the latter (but not the f ormer!) can be justified in the large Nc 
limit. If the even glueball spectrum starts at 1.5 - 2 GeV, then Eq. (15.7) suggests that the lightest 1 glueball 
may be near the '4'' states. Its putative mi xing with the latter could then modify the perturbative approach to 
^ 3g decays and/or other matrix elements [156|. 



16. QCD INEQUALITIES IN THE EXOTIC SECTOR 



So far we have focused on correlation functions of QaQa bilinear and e'^'"^qaqbQc trilinear currents. In the following 
we utilize {J°'^{x)J^^{y)) with quartic "currents": 



q,{x)qj{x)qk{x)qi{x) , 



(16.1) 



to probe the "exotic" M^^j^i sector. We also address the possibility that the lowest lying intermediate state in {J{x)J{y)) 
is not a stable one particle state, but rather a two particle threshold at mi + 7712. 

This last issue has been encountered already for non-exotic currents. Thus, the lowest lying state in the vector 
1~ ud channel is not the p(760) but the tttt threshold at ~ 280 MeV. Since this happens because mp > 2™^, the 
Weingarten inequality nip > is not invalidated. 



Let us next analyze the general case 
both channels are resonances of widths ra(rf,) 
inequality mi°'' 



Ja(x)jt(y)|0) > (0|Jf,(a:)J^^(y)|0), when the lowest. 



If the states are narrow, i.e. 



Fa < 



non-exotic mesons in 
then a mass 



< mf''' can be derived. The point is that the correlation function inequalities are derived for all 
\x — y\. By varying — y| we effectively scan the spectrum since we keep changing the relative weight of different (fj?) 
regions in Eqs. (3.2) and (8.3). In the |a; — ^0 limit the high /i^ region dominates. The perturbative expressions 
for the correlation functions are adequate and indeed conform to the inequalities. In the other extreme, |a:: — ?/| — *■ 00, 
the threshold region dominates and we have, e.g., {VV) ~ exp{— 2m7r|a; — y\}. 

However, following lattice calculations of hadronic masses, let us consider "intermediate times" during which the 
local u{x)^i^^d{x) state, for example, evolves into a true p, qq bound state, but has not yet evolved into the final decayed 
form of two pions (see Fig. 13). In the range F^ ^ ^ \x — y\ ^ ™p from which meaningful mass values (and mass 
inequalities) could be extracted for the p resonance, (0| J^(a;) J^(0)|0) ~ exp{— mp|x — y\} is a valid approximation. 

In this section we focus on the threshold regions and the |a; — j/j —^ 00 domain in the correlator inequalities. Since 
near threshold the kinetic energies are much smaller than the masses of the stable particles, nonrelativistic kinematics, 
and some aspects of a potential model (in the tttt, ttK, and KK channels) may be applicable. 



F, < mi") 



The inequalities to|y ^ > 



2m j) [or mx 



,(0) 



often become in the continuum limit trivial equalities mi 



■ mi 



5 (2m, 



i(2m/f + 2m„) when we have confined i — u, j = s quarks and a meson- meson continuum]. 



This is analogous to the trivialization of ms < 2m7r into m27r < 2m7r mentioned above. 

If we consider systems confined to a volume of diameter R, we expect the various inequalities to be satisfied with 
a margin of order AE ~ l/2mR^, the level splitting for a system of size R. For the continuum, R — > 00, and AE 
vanishes. However, the relevant qua ntities become the phase shifts. We wish to interpret the inequalities as nontrivial 
statements about the latter 



To this end, let use consider the "tr^^ft,i2" inequalities (2.14) 



N N 
n=0 n=Q 



E. 



E 



(n) 



(16.2) 



with the sum extending over the first N excited states. If we put our system in a large box, then each time a "bound 
state" is generated (by increasing the strength of the attractive potential), the phase shift S in the relevant channel 

changes by tt. Levinson's theorem [158| suggests replacing the discrete sum J2n=o i / ^^'^ (1^) ^'^^ 
becomes 



i J^dSE,,{6) >liJ^dS [E,,i6)+E,,iS)] . 



(16.3) 



Evidently (16.3) reverts back to (16.2) in the narrow resonance approximation with ^ localized near the resonances 
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Near threshold S = ka, with k the center of mass momentum and a the scattering length. Using NR kinematics for 
kij {E) we can obtain the relation 

1/1 i\ /ii ii\ 

(16.4) 

Quarks are confined and for i — s,j — u in 0++, 1 , 2++, etc. channels, the lowest states are tt+tt", Ktt~, and KK, 
and one might use this inequality with i,j — > K,Tr. Scattering data for tttt and ttK ca n be ana lyzed by considering 
ttN — !■ 27r7V and KN KttN scattering and extrapolating to the one pio n exc hange pole 1 159, 160 1. This extrapolation 
is much more difficult fo the K exchange case, making a direct test of ( |16.4D difficult. Howe ver, t he existence of the 
0+"'" state in the KK system slightly below threshold is expected to enhance so that ( 16. 4| ) would most likely 
hold.g 

Consider next the exotic current with 0^^ quantum numbers: 

= Jij{^)Jli{^) = 4'i{x)-^z^.j{x)ipk{x)^zi)i{x) . (16.5) 

If all flavors are distinct, we have only one possible contraction, illustrated in Fig. pj](a), contributing to JsJs which, 
using Eq. (S.20), can be written as 

{0\Js{x)jl{ym = J d^iiA)tT[[S\{x,y)]^ S'^ix,y)}tT[[S'X{x,y)]^ S'Ax,y)} ■ (16.6) 

If m^*^^ = m^^\m^^'^ = rnf \ then the integrand in Eq. ( 16. 6D becomes the product of two perfect squares and (JsJs) 
maximizes all other exotic correlation functions. Espriu et al. |2^ speculate that this is associated with the fact that 
the lowest exotic qq qq states found in bag model calculations [ [161[ are 0++ states. 

Even more dramatic results follow ||2^ if we take all four, distinct, flavors, i,j,k, and I to be degenerate. In this 

case we have {JsJs) — /d/i(v4) |tr \^SA{x,yy SA{x,y)] } , which by the Schwartz inequality is larger than the square 
of the pseudoscalar propagator: 

{Js{x)Js{y)) > \{JP'Jn\' ■ (16.7) 

While the last inequality is consistent with having a bound S state in the 7r+7r+ (tt = ^75^, vr' = u'j^d') channel, 
it does not require such a state. The "5" state could simply be a tttt' threshold state and the inequality would be 
trivially satisfied as 

r7i(7r7r) threshold < + 771^^1 . (16.8) 

We would like to interpret this, in analogy with our above discussion, as a statement that the low energy (threshold) 
interacts attractively. This in turn enhances the density of threshold states relative to the non-interacting, free case, 
and enhances the long distance euclidean correlators. 



^°The inequality 

with J^+^/_ = Jk+Jk'-' '^k+ ~ ^a{x)'yzipu{x), etc. is readily derived, if we introduce additional degenerate flavors 
u' ,s' ,m^i = mu,mgi = rus. As \x — y\ —> 00 the various correlation functions are dominated by the respective meson-meson 
thresholds. If we neglect the interaction between the propagating mesons, the above inequality becomes a trivial equality since 
each of the two-point functions factorizes, e.g. (0| (x) J^_/ (x) (y) Jj^._' (i/)|0) = fnfKDn{x,y)DK{x,y) with = (0|J,r|7i") 
and Dtt the pion propagator. In this long distance, low energy limit, we can effectively treat the pions as elementary with 
interactions X-KnTr"^ , XkkK* , and XnKTv^ . To first order in the As these interactions modify the two-point correlation func- 
tion as follows: {Ji{x)J2{x)jl{y)jl{y)) = Dmi{x,y)D,Ti2{x,y) - A12 / d^x D,n^{x, z)Dm2{x, z)Dnii{z,y)D„i2{z,y). The 
second term represents the effect of one mim2 collision. By going to momentum space the z integration can be done and an 
inequality between the As is obtained [^. Unfortunately, it contains, besides the masses, also a subtraction point fio and will 
not be reproduced here. Note that the minus sign in front of the first order euclidean perturbative contribution reflects the 
expansion of e^^'. Thus a negative A12 - corresponding to an attractive nn interaction - enhances the joint propagation. 
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The connection between an attractive potential in the ij channel and the propagator {0\j}j{x)Jij{y)\0) can be 
directly seen in a potential model limit. With one of the particles infinitely heavy, the ij propagator is essentially 
that of the other particle moving in the attractive potential. The path integral expression for such a propagation is 



paths 



/-^^^exp/-/ 



paths 



"d7 



v{xn 



(It 



(16.9) 



with r a "proper time" and a;'^(0) = x,x^{t) = y. Evidently an attractive potential {V < 0) will enhance the positive 
contributions of the individual paths. This is indeed expected from the interpretation of the free euclidean propagator 
as a diffusion kernel: the probability of the diffusing particle to return to the origin is clearly enhanced by an attractive 
potential. A more rigorous and systematic approach directly relating euclidean lattice corr elators to phase shifts and 
scattering lengths was suggested by Liischer |162 . It was utilized by Gupta et 
scattering length is positive. 



ISTfl to prove that the tt — tt 



The idea in Liischer's approach (see also Neuberger's 163 suggestion of a calculational lattice method for deter- 
mining /tt) is to use the finite size corrections to the correlators. The latter are the conf iguration space analog of the 
1/R energy shifts discussed here. We refer the reader to the original Liischer work |l62| for further details. 

We have argued that m(7r7r') < (ttItt + m^r' ) can be fulfilled by an interacting 7r+7r+ threshold state if the scattering 
length is attractive (positive) . The tt and tt' can only interact via gluon exchanges since they ar e com posed of different 
quarks. The interesting fact that this interaction is attractive is in accord with an old result |164] that the retarded 
Van der Waals forces between systems of identical polarizabilities are always attractive. 

Indeed the Casimir-Polder two-photon exchange interaction: 



Vcp{r^ 
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remains attractive so long as the ratios of electric and magnetic polarizabilities aEjotu for the two neutral systems 
(1) and (2) in question are similar. Our above result generalizes this attractive nature of the photon exchange to the 
full non-perturbative case, and applies also for nonabelian gauge interactions. (If, like in QCD, the theory is confining, 
then there is a mass gap in the pure glue sector, and we expect an exponential rather than power law falloff of the 
interaction. The attractive nature does however persist.) 

We note that two independent lines of argument can suggest that the two-photon exchange interaction is attractive. 
The first nonrelativistic, second order perturbation theory argument applies if the systems (1) and (2) considered 
are in their respective ground states. The second, more general, relativistic field argument uses t-channel dispersion 
ifzol and the positivity of the corresponding spectral functions when (1) = (2). Amusingly the conditions required 
for proving the inequality (16.7) simultaneously conform to both arguments, since the «j, kl pseudoscalars are indeed 
the lowest states in their channels, and taking — rrik , rrij — mi makes them (dynamically) identical. Finally we 
would like to mention that the second argument suggests |165| that contrary to some lore, two halves of a conducting 
spherical shell attract rather than repel. 

In the real world the tt^tt^ scattering length is repulsive |166|. Indeed, with u — u'^d — d! we have the alternate 
contraction of Fig. |ll|(b). It has one fermion loop and makes hence a contribution of the opposite sign: 



{JsJs) = / d^Ji{A) \ \tv{S\SA)\ -tr \{S\Sa)Hs\Sa) 



(16.10) 



so that Eq. ( 16.7 ) can no longer be proven. Indeed at short distance, when t he me sons' wave functions overlap, 
we expect a repulsive Pauli effect which is reflected in the minu s sign in Eq. (16.10). We have not been able to 
show that this extra negative term reverses the sign in Eq. (16.7), so that no tt+tt"*" exotic bound state exists and 



r(/ = 2) is repulsive. 



17. QCD INEQUALITIES FOR FINITE TEMPERATURE AND FINITE CHEMICAL POTENTIAL 

Recently there has been much interest in QCD at finite temperature and finite baryon density, i.e. finite chemical 
potential. This interest is partially motivated by the desire to better understand compact neutron / (strange) quark 
matter stars, and by the prospect that heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) can indicate 
the expected finite temperature phase transition. 

In the following we would like to comment on the possible relevance of QCD inequalities in these cases. First, we 
note that all the correlator inequalities are maintained for T > 0. In the euclidean formulation, introduction of finite 
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temperature is simply equivalent to the imposing periodicity /3 = 1/T in the time direction |167]. The restriction 



of the gauge field configurations in the euclidean path integral to such periodic configurations clearly does not spoil 
the positivity of the measure. Also the relation 755'^(a;, 0)75 — S']/'(x, 0) and the ensuing positivity of the fermionic 
determinant and the integrand in the pseudoscalar correlators are maintained at finite temperature - and therefore 
so are all the correlator inequalities. 

Precisely because of the periodicity in time, we cannot use the Hamiltonian and its i — > 00 limiting e^™* behavior 
in order to infer bounds on masses smaller than T = 1/(3. Therefore we cannot attempt, when T 7^ 0, to prove that 
the axial global symmetry spontaneously breaks down, as this feature is closely tied to the massless pseudoscalar 



Nambu-Goldstone bosons. Indeed, numerous thoeretical arguments [16J,169] and lattice simulations have virtually 
proven that in QCD there is in fact a phase transition corresponding to axial symmetry restoration (and quark 
deconfinement) at a temperature ~ Aqcd- The exact character of this phase transition and its dependence on 
Nf, the number of quark flavors, was for a long time unclear However, we believe that the Vafa-Witten results 
regarding the non-breaking of vectorial global symmetry [which rely on Sp{x) < Sp{x)] and even more so the one 
concerning the nonbreaking of parity (which depends only on bulk, i.e. free energy, properties) continue to hold for 
finite T QCDJ^ 

The QCD inequalities techn ique may be even less useful in discussing transient varying T phenomena, such as 



disoriented chiral condensates [|172| , suggested to occur in domains of cooling quark gluon plasma. 

The introduction of a finite chemical potential, i.e. consideration of QCD in a background of uniform baryon 
density, has a much more drastic effect. It amounts to changing the fermionic part of the euclidean Lagrangian to 

Nf _ 

-^(m) = X! [V''7/i-D^V'i + ^V'OoV'i + rriiipiipi] , (17.1) 

so that the new /i 7^ propagator no longer satisfies 755*^75 = Sp, but rather 

75S'i=^75 = 7577- — 75 = — ; 7^ Sl , (17.2) 

and hence the positivity of the fermionic determinant can no longer be inferred. This not only prevents the proof of the 
QCD inequalities, but also excludes the utilization of lattice numerical simulations in which the statistics of occurance 
of lattice gauge configurations prescribes their (positive!) weights. The vario us d ramatic speculations concerning the 
high /i phase (involving parity breaking and even "color superconductivity" [ |l73t ) are therefore not excluded. 

The special case of = 2, i.e. SU(2) gauge theory, is a notable exception, and one can show mea sure positivity 



in this case even for ^jl ^ 0. This result, which has been utilized for some time in lattice simulations |174| is readily 
proven by using 

I5CI2SFI5CI2 = D* . (17.3) 



In Eq. ( 17.3 ) C is 17072 and I2 is the generator of the SU(2) color isospin, and it holds for arbitrary fi. This feature, 
which is due to the pseudoreality of SU(2), is essentially the same one used by Anishetty and Wyl er |80| and Hsu 
to extend the inequalities to chiral SU(2). It has been used by Kogut, Stephanov, and Toublan |l74[ in order to 
prove that the correlator of the 0+, / = {i.e. antisymmetric in flavor) tpip combination 

can serve, for fj. ^ 0, as an upper bound for any other correlator. This implies that for n =/= the 0^ ipip diquark is 
the lightest boson. This complements the claim that for /i = and Nc = 2, the 0+ diquark is degenerate with 
the 0~ ip j5tp pion (see the end of Sec. ^), and nicely fits with the unique patterns of symmetry breaking suspected in 



this case [174] 



and could be Nf dependent. We will later make a conjecture that Tc 



^^It is generally believed to be weakly first order |170 
monotonicall y dec reases with Nf. 

^^T. Cohen 171| utilized the QCD inequalities to suggest that not only is SU(3) axial symmetry restored above Tc, but also 
\]{1)a. Since the latter is broken via the QCD anomaly and not spontaneously, this result is rather surprising. Indeed as noted 
by Cohen the proof involves a technical assumption which may fail. 
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External electric and magnetic fields modify the hadronic spectrum. Indeed lattice calculations utilized external 
magnetic fields to prove, via the jl-B interactions, the hadronic magnetic moments. Also very strong fields (the analogs 
of supercritical fields in superconductors) can even modify the phase structure of spontaneously broken gauge theories. 
In the context of QCD and other vectorial theories the measure postivity is clearly maintained in the presence of the 
external B field. However, QCD inequalities such as BE (parapositronium) > BE (orthopositronium) or rrip > TOtt 
may be modified since the vector (triplet) state has a magnetic moment, and for a sufficiently strong external B 
field {\B\ > Aqqp), the state with jl antiparallel to B may become lower than the singlet pion. In this connection 
we recall the amusing suggestion [ 175 | that in strong enough magnetic fields, existing in appropriate astrophysical 
environments, that this effect can reverse m„ > mp to mp > m„ with an inverted (3 decay! 

Unlike the effect of high temperature and/or chemical potentials, we believe that strong magnetic fields do not 
cause deconfineme nt o r chiral symmetry restoration. Indeed in the limit B ^ h? we expect, in analogy with the case 
of atomic physics |17(:], that in strong magnetic field stars the motion of the quarks becomes one dimensional, along 
the B field lines. It is well known that in such cases confinement only gets stronger (even U(l) theories confine!), 
and so should SxSB. It is an interesting, open conjecture which we would like to make here, that even moderate B 
fields tend to enhance S^SB, e.g. by enhancing the density near A = eigenvalues of the Dirac operator and thus, 
according to the Banks-Casher criterion, enhance the (V'V') condensate. 



18. QCD INEQUALITIES FOR QQ POTENTIALS, QUARK MASSES, AND WEAK TRANSITIONS 



The main application of the techniques developed above is to obtain inequalities between directly measurable 
quantities such as hadron masses. H owever, there are several calculational approaches to QCD such as the potential 
model for heavy quarkonia [177-180|, chiral perturbation theory for the low-lying mesonic sector 181 1, and QCD sum 
rules Each of these schemes depends on a few input parameters and makes many predictions. Applying the 

techniques of QCD inequalities to these input parameters could therefore yield a very large body of suggestive results. 

The area in which most research along these lines has been done (in a large measure prior to, and independent 
from, the introduction of QCD inequalities in 1983) is that of potential models. The point is that given some general 
properties of the potential such as convexity of V{r) or other feature s, we hav e a wealth of information conce rnin g 
the level ordering, thanks to the work of Martin and collaborators 182, 183 1; Lieb |l^; Fulton and Feldman [184|; 
and others. 

An early important observation was that the P-wave excitation in cc (or bb) systems is lower than the first radial 



excitation: E,, 



> En^.i+i- This deviation from the famous degeneracy in the pure Coulombic case is related 



to the fact that the QCD potential Vqcd('') does not satisfy V^Vqcd 
assumpt ion of a monotonically increasing, ^ > 0, convex potential, 
Martin ||l8|l (BGM) to prove En,i < E^-ij+l 

> (or if e""^^ has a positive Fourier transform), then the baryon mass relations Eqs. (5.11a) and 



= but rather V^Vqcd > ||l8j. Also the 
< allowed Baumgartner, Crosse and 



In addition if 



( ^.llb| ) can be proved (see ||T8|). These relations, which prescribe the sign of deviation from linearity of masses in the 
decuplet or from the Cell-Mann - Okubo relation in the octet, were proven above only under the explicit additional 
assymption of flavor symmetric wave function s, a point emphasized b y Rich ard and Taxil . [Indeed as shown by 
Lieb |l^ and by Martin, Richard, and Taxil |185], Eqs. ( ^.lla ) and ( 5.11b ) are violated for potentials V = r"^ with 
sufficiently large a.] 

It is not clear at the present what are all of the model independent statements about Vqq (i?) that can be proven. 
We would like to mention, however, the very elegant result of Bachas [186| on the convexity 



ViR)> 



V 



R 



V 



R 



(18.1) 



of V{R). Together with earlier work by Simon and Yaffe |187| which shows a monotonically increasing V{R), this 
promotes the BCM level ordering into a QCD theorem, which is confirmed in the cc system (and partially tested in 
the bb system). 

Let us next reproduce Bachas' proof. The static potential is definted in terms of a rectangular Wilson loop W |7^] 
in the {t, r) plane of height T — > oo and width R: 



V{R) = lim 

1 -^OO 



-iln(trC/(W^)) 



(18.2) 



with U {W) the ordered product of [/unks around the loop W. We have the path integral representation 
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{ivU(W)) = 



/niinksdM(^)trL/(W-) 

/ Hunks dM(C/) 



(18.3) 



with dfi(U) ~ Y['^[U] exp (^—^^ J2n tr^^p the positive lattice measure. We can write W — Wi{4>W2), where (fAV2, 
the line reversed version of the reflected path, denotes the reflection of the portion of W to the left of the hyperplane 
indicated by the dotted hne in Fig. |l|. Obviously U{x) = W{x) and thus tiU{W) = tr [U {Wi)W {W2)] . We can 
compare {trU{W)) with the corresponding expressions for the symmetric Wilson loops of size R + r, R — r obtained 
by joining Wi and (pWi or W2 and (j)W2. By dividing the d/i(f7) integration into variables to the left of, right of, and 
on the reflection plane, the expectation values of (tr(VFi0Wi)) and {ti {W2(j>W2)) can be written as perfect squares, 

e.g. {tiiWi(j)Wi)) ^ J d^iontr [(/ dfiRWi) (/ duRWi)^^ = J d^o„ tr [(/ dfiLWi) (/ d^iL W^i)^] and the Schwartz 

inequality implies {tiU iW2)U {(jiWi)) ■ {tTU{W2)U {<j)W2)) > \{tTU{W2)U{(j)W2))\^ . The definition ofV{R), Eq. ( pA2| ), 
then readily yields the desired conve xity ( |18.1 ) 



Similar arguments were applied [188| to the Eichtcn-Feinberg |139| representation of the tensor and spin-spin 
potentials. The inequalities obtained are basically in accord with a scalar long range confining potential [189|. 



The masses of the different quark flavors in the QCD Lagrangian (or Hamiltonian) H — Ho + J^fji "^f'^i^i ^'^^ the 
only explicit dimensional parameters. "Dimensional transmutation" generates, however, an additional scale Aqcd (« 
a few hundred MeV). Thus, unlike QED with rBoiu- — (we)^^ and Coulombic binding ~ me, we expect that scaling 



all quark masses m^°^ cmi' changes physical masses (lengths) by less than a factor c (1/c). As we next argue, the 
simple linear dependence of -ffqcD on rnf^ further restricts the variation of mass parameters as a function of Tm^f^ . 

In quantum mechanics the ground state energy E'^^^ (A) is a convex function of any set of parameters that the 
Hamiltonian depends upon linearly 190 |. If H{\) = Hq + ^iHi = Hq + X ■ H, and A = afl + (1 — a);?, then 
£;(°)(A) < aE^^^fl) + (1 - a)E^°^{iy). [This result follows immediately from H{X) = aH{fl) + (1 - a)i?(i7) by taking 
expectation values in ^/'('^^(A), the ground state of -ff(A), and using the variational principle.] 

We wish to apply this result to the masses m[^\j^'^) of ground state mesons with different flavor and Lorentz 

quantum numbers. We are hindered by the fact that we are not free to vary 771^-"'' and by the need to subtract the 
vacuum energy. However, in a large Nc approximation where QiQi do not annihilate and effects of closed quark loops 
are neglected, it is meaningful to describe each sector rriij by a separate Hamiltonian Hij — E[q^y) + ■ 

The knowledge of the masses of the , 1 , 2++, . . . flavor multiplets can then be used to constrain ratios of quark 
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masses like R = (m, 
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) 1 191 1 . The result obtained is somewhat large but still consistent with 



the estimates of Gasser and Leutwyler |181| , and correspond to a lower bound mi"'' > 80-100 MeV.p^ 

Over the last twenty years there has been an ongoing effort to address in a systematic way the low energy sector 
of hadronic physics, and in particular the sector containing light quarks only via chiral perturbation theory (xPT). 
The idea is to incorporat e SySB, the Goldstone pions and current algebra, and the ensuing low energy theorems 
via effective Lagrangians [194, 195 1 which will manifest the desired symmetries and which are written in terms of the 
pionic fields only. Then one uses these chiral Lagrangians to perform a systematic expansion in the external momenta 
of the problem, and/or the pion mass divided by some "hadronic scale" usually taken to be 47r/^. These effective 
Lagrangians contain a series of terms £±,£2, ■ ■ ■ which in general are ranked according to the number of derivatives 
appearing in each term. One then can also systematically compute higher loop processes. A particularly simple and 
elegant form of such an effective Lagrangian is the Skyrme model |196[ , which even incorporates the nucleon as a 
soliton state. 

In principle all the coefficients of the terms in the effective Lagrangian are computable from QCD. In practice they 
are often fixed by fitting some low energy data. In any event, QCD inequalities can often be applied to constrain 
the range of these param eters . While no systematic program of this kind has been completed, some steps have been 
taken by Comellas et al. [197[. By expressing the currents in the inequality relating vector and pseudoscalar currents 



^^It is amusing to note that the relatively "large" e'/e ratio (as compared with previous theoretical estimates [192|) recently 
found in high precision experiments on CP violating kaon decays [193] do indeed suggest a relatively low mi'^' ~ 80 MeV. The 
fact that we have so far been able to obtain only the above modest lower bound, rather than, say, mi°' > 190 MeV, makes it 
easier to accomodate the measured e'/ e in the standard model. In passing we note that the smaller mi''-' value would also favor 
the stability of strange quark matter |192,L9S[. Note however that if the new lattice calculation of e'/e (see next Footnote) is 
adopted, the above remark need not apply. 
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as a xPT expansion in the pionic field, and also by using the momentum space version of these inequalities, bounds 
on parameters in the effective Lagrangian were obtained, and are well satisfied. 

Over the last decade the heavy quark approximation was often used in connection with Qq or Qqq systems. This 
is based on the realization that in the infinite mg limit the heavy quark simply becomes a static color source Icadini 



to some universality relations, heavy quark symmetry, and a systematic expansion in inverse powers of toq [198 200|. 
The Witten pseudoscalar mass inequalities such as 2m(g, Q) > m{q, q) + m{Q, Q) become trivial in this limit where 
the QQ system is essen tially Coulombic, wi th in finite binding a^niQ/r. Despite an interesting effort on the part of 



Guralnik and Manohar 201 |, it is not clear [p02[ how this can be amended to yield useful inequalities. 

Precise information on nonperturbative QCD parameters of weak decay is of particular importance and may decide 
the fate of the Standard Model with the three generation KM scheme [^03[ . It is worth pointing out that some 



inequalities between such matrix elements can be motivated [204|. Let us then consider the K — )■ tttt decay. In 
addition to the standard left-handed four fermion operators [say A = s{x)j^Lu{x)u{x)"/^Ld{x), with L = 1 — 75 
the left projection operator], we can extract from "penguin diagrams" additional mixed terms of the form B — 
s{x)j^Lu{x)u{x)j^Rd{x), with R = 1 + 75. After using soft pion and current algebra techniques to reduce one pion 
we need to evaluate {K\A{or B)\tt) matrix elements. The latter can be related to the asymptotic limit when \x — y\ 
and \y — z\ —^ 00 of the three-point function {0\s{x)jzu{x)A(jj)u{z)j5d{z)\0) . The latter has the path integral form 

{K\A\tt) = / dfi{A)tT["f5SA{x,y)jf,{l-j5)SAiy,z)"fr-,SA{z,y)"f^,{l-j5)SAiy,x)] , 



corresponding to the unique contraction in Fig. ^ A similar expression with the second 1 — 75 — > 1 -f 75 applies 
when A B. We have assumed a flavor symmetric limit and used the same propagator Sa for all the quarks. Using 
^^SAix, y)^5 = S^{y, x), the hermiticity of the euclidean 7^, and (757/^)^ = —757/^, we can show that the mixed L — R 
expression (case B) has a path integral with an integrand which is an absolute square. It is therefore larger than the 
expression for the "pure" L ~ L case {A) and allows the conclusion that 

{K\B\tt) > {K\A\n) . 

While the soft pion limit and flavor symmetry assumed in the derivation are considerably weaker than this result we 
do still flnd it interesting and suggestive]^ 

19. QCD INEQUALITIES BEYOND THE TWO-POINT FUNCTIONS 

Most of the preceding sections, and of the work on QCD inequalities to date, has focused on euclidean two-point 
functions. These correlation functions arc sufficient for obtaining hadronic masses via the spectral representation. 



An obvious advantage of two-point functions is that the path integrals expressing them [Eqs. (8.10a) and (8.10b)] 
contain products of quark propagators SAix,y) between the same points x and y. Thus it is easy to prove positivity 
of certain combinations, e.g. tr[S]^{x,y)SA{x,y)] and the ensuing inequalities between the integrands of the path 
integrals for various two-point correlation functions. These algebraic inequalities which are true "pointwise" for each 
external Afj,{x) configuration do survive the path integration with the positive measure dii{A). 

Can we make general statements in the form of inequalities also for three-, four-, and higher point euclidean corre- 
lation functions? Consider a generic four-point correlation function {Ja{x)Jb{u)Jc{y)Jd{v))^ with a, 6, c, d referring to 
Lorentz and flavor quantum numbers. By appropriate contraction we obtain an expression for correlation functions of 
the following general type, where for simplicity we assumed degenerate quark flavors and hence the same propagator 
Sa: 

{Ja{x)Jb{u)Jc{y)Jd{v)) = d^{A)tr[raSA{x,u)TbSA{u,y)rcSAiy,v)TdSA{v,x)] . (19.1) 



All propagators Sa in the path integral refer to different pairs of points and it is not obvious how to construct positive 
definite combinations for each external Af^{x) conflguration separately. 



■^''it is worth noting that the B matrix element with the "8" construction is one of the ingredients fixing the positive sign of e'/^ • 
The QCD inequalities are generally not operative when we have disconnected flavor cont ractions as in the "eye" contraction of 



Fig. hq. Since the latter appears to dominate the lattice calculated matrix element [205|, we unfortunately cannot fix the sign 
of this very important quantity by QCD inequalities alone. 
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However, if we wish to study hadronic properties beyond the mass spectrum such as couphngs, scattering amphtudes, 
weak interaction matrix elements, or wave functions, we need to c onsid er more than two-point correlation functions 
[ p06| . In particular correlation functions of the form given in Eq. (19.1) have been used in order to study the charge 
distribution of mesonic states. Let the current 



J1 = ■4^t{x)Ta^j{x) , 



(19.2) 



with X — (— T, 0), create at the origin, at some remote past instant, a quark and antiquark of flavors «, J in a specific 
Lorentz state. As the system evolves under the QCD Hamiltonian it settles into the wave function of the ground 
state meson in this channel. Specifically, we can infer from the spectral expression Eq. (B.3) that after time T the 



components of excited states in the wave function of the system are suppressed relative to the ground state amplitude 
by e~"^^™, with Am the mass gap between the ground state and the first excited state. If we wish to find the relative 
separation of the quarks in the ground state, we can probe the system again with two external currents Jb{u) Jd{v) 
at u = (0,ri), V = (0,r2), with ri,2 = (i?±r)/2 and r the relative separation. The currents Jb,Jd should refer to the 
flavors i and j of the quarks respectively 



Jb = ^i^^bipt 



Jd = iijTdiij , 



(19.3) 



and we will take F;, = = F. Finally, in order to obtain the complete gauge invariant correlation function the two 
quarks are propagated back to the origin where they are annihilated at w = (T, 0) by the current Jc = J\- Various 
attempts have been made to meas ure, via lattice Monte-Carlo calculations, charge distributions (or form factors) for 
the ground state mesons [207 208 1. In the following we will focus on the particular case of the pion, i.e. using 



JI = V'«75'0d(a;) = Jc 

Jb = ^u^^u, Jd = ^dT^d 



we define 



F{t, f ) = d^R (O 



jt(-r,o)jJo, 



R 



Ja(T,0)Jd 0, 



R 



(19.4) 



This function has the following general properties [105|: 

1. F{T,r) < F{T,0), i.e. the "charge density" is maximal at the origin. 

2. The Fourier transform {i.e. the "form factor") 



L{T,p) = / e'P-^F{T,r} 



is positive for all T and p values 



L{T,p)>0. 



(19.5) 



(19.6) 



Evidently (2) implies (1): F{t,r) is also the Fourier transform of the positive L{t,p), and hence has its maximum 
value / L{t,p)d?p at the origin (r = 0). 

In order to prove (2) let us consider the four-point correlation function 



(19.7) 



F= / d^i{A)tY[^r,S\{x,u)TS\{u,y)^r^S\{y,u)TS\{u,x)\ 

d^i{A) tr { [S\{u, x)] ^ f [S\{y, u)] ^ S\{y, u)rS\{u, x)] , 



where we have used the "charge conjugation" property '^^SA{x,y)^5 — S\{y,x). f is defined as 

f = 75r75 - ±F = ±Ft . 



(19. 



Using u,v = (0,fi_2) we take the Fourier transform with respect to r = ri — r2, assigning momentum p{—p) to the 
quark (antiquark) at ri{r2). Exchanging the order of the d/i(A) and dridr2 integration the Fourier transform can be 
written as 
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L{T,p)= f dn{A) f dn f dr2e'P<^'--^^hr\^S{y,u)rS{u,x)[S{y,v)rS{v,x)]''^ 



dn{A) tr 



dfie'P-^'S{y,u)TS{u,x) 



(19.9) 



and the manifest positivity of the integrand persists after the dfj,{A) (> 0) integration, yielding the desired result 
L{T,p) > 0. 

While F{t,f) and L(t,p} are well-defined, gauge invariant, and, in principle, measurable quantities, the sugg estive 



interpretations as "charge density" and "form factor" are more heuristic. In particular it is not evident from Eq. (19.4) 
why F{T, r) is positive definite as T ^ oo. Here we interpret F as |-07r(?')P, with -0^ the "pion wave function" . Recall, 
however, that, as emphasized in the conclusion of Sec. ^ the positivity of norms [and spectral weight functions in 
Eq. (|8.3| )] emerges only after the d/i(A) integration has been performed and need not be manifest for each external 
Afj_{x) configuration separately. 



As f2 approaches fi, S{y, v) and S{y, u) [and likewise S{x^ u), S{x, v)] in Eq. (19.1) connect one vertex to two nearby 
points. For a typical smooth A^j_{x) configuration we expect S{y,v) S{y,u) and the mixed products will gradually 
become squares (see Fig. |7|). This suggests that F{T,r) is not only maximum at r = but is also monotonically 
increasing towards r — Q. This property does not hold for all gauge configurations, but only after averaging, since we 
expect smooth A^{x) configurations to dominate in J d/i(A) . Thus a proof of the conjectured monotonicity of -^(r, r) 
requires understanding the correlation between (x) along neighboring paths and would depend on the specific form 
of the Yang-Mills action Sy m ■ This action prefers small variations of the gauge fields over a small region [gradients 
of yield E and B and Sym = J d^x {E^ + B^)]. 

If we identify F(T = oo, r) with the pion's charge density, then it is a musi ng to note that the monotonicity of the 
ground state wave function can be proven in nonrelativistic quark models |209|[ wh en the potential is purely attractive. 
The fairly simple argument employs the spherical rearrangement technique [ ^10[ , showing that given any trial wave 
function will always lower the energy by "rearranging" its values into a radially monotonically decreasing set. 

The present discussion of four-point correlation functions applies only to pseudoscalar currents Jc,[x)Jc,{y). Indeed 
it is precisely for the spin singlet S-wave "pion" state that all components of the nonrelativistic quark model potential 
- the confining linear part, the Coulomb force at shorter range, and the very short range, hyperfine color magnetic 
interactions - are attractive. 



20. SUMMARY AND SUGGESTED FUTURE DEVELOPMENTS 



We have presented above many inequalities for hadronic masses and analyzed the possible theoretical and phe- 
nomenological aspects. Many results follow essentially from the positivity of the measure in the functional path 
integral for euclidean correlation functions in QCD (or other vectorial theories). This may be on occasion comple- 
mented by fairly mild assumptions on the number of light degenerate flavors or the Zweig rule (large N^) suppression 
of qq annihilation. These results include the Weingarten mass relations, the two Vafa-Witten theorems, the convexity 
of the QQ potential, 171^^+ > 771,^0, and Witten's interflavor relation for pseudoscalars. 

The fact that so many results, which have far-reaching implications, can be proven with such minimal input, is 
truly fascinating. We believe, however (and will try to make slightly more concrete conjectures later), that many 
more results would follow if we appeal to the specific form of the QCD action and in particular to its "ferromagnetic" 
character. 

A large class of meson-meson and baryon-baryon mass relations follow from the flavor independence (apart from 
the explicit mass terms) of Hqqd ■ This allows us to prove operator relations for Hqcb restricted to different flavor 
sectors, from the mass relations following - though only for flavor symmetric wave functions. Under fairly general 
assumptions we can also use Hamiltonian variational techniques to prove detailed baryon-meson inequalities. All 
of the above together with some general level ordering rules for the different J^*^ states (again motivated by QCD 
inequalities applied to qq potentials) can serve as extremely useful "Hund like" rules of thumb in the hadronic domain. 
This in turn can restrict the masses (or J^*^ quantum numbers) of new flavor combinations or radial excitations. 

In many cases the inequalities are more sophisticated versions of relations suggested by a naive quark model. 
Eventually we hope that much of the vast information available on hadronic parameters (including scattering, wave 
functions, etc.) will be constrained by such inequalities. 

On the more theoretical side the generic properties of QCD and all vector QCD-like theories should be analyzed. 
The constraints imposed by the inequalities on composite models of quarks and leptons, together with the anomaly 
matching conditions, are particularly interesting. The fermion-boson inequalities, m^"* > m^^\ exclude the protection 
of small masses for composite quarks and leptons via an unbroken global axial symmetry, in all cases when we 
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have vectorial, QCD-like, underlying dynamics. Such constraints are avoided by going to chiral gauge models or 
(supersymmetric) models with scalars where the measure positivity dfi{A) is lost. (In passing we note that large 
classes of scalar and scalar plus fermion theories are ruled out by "triviality" difficulties |211|.) 

We have often referred to the euclidean correlation function approach to deriving the QCD inequalities as "rigorous" 
and to that based on the Hamiltonian variational approach as more "heuristic" . It should be emphasized that this does 
not reflect any true objective distinction. Clearly the Hamiltonian and Lagrangian approaches to classical mechanics, 
and to quantum field theory, are equivalent and equally rigorous (or not) depending on the practitioner. 

The Hamiltonian variational approach allows inputting various approximations and/or physics intuition much more 
readily, however. These inputs can motivate certain restrictions on the wave function(^l)s allowed. An example is 
our restriction in the nonperturbative derivation of the baryon-meson inequalities (Sec. |6|) to baryonic configurations 
with only one junction point (see Fig. |4|). 

Clearly any restrictions on the class of wave function(al)s allowed potentially pushes the energy of the ground state 
higher. Hence, in particular, our restricted baryon is not the true ground state baryon. If, however, we can assume 
that any network containing some extra junction (and in this case, also extra anti-junction) points corresponds to 
"massive" components in the wave function(al) (with essentially extra virtual baryon-antibaryon pairs), then such 
components are likely to be small in the true ground state, and thus fairly safe to neglect. Similar comments apply to 
the neglect of possible sextet internal fluxons in our ansatz functionals for the quadri- and pentaquarks (see App. ^). 

Much of the intuition for hadronic physics stems from naive quark models. In these, the baryon-meson wave 
functionals are approximated by gcon9con9con and gcon^con wave functions. Here gcon refers to "constituent quarks", 
namely some effective quasiparticles obtained when some short distance modes are integrated. Clearly this concept 
would be on firmer ground if we could indeed show that the mass of the constituent u and d quarks - which in the 

mi"' — m^"' — axial SU(2) symmetry limit is purely dynamical - is indeed generated by physics (instaneous, or 
other) operating on distance scales smaller than the radius of the baryon or meson.^ 

Because the constituent quarks are extended, complex entitities, the "potentials" acting betwen them are in general 
rather complex - with spin-, flavor-, and possibly energy-dependence and non-locality. Still, the naive expectations for 
potentials generated via gluon exchange seem to lead to very successful predictions. Overall the naive quark model, 
which predated QCD by almost a decade, still provides more insight and results than any of the more sophisticated 
effective Lagrangian approaches. Justifying its usage from first principles would therefore constitute a major triumph 
and is certainly a worthwhile effort. In our discussion of the QCD inequalities we have encountered on several 
occasions [e.g. in discussing the pion mass and wave function, the interflavor relations, and the positivity of EM 
and more general vectorial interaction energies (see Sees. ^ and ^ and App. |^)] a remarkable coincidence between 
the "naive" quark model and rigorous QCD inequalities. Thus the efforts to extract constraints on the interquark 
potentials from QCD, and then to utilize this rather broad class of potentials allowed to derive level ordering theorems 
and/or baryon-baryon mass relations (see App. ^), is quite worthwhile. 

Even when we have mass dependent 9con9con interactions, the system does retain its symmetry under q™'^ <-> Qj""^, 
and one may thus wonder if just this single feature, abstracted from the constituent quark model, is sufficient to 
supply the flavor symmetry assumption which was the missing link in our proof of meson the interflavor relations 
in Sec. |^. Unfortunately this is not the case. Essentially the operation ipl^\d) ip^^\u) in Sec. || corresponds 
to exchanging the "fundamental", pointlike, flavor carrying, entities inside the extended constitutent quarks. For 
sufficiently heavy quarks the "cloud" of light degrees of freedom (uu, dd, and possibly ss and gluons) are universal 
and flavor-independent. However, the clouds around s and u quarks may differ significantly. Hence the above 
symmetry assumption cannot be justified in this case when one quark is heavy, and indeed 2m,K* < nrip ~\- m^, which 
relies on an s^°' u'-"-' exchange, does marginally fail. 

Clearly the c u and b u exchanges are even less justified, yet since the hyperfine mass splittings are very small 
(ss 1/toc or 1/mf,), the naive, spin-dependent, potential model derivation of the inequalities holds. 

The rigor - and corresponding paucity of the euclidean correlation function QCD inequalities - reflects, in our view, 
the lack of intuition as to which ^^(a;) field configurations are more important in the functional path integral.^ This 
makes justifying the approximation of keeping certain field configurations more difficult, and is more of a handicap than 
a virtue. By developing such an intuition akin to that inspired by the quark model or strong coupling approximation 



The celebrated example of the extended yet most useful Cooper pair quasiparticles indicates that this may not be a necessary 
condition. 

^''This is clearly a subjective statement, as such intuition might have been developed by lattice gauge calculations. Indeed 
the positivity of contributions of field configurations to the pseudoscalar propagators was known to practitioners in this field 
prior to the advent of QCD inequalities. 
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(which suggested the more important components of wave function(al)s in the variational approach), we could obtain 
many more heuristic, yet very important and useful inequalities. In delineating directions for future research on the 
subject of this review, an obvious goal is to try and find more rigorous proofs for likely correct inequalities, whose 
present derivations are lacking. These inequalities would include the detailed, flavor-dependent baryon-mcson mass 
inequalities; relations between masses of different mesons such as mai(l^) > mp(l"); and mass relations for glueballs 
such as JTig^^ < (mass of any glueball) . 

We believe, however, that there are additional, promising, richer avenues for further research. These involve two 
general (although unfortunately not completely well defined) conjectures that we would like to make next, and the 
application of the QCD inequality techniques to observables other than hadron masses. 



20.1. A conjecture on the Nf dependence of the QCD inequalities 



We have seen that many of the QCD inequalities such as > ttItt and m^^ ■* > m'"^^ ^ are related to observed 
symmetry patterns, namely the spontaneously broken global axial and the conserved vectorial symmetries. The 
correlator inequalities hold for arbitrary Nf. Indeed for Nf degenerate quarks, Nf only appears in the positive 
determinantal factor via [Det(^ + m)]^^ . In particular, all the inequalities hold in the "quenched approximation" 
in which Nf — and the determinantal factor disappears altogether. It is well known that in the quenched limit, 
lattice calculations of QCD, such as hadronic masses, couplings, etc. simplify enormously and one has (within this 
approximation!) completely reliable results. 

The QCD inequalities, for example m^r > to^ and nip > m^, are often satisfied by a large margin. If we had, 
however, a systematic trend for the inequalities to become weaker in the sense, say, that m-^ jrfiT^, rfiplm-^ decrease 
as N f increases, that could furnish very useful information such as 



TOJV 



cxpt "^^r 



(20.1) 

quenched 



We would like to conjecture that this is indeed the case. Our motivation is the likely decrease of the "strength" of 
the SxSB as manifested via the magnitude of (gg) (measured in an appropriate way) with N f. Indeed increasing N f 
weakens - via enhanced screening - the qq gluon exchange interactions and consequently the expected value of (gg) . 

Consider (V'V')a = (trS'F)A, with {'tl^ip) a evaluated for a particular gauge background A [here {0)a = 
J dfi{A){0)A/Z]- We can write 

{trSF)A = V TTT-T-n T = TtTT? 9 > (20.2) 

where in the last expression we paired together the 75 conjugated, non-vanishing eigenvalues ±iAj of p. For the 
m — > limit of interest, the Lorentzians become (5- functions and hence we have 

{qq) ^0^ limp(A) ^0, 

A— !-U 

with p(A) the density obtained after the ^ / A^i{A) averaging of the eigenvalues of ^. 

If d/io(^) = / d[^p(2:)] exp [—Sym{Ap)\ is the Nf = measure, the conjectured decrease of when Nf = ^ 
Nf ^ amounts to: 

/l[[Xf{A)+m']Y^[X'^{A)+m']-'\ < (l[[Xf{A) + n,']\ (j^i^'M) + ^T") > (20-3) 

with (/)o indicating averaging with the dfio{A) measure. Since the factors Yiii^f ~^''^^) ^^"^ X^jl-^j +"^^)^^ appearing 
in the conjectured inequality (|20.3|) are, respectively, monotonically increasing (decreasing) with each A, the conjecture 
is very suggestive. Indeed for any two positive functions of one variable f{x) and g(x) which are monotonically 
increasing (decreasing), {f){g) > (fg)- However, as pointed out to us by O. Kenneth [^12| , this inequality does not 
generally hold for functions of many variables. In t he pr esent context it implies that measures d^(Ai, . . . , Xn) can 
be constructed for which the conjectured inequality ( |20.3| ) can be reversed. This inequality is therefore heuristic and 
depends on additional assumptions. Recalling our comment on functions of one variable, the inequality would apply if 
there were one dominant variable in the measure d/x [^(a;)]- Color confinement and asymptotic freedom suggest that 
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the overall scale R of the {x) fluctuations could serve in such a role. The enhanced weight of larger fluctuations 
coupled with the expected smallness of the corresponding Xi{A) indicates that this may indeed be the case. 

Changing Nf changes the /3 function and hence, for fixed g^, also Aqcd- This in turn implies that the mass/length 
scales used in two lattice calculations with different Nf values should be appropriately changed to allow for meaningful 
comparisons. 

In general hadron masses [and also y/a, with a the string constant (w the coefficient of the linear potential in heavy 
quarkonia, and (qq)^^^)] scale linearly with Aqcd and comparing ratios of such quantities for different Nj values is 
straightforward. However, m7r+ ~ V/tt^tt ~ ^AqcdWq, with mo the {u + d)/2 bare quark mass. Hence instead of 



Eq. (20.1) we should use 

m ivr f _ m f_ 

(20.4) 



quenched 

In conclusion we recall yet another heuristic supporting argument for the conjectured decrease of (qq) / ( Agcp )'^, 
i.e. of the quark condensate, with Nf. It is the restoration of Qs symmetries in exactly SUSY QCD theories ]213| ] at 
moderate Nf /Nc ratios, for which asymptotic freedom (and even confinement) may still hold. Since real QCD is not 
supersymmetric, and further, since QCD inequalities may not apply for SUSY theories, the significance of this fact 
in the present context is not clear. 

20.2. Conjectured inequalities related to the "ferromagnetic" nature of the QCD action 

The monotonic pion (configuration space) wave function was motivated in the conclusion of Sec. ^ above by the 
ferromagnetic nature of the Yang-Mi lls a ction. In some ferromagnetic spin systems this feature is embodied in a 
rigorous set of Griffith's inequalities |214| for the spin correlations. The analogous correlations are not gauge 
invariant but some gauge invariant versions can be conjectured. 

The ferromagnetic character of QCD can motivate many other conjectured inequalities, of which we will mention 
just a few. An (a;) excitation at x (with 1 and 2 spatial indices of E"", the color electric field) and an EfES^d/) 

excitation at correspond intuitively to parallel "spins" (the corresponding plaquettes at x and y in the euclidean 
lattice are indeed parallel). Likewise -Bf-BKx) and B^Bi^iy) are "parallel". However, E^[x)B2{x) and _Bf(y)i3"(y) 
are not "parallel" . The preference of parallel configurations and the stronger positive correlation between such 
configurations suggests therefore that 

{El{y)E^2{y)E'l {x)E^^ [x)) > {El{y)B^,{y)E^^ {x)B'( [x)) 
{B'^{y)B^{y)B'^' {x)B-' (x)) > {El{y)B^,{y)E-^ {x)B'i [x)) , 



{EiE2{x)EiE2{y)){BiB2{x)BiB2{y)) > \{EiB2{x)E2Bi{y))\^ . (20.5) 



and in particular 



which has been suggested by Muzinich and Nair |2l|Jl52( | along with many other inequalities. Since E1E2 (or B1B2) 
and E1B2 acting on the vacuum creates 2++ (2"^) states, the last inequalities suggest that the lowest lying 2++ 
glueball is lighter than the lowest 2 state: 



m^2+^ < m^°_V . (20.6) 



We believe that Eq. (20.5), many other Muzinich-Nair relations (some of which involve J dependence of masses as 
well), and other yet to be discovered relations are true, but also furnish tests not just of the measure positivity in 
QCD, but rather of the ferromagnetic nature of its action. 



'^'^i.e. the state of the system obtained immediately after E^{x)E2(x) operates on the vacuum. Here a = 1 . . . 8 is a color 
index, and 1 and 2 indicate the x\ and X2 components (iJi — etc.). 
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20.3. Inequalities for quantities other than hadronic masses 



The inequalities between two-point correlation functions and the variational Haniiltonian techniques naturally tend 
to yield inequalities between masses, rather than between other hadronic observables. We have seen, however, in 
Sees. and and at the end of Sec. other applications of the QCD inequalities techniques involving scattering 
lengths, form factors, and weak matrix elements. In this last subsection we would like to consider yet another 
quantity, namely high energy hadronic cross sections. These attracted much attention in the sixties and seventies, in 
connection with "Regge pole" exchanges in the crossed ^-channel, and the approximate "quark counting" suggested 
relation cr^jv ~ f fATAr was one of the early indicators for the relevance of the quark model. 

In retrospect asymptotic cross sections may reflect physics which is quite different from that controlling meson and 
baryon masses. Hence while the quark counting rule also suggested ms ^ ^fnnj, which has been "transformed" here 
to the baryon-meson mass inequalities, it is not clear that we have analogous inequalities for hadronic cross sections: 

3 

<^NN ^ -^O-MN ■ (20.7) 

However, this is precisely what we would like to conjecture here (along with the more detailed flavor dependent 
variants). We have attempted to m otivate such a conjecture by using intersecting chromoelectric fluxons as a model 
for high energy cross sections |215| and the stringy trial wave functionals of Sec. ^ While we have not succeeded 
(mainly because of "shadowing" - multiple intersections which are more prominent in baryon-baryon collisions) , this 
model is highly inadequate as it does not account for the rising cross sections. 

We note that if J^^ Ef d^x and J^^ i?| d^a; roughly reflect the masses of hadrons 1 and 2 extending over fii, 5^2, 

and J El ■ E2 d^x generates the Born scattering amplitude, we expect some relation of the formp^ 

C 

0-12 ^ ji mim2 

^^QCD 

Clearly the case of the Goldstone, almost massless pion is again very special.^ The arguments of Sec. |l^ could be 
formally extended to the (non-planar) exchanges pertinent to high energy cross sections. This would then suggest 
relations like {(JTnr)io'NN) > (ctu-at)^, which unfortunately would be difficult to test. 

Physics similar to that motivating the QCD inequalities could be used at many other length scales, and not only for 



underlying composite models. Thus some variant of the Schwartz propagator inequality (12.4) may have implications 
for the frequency of occurence of like-sex non-identical twins, and the more heuristic interflavor inequalities of Sec. |^, 
can suggest many inequalities between bindings of polar molecules. We discuss these applications in App. 
We believe that much more can and will be done on the subject of inequalities, both in and out of QCD. 
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APPENDIX A: LIEB'S COUNTEREXAMPLE TO EQ. (4.1) 

Following Lieb [p^ we prove that for the case of m = oo and two-body potentials which arc infinite square wells: 

V^.in - r,) = U..,...(r) (Al) 
the conjectured "convexity" relation 

m,m,M) (A2) 

is violated. 

To achieve minimal (and indeed even finite!) ground state energies, the wave functions have to vanish when 
— ?^"| > '"o for any quark pair. For heavy quarks this can be done with no kinetic energy penalty. Thus let us 
fix ri = in all three wave functions ^p'^^\m,m,m),4'^^\m,m, M), and ip''^\ra^M,M), thereby disposing also of the 
overall translational degrees of freedom. For ^^^"^(171,171,171) we could take all fi — 0, achieving the minimal possible 
energy: E'^^\m, m, m) = 0. In ip'-^'' (m, m, M) it is clearly advantageous to put also r2 = = so that the condition 
1^3 ~ ^2! < ''o is satisfied automatically once |r*3 — rij < r^. Thus finding E'^^^ {rn, m, M) reduces to finding the minimal 
energy of the one-body Hamiltonian p'^/2M + 2V{r), which for V = Vi,s.w. is h = p^/2M + V{r). 

On the other hand, E'^^\m, M, M) is the minimal energy of 

h = [pi/2M + U(r2)] + [pI/2M + V{r^)\ + V{r2 - r^) ■ (A3) 
If V{f2 — fa) were absent, h separates into two one-body infinite square well problems and 

E'-"Hh) = 2E^^\h) . 

However, since (V23) > we have 

{h) = E^^^ (m, M, M) > 2E^°"> {h) = 2E^°^ (m, m, M) , (A4) 

and the desired inequality ( |A2[ ) is violated. Indeed V'mMAfl'^i ~ 01^2,7^3) has to vanish now in all of the region 
1^2 — J^l > ''o in addition to the vanishing for ^2 > ro, ra > rp. Incorporating this extra constraint reduces the allowed 
six-dimensional volume in (r2,?'3) space from (47rrQ/3)^ to a fraction thereof. Obviously this will increase the kinetic 
energy from 2E^^\h) to (2 -I- a)E'^°'>{h), with awl. 

By continuity we therefore expect Eq. ( |A2| ) to fail already when we approach the limit m 00 and V{r) —> Vi,s.w.i^)j 
e.g. via 

y(r) =c„(r/ro)", n ^ 00 . (A5) 

In passing we note that exp [—PVi,s.w. {''')] — ©(^^^o) is indeed not a postive semidcfinite operator, since its Fourier 
transform is not positive. This should be the case since we have, as indicated in the following Appendix, Lieb's 
theorem that the desired inequality holds when exp \—0V] is positive semidefinite. Also for V{r) = r", n > 4, a case 
for which another counter-example was produced first |39| , exp [— /3U(r)] is not positive semidefinite - consistent with 
the fact that it fails the convexity conditions. 
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APPENDIX B: DISCUSSION OF LIEB'S RESULTS FOR THREE-BODY HAMILTONIANS 



Following Lieb [|18[, we prove the inequality (|A2|) when the potentials are flavor independent: 

V{m, m) = V{M, M) = V{m, M) = V , 

and exp [—(3V{x — y)] positive semidefinite. 

The propagation over (imaginary) euclidean time /? of the three-body system is given by (X|e~'^^|X'), with X — 
{xi,X2, xzjiniUauX' = {xi , X2, X3) f.nai ■ We compaie three systems with 

Ha = Ti{xi) + T2{X2) + T'2(.T3) + V2{xi,X2) + V2{xi,Xz) + Vi{x2,Xz) 

Hb = Ti{xi) + T2{x2) + Tz{xz) + Vs.{xi,X2) + V2{xi,xz) + Vi{x2,xz) 

H, = Ti{xi)+T:i{x2)+T3{:xz) + V:i{xi,X2) + V:i{xi,xz) + Vi{x2,X3) . (Bl) 
For X = X' , Zp{x) = {X\e^P^\X) is dominated, for (3 oo, by the lowest energy state 

(Xle-f^^lX) = Y,\{X\n)f e-^^" pZL |(X|0)|' e"'^^" , (B2) 

where we use the completeness sum over energy eigenstates. 
Thus to prove the desired inequality 

it suffices to show that 

Z0{a)Zpic) > Zf,{b) . (B3) 

In general the path integral, and Zp in particular, is obtained by dividing the total evolution into many consecutive 
evolutions over small time steps, so that 



Zf3 = hm ZpiN) 

N^OO 



N 



with T and V the total one-body kinetic and two-body potential parts of the Hamiltonian H = T + V . Clearly it is 
sufficient to prove the inequality for each N . 

We next insert a complete set of X-space states between each pair of c^l^"^/^ q-0v/n fg^gf^Qj-g This gives Zp as a 
path integral. Each such path consists of three (spatially) closed polygonal paths Xi,X2,X3, with Xi consisting of 
the iV 4- 1 points 

Xi(0) =Xi,Xi(/3/Ar),Xi(2/3/7V),...,Xi(Ar/3/iV) =Xi(/3) =Xi , 
and likewise for X2 and Xy, (see Fig. ^l|). Since each factor 

{X,{j)X2{])X3{j)KP^"' e-''^/^|Xi(j + l)X2{] + l)X3(j + 1)) (B4) 



involves an evolution over an infinitesimal "time" P/N , the non-commutativity of the kinetic and potential parts of 

7N 



H is neglected. The kinetic single particle operators contribute to Zg (a) a factor of 



and to Z^{b) and Z^(c) 



Fa ^ Fi{Xi)F2{X2)F2{X:,) , 

Fb = F,iX,)F2{X2)F3iX3) 
F, ^ Fi{Xi)F3iX2)F3{X3) 



with Fi{X) corresponding to Ti in Eq. (Bl) above. 

The potential two-body operators contribute to Zp' a product of three terms depending on the three pairs of paths. 
Specifically, the contributions in cases (a), (b), and (c) are 
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{^It^2)G^ {Xi,X'i)Gi {X2tXz) , 



-0 



where the Gi{Xk, Xi) are generated from Vi{Xk,Xi) via 

AT 

G,{Xk,Xi) = l[cxp[~{f3/N)mXk{f3j/N),Xi{/3j/N))] . 

For the particular case of Vi, the fact that e~^'^'-^-^ is positive semidefinite ensures that the iV-fold tensor product 
defining Gi{X2, X3) is also positive semidefinite. Collecting all terms and separating the d^^^~^-'Xi integrations over 
the — 1 intermediate points X^, . . . , X^ along the polygonal path Xi we have 



Z^{a)^J S^^-'^X,F,{X,) J d'^^-^^X2 j d3(^-i)X3 

X G^{Xi,X2)F2{X2)G^{Xi,X3)F2{X3)G^ {X2,X3) 
Z^{b) - / d3(^-l)Xi^^i(Xi) / d=^(^-l)X2 / d3(^-l)X3 



X G^{Xi,X2)F2{X2)G^{Xi,X3)F:i{Xs)G^ {X2,X3) 

Z^{c)=J d3(^-l)Xi^^i(li) J d^^^-'^X2 J d3(^-l)X3 

X G^{X,,X2)F3{X2)G^{Xi,X3)F3{X3)G^ {X2,X3) . (B5) 
Defining next the N x N matrices 



V^^Xi,X2) = x^F,{Xi)G^iX,,X2)F2{X2) 



V^3^(Xi,X2) = JF,{X,)G^{X,,X2)F3iX2) , (B6) 



we can write Eq. (B5) in the concise form 



Z^{a)^triV2^ ■G,-V2) 



Z^(6)=tr(F/.Gi-F3) 

Z.^(c)=tr(F3^.Gi-F3). (B7) 



We note that the positivity of Fi{Xi) was implicitly assumed in taking the square root in Eq. ( p36|) . For the case 
of interest with Ti — pf/{2m.i), {Xi\e^^'^\Xi) is the probability for returning to the initial point in a Gaussian 
random walk after a time /3. F^ (Xi) is the probability that this happens for a specific path Xi with — 1 specific 
intermediate steps, and hence is clearly positive. This positivity of the "hea t kernel" can apparently be generalized 
to other relativistic forms of the kinetic ener gy, for example Ti = \/ + . 



To complete the proof we note that Eqs. (B7) define a "scalar product" of V2 with itself, V2 and V3, and V3 with 
itself. This can be most clearly seen by transforming to the basis in which Gi is diagonal with all diagonal elements 
real and positive. This change of basis leaves tr (V^GiV^) invariant but casts it in the form 



J2iV2)an{Gl)nn{V3) r. 



which is bilinear in V2 and V3 and positive whenever two identical "F-vectors" ( of length [3(A^ — f)]'^ ) are used. 
Hence these scalar products satisfy the Schwartz inequality which is precisely the desired result ( |B3[ ). 

It has been noted by Lieb that the assumption of having only two-body potentials is rather restrictive and we 
can allow also genuine three-body potentials V{Xi, X2, X^) as long as e~^^ is positive semidefinite as a function of 
X2 , X3 . 
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APPENDIX C: PROOF OF EQ. (5.2) 



We prove Eq. (5.2) using a concrete basis of states in which it holds for all matrix elements. We will use the 
Hamiltonian version of QCD with a spati al cu bic lattic (sites n, unit vectors fj along the positive x,?/, z axes). 
The QCD Hamiltonian is written as [ ^217| ] 



^QCD = y E ^<En,fi+f,f ^^(UUU^U^ + h.c.) 

n , ?) > plaquct tcs 

+ E E ijl{n)Un,n+fiA{n + fj) + h.c. + YYmiilji{n)-ilji{n) ^ . (CI) 

i=l n,rj>0 i=l n 

It depends on the SU(3) "connection" matrices C/f;,n+?j which live on lattice links between n and n + fi, the canonically 
conjugate generators Efi fi^fj, and the fermionic spinors ipi{n) [and conjugate V'i(n)] at each lattice site. The first 
three terms are analagous to the £"^,8^ and ^(a:) pip{x) terms in the continuum Hamiltonian. 

We adopt the "strong coupling" basis in which all tr(£'.,-j^fi+jj)^ are diagonal. At each lattice site we specify the 
spinors ipi{n). The mass term is 

Nf 

J2J2'^Uf^homMn) . (C2) 

n i—1 

Each of the plaquette terms creates a minimal, closed, E flux line going clockwise (or anti-clockwise) around an 
elementary plaquette. This will, in general, change tr(i?fj 3+^,)^ on each of the four adjoining links and hence gives 
rise to non-diagonal elements in the representation considered here. Also the third ['ip{x)pip{x)] term changes the 
lattice configuration by allowing the quark to "hop" to a neighboring site (dragging along a flux line) , and by creating 
or annihilating qiqi pairs at neighboring sites. 

We illustrate some of the matrix elements of Hqcb on a small 3x3 two-dimensional lattice in Fig. |^. Each of the 
columns and rows is labeled by a complete lattice conflguration. In the matrix we indicate both the actual value of 
the matrix elements (above dotted line) and the specific terms of iJqcD contributing to it (below dotted line). We 
omit additional four-fold spinor indices for each of the occupied sites. 

Local gauge invariance 

V'j(n) V{n)^pi{n) 
Un,n+fi V{n)Uri,f{+fiV^ri + ff) 

En.n+f, ^ V{n)En,-i+ny^ (n + v) , (C3) 
is respected by all terms of -ffqcD- It constrains the physically allowed states via Gauss' equation 

stating that the sum of all outgoing Efi fi^jj flux lines [which generate T^(?i)] vanishes at all lattice sites n except those 
n with external color charges Qi. (The Qi are due to quarks. Up to 2Nf quarks and 2Nf antiquarks are allowed at 
each site by Fermi statistics.) 

It is precisely via this Gauss condition that the specific (mesonic, baryonic etc.) sector becomes relevant. Figure 
^ pertains to the M„j sector. The simple configurations in Fig. ^ have just one u and one d in the lattice. Eq. ( |C^ ) 
then simplifies into 

J2 ^n^n+v = ^i^' ^u)X + 5{r^ rd)At , (C5) 

where r„(j-) are the locations of u((i), and A, are the corresponding color matrices. 

Since Hqcu commutes with the gauge transformations, the evolution of the state on the lattice, allowing for creation 
of qq pairs at one lattice point and their subsequent separation and possible annihilation with other qq, maintains the 
Gauss condition (p4) [though not, in general, its simple "valence" version (pq)]. Since the Gauss l aw c onstraints are 



common to all the mesonic sectors, they do not interfere with the suggested operator relation Eq. (5.2). 

We start by showing that ( ^.2| ) is satisfied when we switch off the quark creation and annihilation terms and then 
work our way gradually to the general case. 
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1. In the simple case where consists of just the qigj quarks and not other pairs, it it easy to see that all the 
non-mass terms in Hqcb have identical matrix elements in any Hij sector. Also the Gauss law constraints are the 
same. The matrix elements depend only on the C/'s and -E's and the generic quarks or antiquarks but not on the 
specific quark flavors. Taking for concreteness ijkl — udsc, we have flavor dependence only in the mass terms ('0'^ 
denotes a charge conjugated spinor of an antiquark): 



^mass 
ud 




hmdi^lJ^y{fq)-fo^^{rq) 


(C6a) 


^mass 


= mcip^rq)jQip{fq) - 


^msii;^yirq)jo^^{rq) 


(C6b) 


rrinass 

us 


= rriuij'' {rq)'jotp(rq) - 


f msiij^)'' {fq)'j0lp'^ {fq) 


(C6c) 


rrinass 


= mctp'' {rq)jQ^{rq) 4 


- TUditp^y {rqhoiJ^ {rq). 


(C6d) 



Since the same four terms appear in H^^ + ^cs as in Hag + H^^^ Eq. (|5.2| ) is valid in this valence approximation: 

Hus + -ffcd — -^ud + (subspace with two quarks) . (C7) 

2. The above argument is not affected by the creation/annihilation terms of all pairs xx uu, dd, ss, cc. The created 
quarks or antiquarks cannot be Pauli blocked or annihilated by the valence quarks or antiquarks. Their effect, just 
like that of gauge flelds, is common to Hcs, etc. 

3. Finally, we turn on also pair creation/annihilation of the valence flavors. We would then generate configurations 
such as q^{uiUi) . . . {xpXp)q^ in which ui is at the same lattice site, and same spin and color state as the valence quark 
g". If = u, such configurations should, by the Pauli principle, be disallowed. Also ui could subsequently annihilate 
the valence quark. Both effects occur in M^^ and ^^us but not in Mcs and M^^- To avoid a possible difficulty in 
deriving Eq. ( |5.2| ) we will define the configuration corresponding to 

Uv{uiui) . . . {xpXp)dy in M„j [or Uy{uiUi) . . . {xpXp)sv in M„j] (C8) 

to be 

c„(cici) . . . {xpXp)dy in M^^ [or Cv{ciCi) . . . {xpXp)sy in Mcs] ■ (C9) 

This, rather than the definition with a common set of pairs, ensures that configurations excluded in Af^j or M^s 
will be excluded in Af^j or Mcs', and that for any subsequent uiu^ annihilation induced by H^^ there will 

be a corresponding annihilation of ciCy in iJ^j or Hcs- The diagonal mass terms which contributed to + Hcs 
and iJ^j + Hys in such configurations will still be the same. In each case we have Smu, 3mc, rud and nis terms and a 
common set of non-wscd quark masses due to other pair creation. 

The construction of corresponding configurations is readily generalized to the case of arbitrary number of pairs with 
the flavors of the valence quark. A generic state in Af^ is 

g,(Q/^)g-/^)) . . . (g^)g-^))(x^)-(gf . . . iqf^kj'^% • (CIO) 

The superscripts py, n^, ■ ■ ■ simply count the number of pairs of type (qiqi), {x, x), etc. Altogether we have py pairs of 
t he valence quark flavor, py pairs of the valence antiquark flavor, and additional xx pairs. With the specific relation 
( ^.2| ) in mind we will explicitly exhibit among those xx pairs the {ququ) and [qiqi) pairs: 

A/.,- = q^{q^q^^- {q,q,f^ qM^kf' (ffi^O^^ (a^^)"- ■ ■ • (CH) 
with X 7^ qi,qj,qkiqi- The corresponding states in the other three mesonic sectors will be taken as 

Mki = gfc(gfeft)''"(gig/)''"gifeg»)''^(gjgj)''^(xa;)"== . . . 
Ma = g.fegO^^ta^r^ife^fcrfe^jfH^^s)"^ • ■ • 

Mkj = gfe(gfegfe)''^fegi)P^"g,(ft%)^n9;ffir(2^s)"^ ■ • ■ (C12) 

All of these configurations have the same total number of quarks and antiquarks: iV = 2 + 2py + 2py + 2pi + 2p2 + 
2j2x^ijki''^x- Also the corresponding quarks (or antiquarks) in Mij,Mj.i, etc. are taken to be at the same locations 
with identical spinor and color states. The Gauss constraints, the antisymmetrization effects, and (except for mass 
contributions) also all matrix elements are the same for , Hj^, Hq and H^j. Finally, we have in Hij + Hf^ the same 
mass terms as in Ha + H^-j. [Altogether there are 2pi + 2py + 1 terms with rrii as a coefficient, the same number with 
coefficient m^; an d 2p 2 + 2p„ + 1 terms with coeffecients mj{mi).] All of these terms make identical contributions to 
both sides of Eq. ( |3.2[ ), which is therefore correct. 
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Viewing states with equal of numbers of (say) cc and uu pairs as "corresponding" states does not imply that heavy 
and low mass pairs occur with equal probability in physical states. This issue is determined by solving H^'^^'^\^ij) = 

]. This is analagous to the relative separation rqq in 
the example of Sec. ^, where r^c and run correspond to the same "generic" degree of freedom r. While solving the 
Schrodinger equations with TOc ^ leads to very different expectation values of r 

the operator relation H^c + Huu — 2iJc« still holds for the nonrelativistic Hamiltonians. 



-^'■^"'^ l^ii) s-i^d heavy pairs are in fact strongly suppressed [218 



APPENDIX D: ENUMERATION AND SPECIFICATION OF THE BARYON-MESON INEQUALITIES 

Altogether we have forty two different flavor-spin baryonic "ground state" combinations consisting of udscb quarks, 
when we do not distinguish between members of /-spin niultiplets and assume no orbital, radial, or gluonic excitations. 
(Electromagnetic mass splittings were discussed in Sec. p^ ). 

We consider first the — (3/2)+ "A-like" baryons. In baryons made of three identical quarks bbb, ccc, sss {^^), 
and uuu = ddd (A++, A^), each diquark combination is in a spin triplet with = (3/2)+. Also, in all other (3/2)+ 
states all quark pair subsystems, and hence the corresponding mesonic systems qiqj{-^ Qi(lj),1jQk{-^ QjQk), Qk^ii-^ 
QkQi), must be in the triplet, 5 = 1 state. Hence only vector meson masses appear on the right hand side of the 
inequalities. If we list all mesons in lexicographic order with b > c > s > u,d, we have the following (3/2)+ states: 
uuu (A), suu (S*), ssu (S*), sss {^~), cuu (S^); csM (S*), CSS (ri*); the doubly charmed baryons ecu and ccs; the ccc 
state; and buu (EJ), bsu (SJ), and bss (fi^). There are also the heavier states bcc, bbu, bbs, bbc, and finally bbb, which 
are unlikely to be discovered soon. 

For the sake of completeness, we list in Table || all the relevant inequalities. The first five can already be tested 
and are satisfied with a reasonable margin in all cases. The next eight inequalities constitute lower bounds predicted 
for as yet undiscovered baryonic states, which hopefully can be verified in the near future. The remaining three will 
presumably be verified only in the distant future. 

We move next to the = (1/2)+ baryons. Let us first focus on the case where we have two identical quarks, 
namely duu (jj), suu (or sud with / = 1) (S+), ssu (S^), cuu (S++), ess (Sc), ecu, ecs, buu, bss, and bee. In all these 
^^yxiX2 - type" states the identical flavor quarks xiX2 must be in a triplet state, namely 

(s.i + sx,f = sl^,^ = 1 . (1 + 1) = 2 . (Dl) 

Writing next the total baryon spin as 

{Sx, + S + {svY - 4 (D2a) 



or, using Eq. (Dl) and Sg = 1/2 



2 -I- 2(5;, • Sy + Sx, ■ Sy) + 3/4 = 3/4, (D2b) 

i.e. 

{^xi ' ^y) B (^^2 ' "^y) B ~ (D2c) 

with ( ) B denoting the expectation value in the baryon state. Since we have Xi X2 symmetry of the above 
expectation values, we finally conclude that 

{Sx^■Sy)g = {Sx^■Sy)g = -1/2. (D3) 

Hence each of the yxi {-^ yxt) subsystems must be a mixture of 3/4 singlet and 1/4 triplet states for which {s ■ s) 
= - 3/4 and -f 1/4, respectively. This then implies, by collecting all the terms in {X1X2) ,{yxi) ,{x2y), that the 
baryon- meson inequalities read as follows: 

-(TOp -I- m^) 





> 


rrip 
2 


h ^(3m^ 


+ rup) ^ 




> 


TUp 

2 




+ niK' ) 


ms 


> 


m0 
2 


1- ^{irriK 


+ rriK') 




> 


nip 
2 


h ^(3m£) 


+ niD') 




> 


TUp 

2 


h ^(3mB 


+ ms*) 



(D4) 
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We have similar relations for the doubly strange baryons: 



TTi 1 

ms, = mess (1/2)^ > + ^(3tod, + mo^) 

TTI 1 

ms, = mbss(l/2)+ > ^ + ^(Sms, +771B.) , (D5) 
and the doubly charmed baryons (hopefuUy to soon be found in FNAL experiments): 

mccu{l/2)+ > —^ + ^{3mD+mD') 

m,,,{l/2)+ > + i(3mD, + mo-J ■ (D6) 

The specific numerical values for all of the above and a few other xiX2y states are listed in Table ||. 

Among the remaining (1/2)+ baryonic states composed of three different quark flavors, we have three "A-type" 
states in which the light quark ud subsystem coupled to / = (and hence to Sud = lud = 0). The "S-type" states 
with lud = 1 have already been discussed above. The "A-type" states are 

s = A c ^J^^d^ ^ b .^iJ^ = At, . (D7) 

7=0, s=0 7=0, s=0 7=0, s=0 

Using 

(Sud) = {Su + Sdf = 

and 

(sa) = {sx + SudT = 3/4 , 
we can now deduce, using also u ^ d exchange symmetry, that 

{sx ■ Su) = {sx • Srf) = . (D8) 



Equation (D8) implies that in e ach xu {-^ xu) and xd (-^ xd) subsystem we have the opposite triplet-singlet mixture as 



in the previous case of Eq. (D3), namely 1/4 singlet and 3/4 triplet. Thus collecting these and the term corresponding 



to the singlet ud — > ud (« it) subsystems, we finally have: 

TIT'T! 1 

TOA > — + -{irriK + mK') 

TTi 1 

JTiAe > + -^{imD + mo') 

TTI 1 

rriA^ > + -{Stub + niB') ■ (D9) 

The remaining baryonic states which consist of three different flavor quarks csu, bsu, bcu, and bcs are either of the 
"A-type" in which the lighter two quarks couple to = subsystems, or of the "S-type" in which the lighter two 
quarks couple to 5* = 1. Because the hyperfine interaction [cx (si • Sj)/{mimj)] is bigger for the lighter quark system 
we expect that the A-type states will be lighter than the corresponding S-type states. This will be particularly so 
for the su subsystem in csu or bsu, which have been observed or are likely to be observed soon, than the bcu and 
bcs states. When we specify the spin q of the su subsystem to be Sq-^q^ = 1 (or Sq-^q^ — 0), and follow the previous 
discussion, we conclude that, for csu for example: 

{sc ■ Su) + {sc ■ Ss) = -1 (for Sud = 1) (DlOa) 



{sc ■ Su) + {sc ■ Ss) = (fors„d = 0). (DlOb) 

However, unlike in the previous cases, we cannot deduce that {sq ■ Su) ^ — i^Q ' ^d) g without appealing to the 
approximate Gell-Mann - Ne'eman SU(3) u ^ s symmetry. If we do this nonetheless, as a first approximation, we 
obtain 
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.(A)> 



rriK 



{niD, + mo) 3(mcj + mo') 



(Dlla) 



fn\- 1 



(Dllb) 



In the real case, we expect the spin-averaged variant, consisting of the symmetric combination of Eqs. (Dlla) and 
( DllbD (this is the value we list in Table ||, with the left hand side called rricus)'- 



(A) + ^CUS 



[rriK + mK.) {mD+mD') {nio, + mo'J 



(D12) 



and in general we expect that the right hand sides of (Dlla) and ( Dlllj) to be respec tively an over- (under-) estimate 
of the combined mcsonic masses. Specifically, we expect that Eqs. ( plla ) and ( Dllb ) t o be sa tisfied wit h a som ewhat 
smaller (larger) than usual margin. The difference between the right hand side in Eq. ( |Dlla ) and Eq. ( Dllbj ) is 



mK) 



ISOMeV , 



(D13) 



We will not pursue this issue further here. 



APPENDIX E: APPLICATION TO QUADRONIUM 

It has been sug gested that the puzzling narrow resonance found in heavy ion collisions is in some sense an e+e^e^e" 
bound state [ ^19| . The need to achieve very strong bindings (BE ~ rrig — 0.5 MeV) in this normally weakly (electro- 
magnetically) coupled system is a serious problem of the quadronium hypothesis. While the experimental motivations 
are questioned, we find that there is a nice, simple, application of the variational technique [ p20[ . It has been sug- 
gested that we cannot really have a two-body dominated interaction in this system, and that new, unusual, four-body 
interactions are called for. 

The following inequality on the binding energy of quadronium in terms of the binding of positronium serves to 
sharpen the issue. Let ipg = 'i/'''(li22) be the wave function of quadronium in its ground state. The Hamiltonian is 

Hq=Ti+Ti+T2+T2 + Vn + Vi2 + 1^21 + V22 + Vi2 + V12 . (El) 

Clearly we expect the repulsive electron-electron {V12) and positron-positron (V12) interactions to lower the quadro- 
nium binding. Hence 

e%{Q)<e%{Q). (E2) 

with e%{Q) the binding of a fictitious Hamiltonian Hq, which is the original Hq with these repulsive interactions 
eliminated: 

Hq=Ti+Ti+T2+T2 + + V^2 + V21 + V22 . (E3) 

If we compare Hq with the two-body Hamiltonians, 

H^l=Ti+Ti + Vn, H^2 = Ti+T-2 + V^2 

H21 = 72 + Tj + V21 , H22 =T2 +T2 + V22 

we see that 

2Hq ^ H[j + + + ^22 ' (E4) 

where H'^^ is obtained from H^i by doubling the interaction: H'^^ = Ti + Tj -I- 214j, etc. This can be achieved for the 



one-photon exchange by doubling the effective coupling strength a ^ a' ~ 2 a. Equation (E4) can next be used in 



the by now familiar way to put an upper bound on e]j(Q), and thus, via Eq. (E2), on eg{Q) as well. The bound is 



e%{Q) < 24° (P) , 



65 



with £'g{P) the binding of positronium in which the strength of the interaction has been doubled by a ^ 2a. Thus 
we have 

< 4iQ) < 24° (P) . (E5) 

The fact that the real positronium spectrum conforms so nicely to QED predictions with radiative corrections included 
strongly suggests that if a = 1/137 is scaled up to a' = 2/237 we can still treat the e^e~e~^e~ system with the 
perturbatively calculated potential. In this case 

e'i>(P)^(aTj = 2a2|=2Ry = 27eV, 
and hence e'siQ) < 2e'g{P) < 54 eV, and the binding falls short by about 10^ of the required 0.5 MeV binding. 



APPENDIX F: QCD INEQUALITIES FOR EXOTIC NOVEL HADRON STATES 

Most known hadrons are qigj mesons and qiQjqk baryons. We will refer to (non-glueball) states belonging to neither 
of the above two categories as "exotic" . Such states have been searched for and discussed for more than 30 years. 
"Duality" arguments predating QCD suggested "gggg" states cou pling mainly to baryon-antibaryon as the most 



likely exotic states. Within QCD hybrids {qq + gluon) ar ise naturally ||22l| . Also qqss "bag states" were suggested as 



candidates for the /o and KK threshold states 222| -224|. Considerations of hyperfine chromomagnetic interactions 



generated via one gluon exchange between the quarks qi and qj suggested a particular "hexaquark" Ti. = uud dss 
below the AA threshold which could be strong interaction stable.^ More recently similar considerations [ 227 , p28| 
singled out the specific "pentaquark" V = csuud as a more likely strong inte ract i on s table, new, exotic, bound state. 
Both states have been experimentally searched for with inconclusive results |22£ -231| in almost all cases. 



The exotic states can be separated into disjoint color singlet hadrons.^ Indeed we have seen in Sec. |l^ that for 
degenerate quark flavors (and most likely in other cases as well) there is an attractive interaction in the qiqjqkqi 
channels. The question of whether these residual color forces between color singlet states can lead to weakly bound 
states, in analogy to the nuclear forces in the deuteron and other nuclei, is of some interest. However we believe that 
it is a detailed, delicate, issue that QCD inequalities cannot decide (much in the same way that the deuteron's binding 
requires detailed calculation). 

Here we would like to focus on the question whether genuine multiquark states - referred to in the literature 
qualitatively as "single bag states" - which have significantly higher than nuclear bindings, exist. The arguments for 
bound hexa- and pentaquark states utilized the ad hoc assumption that in all multiquark hadronic states the quarks 
are in the same "universal" bag as the baryons and mesons, and have the same hyperfine interactions. Here we will 
utilize the complementary, strong coupling picture, where the quarks in the exotic states are connected via electric 
fiux lines into a more complex color network than the qq for mesons or the "Y" topology with one junction point for 
the baryons. 

Let us first focus on quadriquarks. In the limit when two quarks (or antiquarks) are heavy, one can show dire ctly 
that a novel pattern of an overall singlet system which cannot be broken into color singlet clusters will form ]232| | . In 
a qiqjQkQi state, the two heavy quarks always bind coulombically into a color triplet {QkQi)i diquark with binding 
and size proportional to mg and Tig^, respectively. The small {QkQi)'i subsystem acts effectively as a heavy quark. 
Together with the light qiqj it will then form, as a Ag analog, the qiqjQkQi quadriquark. Because of the small 

numerical coefficient in the coulombic energy, Eb ~ \ (^)^ ^he latter does not exceed Aqcd — 0.2 GeV even 
for toq — niB — 5 GeV. Hence the question of the stability of the novel pattern of color coupling - which in a strong 
coupling chromoelectric flux tube picture corresponds to the connected color network illustrated in Fig. |l^ - remains 
open. 



^°The existence of a AA bound state could encourage the further conjecture 
exist and that at large densities strange quark matter is the most stable, which if true has fascinating astrophysical ramifications. 
^^With the exception of genuine "quantum number exotics", for example 0* or 1* states that cannot decay into mesons. 



225 226 j that higher 



"strangelets" 
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In the following we would like to use QCD variational inequalities to address this problem. This we do by following 
the derivation in the same strong coupling, flux tube approximation as the baryon-meson mass inequalities in Sec. g 
above. Thus, let us take for the quadriquark wave functional: 

I*99W)=E^>^H , (Fl) 

with the possible generalizations maintaining the two junction points (1) and (2) (see the discussion in Sec. |^). We 
have used ixi to symbolize the graph of Fig. As indicated in Fig. we can extract from {^qq'QQ') - and the 
wave functional of an anti-quadriquark superposed on it with flux lines reversed - trial wave functionals for the four 
mesons QQ', Qq' ^ Qq' ^ and qq' . In the above Hilbert space, we then have the operator relation 

2HQQ'qq> = Hqq, + HQq> + Hqq> + Hg,Q, , (F2) 

from which we obtain, via the variational principle, the mass inequality 

^^QQ'qq' - ^qq' ^ ^Qq' + ^qq' + ^Q'q' ' 

For the case Q,Q' — c, qq' — u, d, this binds from below the (hypothetical) quadronium mass by the known charmo- 
nium, charmed meson, and light meson masses. Specifically for ccuu (or ccdd), the inequality reads: 

TTT-^cla > \ [mj/^ + rrip + (1 + a)m£,. + (1 - a)m£,] , (F4) 

where nice = nT-j/tp and niuu = follows from the generalized Pauli principle as discussed previously The weights 
of the D and D* depend on the total spin of the ccuu system (see App. ^). For example, consider — 1+ quantum 
numbers. In this case the transition quadriquark DD is forbidden since the two D mesons can couple only to 
0+, 1~, 2+ etc. The lowest state available for decay is then DD* . Performing a calculation similar to those in App. ^ 
we find that a = is the correct weight and the inequality reads 

m-i°L(l'*') > ^ {mj/^ +mp + mo* + mo) 

"^SL(1+) > 3870. 7MeV w mo' + mo ■ (F5) 

Our previous comparisons of many baryon-meson inequalities with data indicate that these inequalities are satisfied 
with a margin of 150-300 MeV. Cohen and Lipkin and Imbo suggest heuristic arguments for this margin. 
Such a margin should a fortiori be valid for the quadriquark inequalities. Thus we expect that the ccuu state 
actually lies well above the D* D threshold and therefore is unstable. The situation is very different for ccud. In this 
case we could choose ud to be in the spin singlet state. This is consistent with the cc spin triplet diquark and the ud 
spin si ngle t anti-diquark, with relative zero orbital angular momenta, adding up to an overall = 1+. The analog 



of Eq. (F5) is then 

(o)/+\l/ 3 1\ 

m^°^^-(l+) > 3590.6 MeV w mD* + mo - 300 MeV , (F6) 

A bound ccud{l~^) state with rn^^^_^ < mjj + would satisfy the inequality with a margin of approximately 300 
MeV, and is therefore quite likely. 
The analogs in the Q = b system 

m["],^(l+) > i (mx + mp + tub* + ms) 

™il!L(l+)> 10416.3 MeV, (F7) 



The identical uu or cc quarks are in the lowest L — state, and hence the uu or cc in the corresponding meson should be 
in the triplet state. 

^^Here we have for some of the Qq mesons worse trial wave functionals as compared with those extracted from baryons, 
namely mesonic strings consisting of three segments going through the two junction points. 
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and 



m'Jad^r) - 2 i ™^ + ™7r + ^mB' + -TUB 
m[°|^,-(l+)> 10113.3 MeV (F8) 

allow bindings in both flavor combinations, though again with a considerably higher safety margin in the bbud case. 

In baryons each quark pair must couple to a 3. In the QQ'qq' system the QQ' and qq' could also couple to a color 
sextet. Thus the segment connecting the junction points (1) and (2) in Fig. |l^ would then carry a sextet flux - a 
possibility that our discussion above neglected. This is justified in the strong coupling limit. The potential energy 

J d'^xE^ then dominates, and the quadratic Casimir operator (to which J E"^ is proportional) of the sextet is 10/3 
times larger than that of 3. Nonetheless, the need to exclude the sextet flux, and the related fact that there is no 
alternative derivation in the weak coupling, one gluon exchange limit, puts the present inequalities on a weaker footing 
than that of the baryon-meson mass inequalities. We believe however the enhanced likelihood for a ccud quadriquark 
state suggested by the inequalities. 

It is amusing to note that an alternative, purely hadronic approach to the problem of a bound DD* state also 
suggests that D^D* (of a ccud composition) is more likely to be bound than D'^D^* . Thus following Tornqvist 
| p33[ (who coined the name "deusons" for these deuteron-hke extended mesonic states), let us consider the one pion 
exchange interaction in the D*D channel. The resulting Yukawa force has, due to the Dtt threshold ~ D* proximity 
(i.e. rriB)* ~ mB — m-r; = e <C m-^), an anomalously large range ~ \/\/2ern^ ~ 3 — 7 fmp| and could generate a 
D*D bound state. However, because of simple /-spin Clebsch-Gordan coefficients, the potential generated by the 7r+ 
exchange in the £)*+_D° <-> _D*°Z?+ system is twice as strong as that due to the 7r° exchange in the D*^D^ ^ D*^D^ 
system, and a D*^D^ bound state is more likely. 

Let us next consider pentaquarks in the same strong coupling flux string approximation. The corresponding novel 
connected color network is illustrated in Fig. This "network" contains three junction points. At point (1), two 
light quark fluxes are coupled to a 3 flux, and at point (2) the same holds for the fluxes emanating from the other two 
light quarks. Finally at the third "anti-junction" point the three 3 fluxes originating from points (1), (2), and from 
the heavy antiquark Q, couple to a singlet. As in the case of the quadriquark we could have coupled the fluxons in 
either the first vertex (1) or the second vertex (2) to a color sextet fiux, which again is neglected in the strong coupling 
limit. For pentaquarks, we furthermore have an alternative (similar) network obtained by exchanging one of the u 
quarks from the uu in vertex (1) with the d quark in vertex (2) - which represents a different, though not necessarily 
orthogonal, color coupling configuration. By reversing the flux lines on the segments (2)-d and (3)-Q indicated by the 
double arrows in Fig. |2^, and ignoring the intermediate fluxon between (2) and (3), we can extract from a "stringy" 
pentaquark wave functional trial wave functionals for a charmed baryon [cuu in this speciflc case) and an sd meson. 

The full Hamiltonian for the quark and fluxon system is 

quark flux segment 

where are the Dirac Hamiltonians for the quark, {E)'^ represents the "potential energy" density for the fluxons, 
and plays the role of kinetic energy for moving and distorting the string segments and the in i/^). "moves" 
the tip quarks. 

Omitting (the generally positive contribution of) the (2)- (3) flux segment will only reduce the energy. This and the 
fact that the trial conflgurations of the baryons and mesons thus extracted do not optimize the wave functionals of 
the meson and baryon at rest, strongly suggest the inequality 

m-p > m^^ + nicuu ■ (FIO) 

Because of the Pauli principle the two u quarks in cuu which are color antisymmetric and flavor symmetric must be 
spin symmetric and hence couple to Suu = 1- Since these subconflgurations originate from the speciflc pentaquark 
state for which the four light quarks' spin has to add to a total spin S* = 0, we must also have S^^ = 1. This 
means that we have to identify the sd meson with K* (980) and the baryon with Ec (2445) with a total mass of 
ms^ + niK' = 3334 MeV, which is way above the threshold at 2906.8 MeV. In general, however, the pentaquark state 



'Because of the P-wave D* Dtv vertex this is compensated by a reduced effective coupling. 



68 



consi sts o f a superposition of various flavor assignments at the tips of the same color network. The right-hand side of 
Eq. ( FIO ) should therefore contain, in general, some weighted average so that 



m-p > a TOsrf (s=i) + rriccu (S=i) + P [msu + rricud] 
+ 1 [^ud + "nicsu] + 5 [m„a (5^1) + nicsd (s=i)] 



(Fll) 



with a-|-/3 + 7 + 5= l and a, (3,j,6> 0. The spins in the last configuration are fixed by the same consideration as 
those used in the case of sd — cuu separation discussed above. The corresponding total mass 

rup + ms^ = 3221.3 MeV 

is 315 MeV above the Dg — p threshold. In the other two configurations we will have in general admixtures of the 
lighter pseudoscalar tt and K states, and of the lighter Ac baryon. The lightest combination is 



= 2605.2 MeV. 



(The other alternative including the kaon. 



rriK + m\ 



2778.6 MeV 



is only 128 MeV below threshold.) 

We need an unlikely large admixture of the particular tt + Sc configuration to bring the right-hand side of Eq. (Fll) 
below 2906.8 MeV, the — p threshold. The baryon-meson inequalities are generally satisfied with a substantial 
margin of ~ 150 MeV. If this is also the case for the relation in Eq. (Fll), then it appears that the present inequalities 
cannot allow for a stable pentaquark bound state. 

It should be emphasized that the present approach, assuming that confinement via the specific mechanism of color 
flux tubes (or strings) is the dominant aspect, is practically orthogonal to the one-gluon exchange potential approach. 
Indeed in the latter, confinement is assumed from the outset by simply postulating that all the quarks in penta- and 
hexaquark states are inside the same universal bag. Yet confinement forces play a minimal role in the pentaquark 
binding. Rather, the overriding consideration is to maximize the overall hyperfine interaction even if that involves 
having large components of the wave function with qq pairs in a color sextet state. It is therefore logically consistent 
that the inequalities derived in the framework of one approach confiict with a prediction of a bound pentaquark in 
another framework. 

Nonetheless, the inequalities do suggest that if the pentaquark is not found, then neither QCD nor the naive quark 
models which so nicely explain the baryonic and mesonic spectra are at fault, but rather a lack of proper treatment of 
the basic confinement mechanism as adapted for these circumstances of pentaquark and/or hexaquark configurations. 

One other class of exotics are hybrid states containing an extra gluon such as {^piipj)sG and {ipiiljjipk)8G. The (• • •)8 
notation indicates that the quarks in the hadron couple to a color octet, which then forms the color singlet hadron 
with an additional gluon. To probe this sector, we should consider correlators of the form: 



If we represent the gluons via a field strength, then H{x,0) can be written schematically as 



(F12) 



H{x,0) - J dA*(A)F^,(x)i^^,(0)5^(x,0)5^(x,0) , (F13) 
and hence a Schwartz-type inequality implies 

\H\^ < f Ap{A)Fl{x)Fl^{Q) f dKA)S\{x,0)^S\{x,0) f dp{A) S\ix,0)^ S\{x,0) . (F14) 



This would imply, if the difficulty with the vacuum expectation value (0|i^^|0) could somehow be surmounted, that 



JTlhybrid < 2 ("^glucball + m^) ~ 0.8 GcV , 



(F15) 



which is a rather weak bound, from a phenomenological point of view. 
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A similar inequality in the context of SUSY models (in the mgquark oo limit )^ has proven useful in assessing the 
symmetry (phase) structure of that theory. Specifically, the possibility that massless, singlet, fermionic composites of 
qiQj and a gluino may manifest unbroken chiral symmetries, has been ruled out. Thi s wa s done by proving a QCD 



inequality between the mass of this state (the "hybridino") and the mass of the pion [121| 



ruhyhrid ^ m'-°\qiqjg) > mf^q^ ^ m-r: . (F16) 

Since unlike for gluons with nonlinear self-couplings one can define here the gluino propagator, the derivation of the last 
inequalities along lines paralleling those of Weingarten's in deriving rriN > (see Sec. is fairly straightforward. 
The hybridino-hybridino correlator (written in a concise, index-free form) can be bound by Schwartz inequalities: 



(•^hyb(^)-^hyb(0)) = y d^iS^S^Sx 

< f / d/. \S^A ( [ dM l^^l') ( / dM I^aI^ ) ■ (F17) 



The first and second factors in the last expression are the square root of a pion propagator. Hence 

(4yb(^)'^hyb(0)) < e-"-l-l , (F18) 



and Eq. (F16) follows 



APPENDIX G: QCD INEQUALITIES BETWEEN EM CORRECTIONS TO NUCLEAR SCATTERING 

In this Appendix we combine the results of Sec. ^ on the essentially positive nature of the EM contribution to the 
energies, and the approach of Sec. in order to apply "electromagnetic QCD inequalities" to scattering states. To 
this end we first elaborate on some of the material of Sec. [ij. 

We have seen in Sec. |l6| that QCD arguments leading to relations of the form 

B{aa) + B{bb) < 2B{ah) (Gl) 

between binding in channels (aa), (66), and [ah) are useful even if there are no bound states in those channels. Thus, 
imagine that the relative coordinate r — ri — r2 of the particle pair is confined to a sphere \r\ < i?, with R an infrared 
cutoff larger than the relevant scales in the problem. Let e° ; and e„^i denote the levels in various I channels without 
and with (respectively) the interparticle interactions turned on. Both sets of levels become dense as i? ^ c». [In 
particular, the free energies are e° ; = (fc° ;)^/2m with /cj^ ;i? = a;° ^ ~ {mr/2 + ln/2 -I- • • •), the nth zero of ji{x).] Yet 
a careful study of the shifts A(n, I) = — e„_; -I- e° ; reveals the full information on the phase shifts in various channels. 

The key observation is that the binding energy inequalities hold in general even when the R cutofi^ is placed, since 
the relevant QCD or nuclear interactions are short range. These imply then that An j{aa) + Anj{bb) < 2A„.i(a6). 
[For the unperturbed part we have trivially, from the definition of reduced masses, i{o.a) + e^^ i{bb) = 2e" i{ab).] 

Therefore the density of levels '^"j^"'''' in the (ab) channel is larger than the average of these densities in the (aa) and 
{bb) channels. 

Levinson's theorem implies that the phase shift Si"'''\k) serves as a "level counter" in the continuum limit. Hence 

dnijab) dSi{ab) 



^^From the discussion at the end of Sec. ^ a decoupling of the squarks and elimination of their Yukawa couplings is indeed 
required if any QCD-type inequalities are to be proven. 
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However, in the k ^ hmit the last quantity is the scattering length for I = 0. The inequality 

dSjaa) ^ dS{bb) (W(a6) 

yields in the fc — s- limit the desired inequality between scattering lengths 

a{aa) + a{bh) < 2a{ab) . (G4) 
We would next like to argue that the statement on the A/ = 2 energy shifts of NN continuum states, for example 

ABE{pp) + ABE{nn) < 2ABE{pn)i=i (G5) 
does imply the corresponding inequality for the respective scattering lengths 

a{pp) + a{nn) < 2a(pn) . (G6) 

The notation ABE{pn)j^i in Eq. ( |G5| ) refers to the shift from the idealized ABE{pn)j^i of some continuum level in 

the 7=1 (flavor symmetrized) np state for the case when aEM = and mi°^ = '^d'^- this exact /-spin limit, all 
{pn)i^i, (pp), and (nn) states would be the same, and hence so would be phase shifts, scattering lengths, etc. 

The AI = 2 combination 2ABE{pn)j=i — ABE{nn) — ABE{pp) and the similar combination of phase shifts, level 
densities, and scattering lengths are effected only by electromagnetism - i.e. by Oem 0- The general feature of 
positivity of such electromagnetic self energies can therefore be applied. This is so since the positivity applies not only 
for ground states but to any set of corresponding states which have equal charge density in the aEM = 0, m ^'^ = rn'^'^ 



limit. This then implies the positivity of all these combinations and in particular the desired relation (G6) 



In applying the suggested scattering length inequality ( G6 ) , we face a "technical" difflculty. The long-range Coulom 



bic interaction yields a (logarithmically) infinite Coulombic phase. In order to extract a meaningful result for a[pp), 
one must subtract this Coulomb phase shift. Our inequality is based precisely on the positive nature of the Coulombic 
self-interaction. Thus subtraction of even a part of this interaction could, in principle, invalidate the derivation of 
the inequality. Indeed, the very approach to continuum phase shifts by quantization in a finite sphere of radius R is 
jeopardized by the 1/r infinite range potential. 

We can, however, avoid the Coulombic pp phase shifts and still have a meaningful relation. This relies on the 
observation that screened versions of the Coulomb potential (i.e. Yukawa potential) still yield positive EM self- 
interactions [specifically the momentum space representation {fi^ + k'^)^^ is positive]. We can choose the cutoff fj,^^ to 
be of the order of the sum of the radii of the two nucleons. Thus the long-range Coulombic phase shift will be screened 
away, yet all other manifestations of EM interactions which occur during the nucleons' overlap will be maintained. 
These include the EM interactions between quarks in the two different nucleons or between their mesonic charge 
clouds, or more subtle indirect EM effects such as the 7r+ — tt" mass difference and the resulting different ranges for 



7r''",7r° exchange potentials |234] 



It is the latter effects which occur also in the np and nn systems, which we wish to consider. We cannot prove 
that the above cutoff procedure is indeed equivalent to Coulomb phase subtraction. However, any cutoff /i [or even 
arbitrary positive superpositions of potentials wit h diff erent cutoffs, such as J dfi^ a{^)e^^' /r] can be used. Hence for 



scattering lengths computed with this cutoff, Eq. (G6) holds 



a"^^^ (pp) + a"^^^ (nn) < 2a'"('') (pn) . (G7) 

We therefore believe that this inequality applies to the quoted "nuclear parts" . Recent measurements indicate that 
the suggested inequality is indeed satisfied wit h a wide margin. 



In principle, the original general relation (G3) contains much more information than that used to derive the 
scattering length inequality by taking I = and going to the fc — s- limit. Thus many other inequalities for other 
scattering parameters in all I waves can be derived for the two nucleon system, although these are not as useful and 
cannot be readily compared with data. 



Inequalities analagous to (G6) should hold for any 1 = 1/2 isospin multiplet. This in particular we should have 

a(3He^He) -I- a(3H3H) < 2a{^lle^K) . (G8) 

The 1 = 2 part can be extracted not only in scattering of two particles which are members of the same /-spin doublet. 
This suggests relations of the form 
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a(/H) + ain^Re) > aip^Re) + ain^Ti) 
a{K+n) + a{K°p) > a{K+p) + a{K°n) 

a{K°n) + a{K-p) > a{K°p) + a{K'n) , (G9) 

and many more. Since many of these relations involve unstable nuclear isotopes, testing them requires radioactive 
beam facilities. 

Finally we note that purely EM AI = 2 combinations of masses or scattering lengths can be found for nuclei of 
higher /-spin. 

APPENDIX H: QCD-LIKE INEQUALITIES IN ATOMIC, CHEMICAL, AND BIOLOGICAL CONTEXTS 

In this last Appendix we present various inequalities inspired by (more or less) justified analogies with the binding 
energy and correlator inequalities. 

1. Mass relations between compounds of different isotopes 

Different isotopes provide an almost ideal example of "flavor independent" interactions. Thus, let there be ni stable 
isotopes of a specific atom {Z'-^\a[^^) . . . {Z^^\A^n^) and n2 of another {Z^'^\A^^'^) . . . {Z'-'^\ aI?^). In the adiabatic 
(Born-Oppenheimer) approximation ]235( |, it is clear that the interatomic interactions are independent of the specific 

isotopes ) and chosen. The arguments in Sec. |l^ then imply that the binding energies Bij of the nin2 

possible compounds thus formed are a convex function of a reduced mass /i sampled at the values 

i{ZWA[^^)+miZ(^)A'-- 



with m(Z, A) the nuclear masses. This includes in particular the analog of the interflavor mass inequalities^ 
BE{x, x) + BE{x, y) < 2BE{x, y), e.g. 

BE{Il2) + BE(D2) < 2BE(BB) 
BEC'^02) + BEC^02) < 2BEC^02 O2) . (H2) 

We have not investigated the availability of data verifying these many true mass inequalities. 

The flavor (i.e. isotopic) independence arguments can be extended to more complex atomic compounds like 
XXZ, XZZ, etc. and the "convexity" relations {e.g. the analogs of those conjectured for baryons in QCD) are 
likely to hold. This will indeed be the case if the two- and three-body interactions between the nuclei satisfy the 



condition of positive exp{y} utilized in Lieb's proof in App. We are indebted to Phil Allen [236|, for pointing out 
to us this rather nice test of QCD-like inequalities. 

2. Conjectured inequalities for chemical bindings 

Atoms with closed shells or even a given {n,l) subshell constituting one Slater determinant of all possible mi, nig 
states are singlets of any relevant quantum numbers. This suggests using such states, say (NejZ — 10), as vacuum 
states upon which we can build "particle" states X — (Na|Z = 11) or "hole" states X — {F\Z = 9). Likewise another 
noble gas "vacuum", say (Ar|Z = 18), yields Y = {K\Z — 19), Y — {C\\Z — 17). If we use t his i dentification of 



vacuums, particles, and antiparticles, we might be tempted to conjecture, in "analogy" with Eq. (2.7), that 



BE{XX) + BE[YY) > BE{XY) + BE(YX) . (H3) 



^^The BE refers to the true molecular ground states or to sums over the first n states. The fact that Bose/Fermi statistics 
can imply even/odd Js in the rotation band of homoisotopic compounds has a negligible effect and can be corrected for. 
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This yields, for example, 



BE{NaF) + BE{KC\) > B£;(NaCl) + BE{KF) 

BE{MgO) + BEiCaS) > BE(MgS) + BE{GaO) , (H4) 

etc. We might, however, consider half-filled shells, e.g. C, Si, etc . with an equal number of electrons and holes to be 
the correct atomic analog of the vacuum, in which case Eq. ( ]H3| ) translates into 

BE{UY) + BE{N£iC\) > BE{UC\) + BEQ^aF) , (H5) 

etc. 

3. Mass inequalities in nuclear physics 



For a while we were excited about the prospect |237] that true QCD inequalities could show that even-even N — Z 



states are - barring Coulomb effects - the most tigh tly bound 



The sophisticated third order Garvey-Kelson |238| difference relations connect masses of many isotopes. However, 
there is no obvious pattern of deviations from the relations, and we have not found a simple motivation for any such 
pattern. 

4. A biological analog for inequalities between correlators 

The pseudoscalar mass inequalities of Sec. |l^ reflect simple Schwartz inequalities for correlators, i.e. for weighted 
bilinears in quark propagators. Generically the latter have gauge interaction induced "phases" . These cancel in the 
tr{SiSl) combinations appearing in the particular case of the pseudoscalar propagators]^ 

In an extremely wide variety of circumstances we may encounter joint propagation of two equal or two different 
entities. Quantum phases are typically irrelevant and we can assign a positive probability for the propagation of 
from "Pi" to "P2" for any given set of relevant influencing "factors" An - the analog of the background gauge fields 
in the QCD case: 

r M(at Pi) ^ ^(at P2)} . (H6) 

The overall propability of A "propagating" from Pi to P2 is then given by a "functional" (path integrated) averaging 
over the distribution of the {An} factors: 

V {Aiat 1) ^ A{at 2)} = J d^i{Ai... An} V {Ai&t 1) ^ yl(at 2)} , (H7) 

with a normalized / d/i {A} = 1 positive measure of {Ai . . . An}. 

The probability of the joint propagation of A{1),B{1) A{2), B{2) is given by the corresponding weighted average 
of the bilinear product of propagators: 

^ { ^(1) Z e(2) } = / ^ {"^(1) ^ "^(2)} |{A„} • V {Bil) ^ Bi2)} . (H8) 

Likewise, the probability of joint propagation of two A objects from (1) to (2) (or two B objects) is given by a similar 
expression involving squares of propagators: 



V 

and 



{1(1)^;4(2)}^ |<^^^{M['P{Ail)^A{2)}\^A„^]\ (H9) 



^'^Schwartz inequalities have often been used in other contexts of particle physics. See, e.g., Ref. 
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6(1) 
6(1) 



^6(2)}= [^{^W 



■m}\{A„}] 



The Schwartz inequaUty readily implies the final desired relation: 



A{1) 
A{1) 



A{2) 
A{2) 



6(1) 
6(1) 



6(2) 
6(2) 



> 



Ail) 
6(1) 



A{2) 
6(2) 



(HIO) 



(Hll) 



There are certainly innumerable applications of the above relation in all areas of science technology and life sciences 
(most of which were very likely well known for some time). In the last page of this review we consider one particular, 
somewhat intriguing, and exotic application concerning the sex of non-identical twins. 

The conception of non-identical twins can be viewed as the joint propagation in the same (or relatively similar 
"background" of two sperms. Hence we expect that for non-identical twins p| 



7'(male, male )7'( female, female) > [7^(female, male)] 



(H12) 



Clearly the probability of any specific single (or double) conception need not be accurately reflected in the percentages 
at birth. We could however still derive the same inequalities if we consider not merely the propagation from "inception" 
to "conception" but also the subsequent nine month long "timelike" propagation to birth. 

The extent to which the inequality ( H12| ) is satisfied, i.e. the deviation of the ratio of the RHS and LHS from unity, 
can serve as a measure of the overall degree of non-parallelism of the male and female birth "vectors" , i. e. as a crude 
measure of the total importance of "variables" known, or as yet unknown, in determining the sex of the fetus, which 
is a quantit y of considerable interest. 

A priori [240| there could be some unique rather surprising effect which would disqualify our proof and possibly 
reverse the sign in Eq. (H12). This would be the case if there were a direct interaction between the propagating 
elements - in this case the fetuses of the two twins. Thus following the biblical story of the rather unrestful pregnancy 



of Rebekah |241| with Esau and Jacob (a clear case of non-identical twins!), we can assume that the competition 
between equal sex brothers (and possibly sisters) extends to the prenatal stages. Clearly if it is too strong it could 
conceivably disrupt (MM) and (FF) births, thus tending to reverse the sign of the effect considered. Hopefully this is 
not the case. 

We have not attempted to verify any of the above non-QCD mass inequalities. In particular, the last twin inequalities 
may require rather extensive statistical analysis .0 We hope that the inequalities will eventually be tested. 
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FIG. 1. Vertex and propagator insertion of the one gluon exchange diagram with the Ai • A2 structure intact. 



FIG. 2. The general two quark interaction is a sum of Ai ■ A2 (octet exchange) and li • I2 (singlet exchange) interactions. 



FIG. 3. A non-separable three body interaction in the baryon due to the three gluon vertex which is shown, however, to 
vanish. 



FIG. 4. (a) The color string picture for a meson, (b) The color network for a qqq baryon with one junction. 



FIG. 5. Illustration of how a given Y configuration in a baryonic wave functional yields configurations for three mesonic trial 
wave functionals. The relation 2//123 = H12 + H23 + -ffai is illustrated by the fact that the term depicted by the small 
string distortion, the tpDip terms depicted by the motion of the end quarks, and the E (corresponding to the weighted string 
length) all appear twice. 



FIG. 6. The triangular inequalities and relevant vectors for the comparison of the potential energies in the strong coupling 
limit of the baryonic Y configuration and mesonic subsystems. 



FIG. 7. Some matrix elements ( configuration' |Hqcd| configuration ) on a small 3x3 dimensional lattice. The upper entry 
in each box is the matrix element with /i ~ 1/a the lattice energy scale. The lower entry indicates the terms in the Hamiltonian 
contributing to this particular matrix element. 
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FIG. 8. Configurations in the generalized baryonic wave functional which still allow separation into three mesonic subsystems. 



FIG. 9. Flavor connected (a) and flavor disconnected (b) contraction contributing to two-point correlation functions of quark 
bilinears. 



FIG. 10. The focusing effect of a parallel B field on charged particles propagating from x to y and its limited effect on 
propagation between regions of size A once B > 



FIG. 11. (a) The unique contraction when all fiavor indices are distinct, (b) An alternate contraction when we have 
permutable quarks in the two currents. 



FIG. 12. The planar duality diagrams representing the euclidean correlations Fi,F2, and Fs. For F2 we use only the specific 
contraction with uu(ss) exchanged in the t{T) channel respectively. 



FIG. 13. The evolution of an initial ud 1 state on different time scales. 



FIG. 14. The original Wilson loop (W), its parts Wi and W2, and the reflection paths used in proving Eq. (18.3) 



FIG. 15. The unique "8" pattern of contraction relevant to the {K\B\n) > (A'|yl|7r) inequality. 



FIG. 16. The "eye" contraction relevant to the (A'^jAlO) matrix element. 



FIG. 17. Illustration of how as r ^ and — > u (and the diamond- like configuration degenerates into a vertical line), we 
tend to get, due to the smoothness of the configuration, similar propagators, and hence a monotonic decrease of V'^d?^) 
with \r\. 



FIG. 18. The color string network for a quadriquark with two light quarks and two heavy quarks. 
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FIG. 19. The color string network for a quadriquark (solid line) plus its doubled network with flux lines reversed (dashed 
line). 



FIG. 20. (a) The color network for the putative pentaquark state with three junction points (1), (2), (3). (b) The string 
picture of the trial baryon + meson states obtained from that of the pentaquark by omitting the (2) - (3) string bit, and 
reversing the flux direction in the (2) - de and (3) - Q/ sections. 



FIG. 21. The three polygonal paths Xi,X2,X3 consisting of the overall periodic propagation of the particles at Xi,X2,X3 
from "t" = /? = to "i" = /3 with Xi{f3) = Xi{0). 
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